Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,3020TG46.2% 134.2 / 18.6 80H219P (CAC→CCC) yghSputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (4/39);  new base G (37/0);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.19e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133153‑3133450
                                                                                                                                                     |                                                                                                                                                     
aaTGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                                                                                                      <  2:626330/151‑1 (MQ=255)
 aTGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtgt                                                                                                                                                     <  2:704630/151‑1 (MQ=255)
     aGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGcc                                                                                                                                                 <  1:887264/151‑1 (MQ=255)
     aGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGc                                                                                                                                                  >  1:1148049/1‑150 (MQ=255)
             tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCCTCCCGt                                                                                                                                         >  1:1048233/1‑151 (MQ=255)
                ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGtt                                                                                                                                        >  2:570367/1‑149 (MQ=255)
                   gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAgg                                                                                                                                   <  1:160581/151‑1 (MQ=255)
                   gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAgg                                                                                                                                   <  2:1126448/151‑1 (MQ=255)
                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGAt                                                                                                                                 <  1:315377/151‑1 (MQ=255)
                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGAt                                                                                                                                 <  1:977793/151‑1 (MQ=255)
                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGAt                                                                                                                                 <  2:390290/151‑1 (MQ=255)
                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTcca                                                                                                                                     >  1:472884/1‑145 (MQ=255)
                           aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTg                                                                                                                            <  2:1030938/150‑1 (MQ=255)
                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAGGATCTTTTCGcc                                                                                                                        >  2:134253/1‑151 (MQ=255)
                                tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCa                                                                                                                       <  1:852856/150‑1 (MQ=255)
                                tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGAGTTTTTCGCCa                                                                                                                       >  1:597319/1‑150 (MQ=255)
                                 aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTTCGCCCtt                                                                                                                     >  1:992765/1‑151 (MQ=255)
                                  aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTGCGCCCtt                                                                                                                     >  2:637636/1‑150 (MQ=255)
                                  aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGTCCATCCAGTTCAAGGATTTTTTCGCCAtt                                                                                                                     >  1:334556/1‑150 (MQ=255)
                                   cAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGATCTTTTCGcc                                                                                                                        >  1:1527/1‑146 (MQ=255)
                                     gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTTCAGGAATTTTTTGCCAtttta                                                                                                                  >  1:799826/1‑147 (MQ=255)
                                      tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTTCAGGGTCTTTTCGCCAtttgat                                                                                                                 >  1:245528/1‑146 (MQ=255)
                                       ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg                                                                                                                <  2:575714/150‑1 (MQ=255)
                                       ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg                                                                                                                <  1:1063049/150‑1 (MQ=255)
                                        tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGt                                                                                                               <  2:884308/150‑1 (MQ=255)
                                        tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGGTTACAGATCTTTGCGCCCTTgga                                                                                                                  >  1:214878/1‑145 (MQ=255)
                                        tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGGTCAAAGGTCTTTTCGCCCTTGaa                                                                                                                  >  2:75659/1‑147 (MQ=255)
                                          agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAg                                                                                                            <  1:1020721/151‑1 (MQ=255)
                                              tggaGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCTCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTg                                                                                                         <  2:721965/147‑1 (MQ=255)
                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgc                                                                                                        <  2:71422/150‑1 (MQ=255)
                                                     cATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                                                 <  2:1053223/151‑1 (MQ=255)
                                                     cATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGGTgcggcg                                                                                                    >  2:443722/1‑148 (MQ=255)
                                                      aTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTTAAGGATCTTTTCGCCATTTAATGTCAGTTGCGGCGCaa                                                                                                 >  1:242970/1‑150 (MQ=255)
                                                       tCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGaa                                                                                                 >  2:171315/1‑149 (MQ=255)
                                                          aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCg                                                                                            >  2:1147173/1‑151 (MQ=255)
                                                           ttCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAAAGGCgg                                                                                           >  2:1139422/1‑150 (MQ=255)
                                                            tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAACGGCggt                                                                                          >  2:19340/1‑149 (MQ=255)
                                                              gTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttt                                                                                        <  2:574414/151‑1 (MQ=255)
                                                                 aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAtttttt                                                                                     <  1:1040283/151‑1 (MQ=255)
                                                                  aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTc                                                                                    <  2:627269/151‑1 (MQ=255)
                                                                     catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTg                                                                                 <  2:870240/151‑1 (MQ=255)
                                                                     catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCCCCCCTCCAGCTCTAAGGTCTTTTCGCCCTTGAAAGTGAGTTGCGGGCTAACAGCGCttttttt                                                                                    >  2:413984/1‑147 (MQ=255)
                                                                      atcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGGGTAACGGGGATTTTTTTCTGt                                                                                >  2:1083027/1‑151 (MQ=255)
                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCGGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACg                                                                            <  2:1074087/151‑1 (MQ=255)
                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACg                                                                            <  1:688722/151‑1 (MQ=255)
                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACg                                                                            <  2:540763/151‑1 (MQ=255)
                                                                            tATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCACCCGTTCAAGGATCTTTGCGCCCTTGTATGGCAGGTTCGGGGTAAGGGGGAGTTTTTTCTGtga                                                                              >  2:663203/1‑145 (MQ=255)
                                                                             aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTg                                                                         <  2:1033253/151‑1 (MQ=255)
                                                                             aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGGTCTTTGCGCCCTTGAATGTCAGTTGCGGGGGAACGGCGATTTTTTTCTTTAAAGCt                                                                          >  1:1163545/1‑150 (MQ=255)
                                                                               ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGTTCAAAGATCTTTTCGCCATTGTATGTCAGTTTCGGCGGAACGGCGGTTTTTTTCTTTTAAGCTG‑cc                                                                      >  1:110432/1‑151 (MQ=255)
                                                                                   tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAAGTCCAGGAACTTTTCGCCCTTTTATGTCCGTTGCGGCGTAACCGCGGTTTTTTTCttttaatccgcccgtt                                                                   >  2:371363/1‑136 (MQ=255)
                                                                                       gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAATTCAAGGGACTTTGTGCCACTTCATGTCGGTTGCGGCGCAACGGCGAATTTTTTCTCTaaaaccgccgggtggt                                                                >  1:391176/1‑136 (MQ=255)
                                                                                        tGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGGTTTTTTTCTGTAACGCTG‑CCCGTtggggg                                                             >  2:1142748/1‑151 (MQ=255)
                                                                                          aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCCAAGATCTTTTTGCCATTTTATGGCAGTTGCGGGGTAAAGGCGATTTTTTTCTTTAACGCCG‑CCAGTtggggg                                                             >  1:800056/1‑149 (MQ=255)
                                                                                           tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTTGCCATTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCTG‑CCAGTtgg                                                                >  2:561331/1‑145 (MQ=255)
                                                                                            cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAg                                                         <  2:321627/151‑1 (MQ=255)
                                                                                              ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCCCCCATCCAGTTCCAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTCCTTTAAAGCTG‑CCCGGTGGGGGGCa                                                          >  2:582521/1‑148 (MQ=255)
                                                                                                aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGttt                                                     <  1:475147/151‑1 (MQ=255)
                                                                                                 aaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCGAGGGGGTTTGCGCCATTTAATGTCCGTTGGGGCGGTACGGCGCTTTTTTTCTTTTACCCCG‑CCCGTTGGTGGGCacggt                                                      >  1:518701/1‑145 (MQ=255)
                                                                                                         cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGAACTTTGCGCCCTTTAATGTCAGGTGCGGCGTAACGGCGATTTTTTTCTGTAAAGCTG‑CCAGTTGGGGGGCTGGGTTTCGCGCaa                                             >  2:1013105/1‑150 (MQ=255)
                                                                                                          gTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTGGCCCCTCCCGTTCTAGGATCTTTTTGCCATTGAATGTTAGGTTCGGGGTAACGGCGCTTTTTTTCTGTTAAGCCG‑CCAGTTGGTGGTGAGGGTTTTGCGCcaa                                            >  1:619126/1‑150 (MQ=255)
                                                                                                          gTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAAc                                           <  2:245528/151‑1 (MQ=255)
                                                                                                          gTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAAc                                           <  1:1187490/151‑1 (MQ=255)
                                                                                                           tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGGTCTTTGCGCCATTTAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCCG‑CCAGTTGGGGGGCAGGTTTTCGCGCaaa                                            >  1:342048/1‑149 (MQ=255)
                                                                                                           tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTTCGGCGGAACGGGGATTTTTTCCTGTAAAGCTG‑CCAGTTGGGGGGCGGGTTTTCGCGCaaa                                            >  1:103916/1‑149 (MQ=255)
                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAAcg                                          <  1:40559/150‑1 (MQ=255)
                                                                                                                 gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGTGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGttt                                    <  2:877703/151‑1 (MQ=255)
                                                                                                                 gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTTGGGCGGAACGGCGGTTTTTTTCTTTAACGCCG‑CCAGATGGGGGGCAGGTTTTCGGGCAAACGCGGgtt                                    >  1:847525/1‑151 (MQ=255)
                                                                                                                      tAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTTCAGGATCTTTTCGCCCTTGAATGTCAGGTGGGGGGTAACGGGGGTTTTTTTCTTTTACGCCGCCCGGTGGGGGGGCGGGTTTA‑GGGCAAAAGCGGgttgtt                                 >  1:272815/1‑149 (MQ=255)
                                                                                                                       aGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtca                               <  1:1026929/150‑1 (MQ=255)
                                                                                                                          gACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  1:344359/151‑1 (MQ=255)
                                                                                                                            cGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATg                         <  1:367302/151‑1 (MQ=255)
                                                                                                                             gCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGGTCTTTTCGCCATTGAATGTCAGTTGCGGGGGAACGGCGATTTTTTTCTTTTACGCTT‑CCCGTTGGGGGGCAGGGTTTCGCGCGAACGCGGGTTGTTCTTCaaag                         >  1:807749/1‑148 (MQ=255)
                                                                                                                               ttGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGCTTAAGGATCCTTTTGCCATTGAATGTCAGTTGGGGCGTAACGGCGGTTTTTTTTTTTTACGCCG‑CCCGTTGGTGGGCCGGGTTACGCGCAAAAGCGGGTTTTTCTTCAATTcc                       >  1:434822/1‑150 (MQ=255)
                                                                                                                                 gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                    <  1:570365/151‑1 (MQ=255)
                                                                                                                                    gAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACg                 <  1:1164432/151‑1 (MQ=255)
                                                                                                                                      aTTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGaa               <  1:7439/151‑1 (MQ=255)
                                                                                                                                         tCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATc             <  2:1048240/150‑1 (MQ=255)
                                                                                                                                             cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca         <  1:1109768/150‑1 (MQ=255)
                                                                                                                                               gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATa      <  1:141550/151‑1 (MQ=255)
                                                                                                                                                cggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     >  2:775600/1‑151 (MQ=255)
                                                                                                                                                   ggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  <  2:334557/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133153‑3133450

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: ‑