Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,3020TG42.4% 108.4 / 16.3 66H219P (CAC→CCC) yghSputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (3/35);  new base G (28/0);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.32e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACG  >  NC_000913/3133156‑3133448
                                                                                                                                                  |                                                                                                                                                    
gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCg                                                                                                                                                  <  2:361224/151‑1 (MQ=255)
   gTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgggttcac                                                                                                                                                >  2:943260/1‑146 (MQ=255)
        gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCa                                                                                                                                        <  2:1102420/151‑1 (MQ=255)
           gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGtt                                                                                                                                     <  2:81716/151‑1 (MQ=255)
            cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGtt                                                                                                                                     <  2:1071494/150‑1 (MQ=255)
             ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCa                                                                                                                                   <  2:523909/151‑1 (MQ=255)
             ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCa                                                                                                                                   <  1:222069/151‑1 (MQ=255)
                gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAgg                                                                                                                                <  2:131150/151‑1 (MQ=255)
                  ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGa                                                                                                                               <  1:202541/150‑1 (MQ=255)
                  ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGa                                                                                                                               >  2:1096949/1‑150 (MQ=255)
                  ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑ccccccccgtttaa                                                                                                                                   >  2:780830/1‑136 (MQ=255)
                    ttGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGG‑CCCAT‑CCAGTTTAAGGGTCt                                                                                                                            >  1:1185032/1‑151 (MQ=255)
                      gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCttt                                                                                                                          <  2:250978/151‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCt                                                                                                                            <  2:177378/144‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTgcgc                                                                                                                      <  1:1239423/150‑1 (MQ=255)
                           aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCCT‑CCCGTTCCAAGGTCTTTTCGcc                                                                                                                     >  2:60871/1‑151 (MQ=255)
                            ttAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCC‑CCAGGTCCAGGATCtttttggc                                                                                                                     >  1:1089819/1‑145 (MQ=255)
                              aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTG‑CCCAC‑CCAGGTCAAAGAACTTTgcggcctt                                                                                                                  >  1:445847/1‑146 (MQ=255)
                              aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCCT‑CCCGTTCCAGGGTCTTTTCGCCAtt                                                                                                                  >  1:670745/1‑151 (MQ=255)
                               aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                 <  2:353856/151‑1 (MQ=255)
                                      gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTc                                                                                                           <  1:1158570/150‑1 (MQ=255)
                                      gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTG‑CCCCTCCCCGTTCAA‑GAACTTTTCGCCCtt                                                                                                                  >  1:1008168/1‑143 (MQ=255)
                                          ttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGtt                                                                                                       <  1:539364/150‑1 (MQ=255)
                                           tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCCT‑CCCGTTTCAGGATCTTTTCGCCCTTGAATGTCCGtt                                                                                                       >  1:343038/1‑149 (MQ=255)
                                            gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCG‑CCCCT‑CCCGTTCAAGGAACTTTTCGCCCTTGAATGGCAggt                                                                                                       >  2:985903/1‑146 (MQ=255)
                                              aGTGCATCCATTCGTTAAATCATCCCCTAACCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcgg                                                                                                   <  2:692207/150‑1 (MQ=255)
                                                tGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCAC‑CCAGTTCAAGGATCTTTTTGCCATTGTATGTTAggttgggcg                                                                                                 >  1:1038310/1‑142 (MQ=255)
                                                 gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTa                                                                                               <  1:705713/151‑1 (MQ=255)
                                                   aTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCC‑CCAGGTCAAAGATCTTTTCGCCCTTGAATGTCAATTGCGGCGGaa                                                                                              >  2:317925/1‑150 (MQ=255)
                                                   aTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCAT‑CCAGGTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTaa                                                                                              >  2:514553/1‑150 (MQ=255)
                                                        ttCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGG‑GCCCC‑CCCGTTCAAGGAACTTTTTGTCCTTTTATGTCAGTTTCGGCGGGACggg                                                                                          >  1:282915/1‑148 (MQ=255)
                                                        ttCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCCT‑CCCGGTCCAGGATCTTTTTGCCCTTTAATGTTAGTTGCGGGGTAACGGCg                                                                                         >  1:334203/1‑150 (MQ=255)
                                                               aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGG‑GCCCT‑CCAGTTCCAGGATCTTTGCGCCCTTGTATGTCAGTTGCGGGGGAACGGCGAttttttt                                                                                 >  1:121539/1‑150 (MQ=255)
                                                                 tcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCCC‑CCGGTTCAA‑GATCTTTTCCCCCTTTAATTTCAGTTTCGGGGGAACGGCGGttttttt                                                                                 >  1:1188118/1‑148 (MQ=255)
                                                                  catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCCGTTTAAAGAGCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAACGGCGGttttttt                                                                                 >  1:577104/1‑147 (MQ=255)
                                                                  catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCCGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAAAGGCGGtttttt                                                                                  >  2:1198235/1‑147 (MQ=255)
                                                                     ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTaa                                                                           <  1:248372/151‑1 (MQ=255)
                                                                      cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAAc                                                                          <  1:462171/151‑1 (MQ=255)
                                                                              aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAg                                                                  <  2:368677/151‑1 (MQ=255)
                                                                              aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAg                                                                  <  1:829610/151‑1 (MQ=255)
                                                                              aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCCGTTCAAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTTTAACGCTGcc                                                                    >  2:398186/1‑149 (MQ=255)
                                                                              aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCCT‑CCACATACCGGATATTTTTGCCACTGAAAGTCCGTTGTGGCGGAACGGGGCGTTTTTTCTGTTAAGCTGc                                                                     >  1:955273/1‑148 (MQ=255)
                                                                                tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGt                                                                 <  2:973109/150‑1 (MQ=255)
                                                                                   ggTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTG‑GCCAT‑TCCGATCAAAGAGTTTTGGGGCCTTTTATGGCAGTTTTGGCGGAAAGGGGATTTTTTTCTGttcagacggccact                                                                >  2:1142055/1‑135 (MQ=255)
                                                                                      gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTggtggt                                                          <  2:1087683/151‑1 (MQ=255)
                                                                                        tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCG‑CCCCC‑CCCGTTTCAGGGTCTTTGTGGCATTTAATTATCGTTGCGGCGTAACAGGGGGTTTTTTTTTTAACGCCGACCGTtggtg                                                            >  1:295388/1‑147 (MQ=255)
                                                                                         cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCa                                                        <  2:948626/150‑1 (MQ=255)
                                                                                         cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCAGTTTCAGGATCTTTTTGCCCTTTTATGTCAGGTGGGGGGGGAAGGCGGATTTTTTTTGTTACACTGCCCGGtggtgg                                                           >  1:837839/1‑147 (MQ=255)
                                                                                           ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGt                                                     <  1:367928/151‑1 (MQ=255)
                                                                                              aaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCCGCTTAAAGATATTTTCGCCCTTTAATGTCCGTTGGGGCGTAAAGGGGGTTTTTTTCTGTAAAGGCGCCCGGtggtgggtcgggtt                                                   >  2:818104/1‑142 (MQ=255)
                                                                                                 cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTG‑CCCCT‑CCAGTTTAAAGATCTTTTCGCCCTTTGATGTACCGTGCGGGGGAACGGCGGTTTTTTTCTTTAAAGCTGCCAGtttggggggtggttt                                                   >  2:519561/1‑134 (MQ=255)
                                                                                                    ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGG‑CCCCT‑CCAGTTCAAGGGTCTTTTCGCCCTTGGATGTTAGGTGGGGGGTGACGGCGGTTTTTTTCTTTTAAGCTGCCCGtttggggggggggttt                                                  >  1:524095/1‑131 (MQ=255)
                                                                                                         gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAAcgc                                       <  1:984930/151‑1 (MQ=255)
                                                                                                          aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCCGTTCAAAGAACTTTTCGTCACTTTATGTTAATTGCGGGGCAAAGGCGATTTTTTCCTGTTAACATGCCAGTTGGGGGGCGGGGTTatgg                                               >  1:1162749/1‑139 (MQ=255)
                                                                                                           aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCgg                                     <  2:1011694/151‑1 (MQ=255)
                                                                                                              gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCG‑CCCCT‑CCCGGTCAAAGATCTTTTCGCCCTTGAATGTTAGTTGCGGCGGAACGGCGATTTTTTTCTGTAAAGCTGCCCGTTGGGGGGCAGGGTTTCGCGCAAAAGCGGgtt                                  >  2:149344/1‑151 (MQ=255)
                                                                                                               cgcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGAGCTTTTCGCCATTTAATTTTAGATGCGGCGGAAAGGCGGTTTTTTTTTGTAACGCTTCCCGTTGGGGGGCAGGGTTACGCGCCAACGCGGgtt                                  >  2:968098/1‑150 (MQ=255)
                                                                                                                gcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGt                                <  2:645324/151‑1 (MQ=255)
                                                                                                                       gACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtcatca                          <  2:1245011/150‑1 (MQ=255)
                                                                                                                        aCGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAAt                        <  1:1077739/151‑1 (MQ=255)
                                                                                                                         cGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATg                       <  2:401744/151‑1 (MQ=255)
                                                                                                                         cGCTTGCAGAATTTCATCGGCGGGGGGGGG‑CCCCT‑CCAGTTCTAAGATCTTTGCGCCATTGAAAGTCAGTTGGGGCGTAACGGCGGTTTTTTTCTGTAAAGCGTCCCGTTGGGGGGCGGGGTTTAGCGCGAACACGGGTTTTttatcca                         >  1:1198686/1‑147 (MQ=255)
                                                                                                                              gCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                  <  1:686602/151‑1 (MQ=255)
                                                                                                                               cAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGa                 <  2:83532/151‑1 (MQ=255)
                                                                                                                                aGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTATGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                <  1:1129321/151‑1 (MQ=255)
                                                                                                                                          cggcggGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca       <  1:77549/150‑1 (MQ=255)
                                                                                                                                              ggggTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACg  <  1:60868/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCG‑CCCAT‑CCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACG  >  NC_000913/3133156‑3133448

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑