Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,3020TG41.7% 182.8 / 18.8 97H219P (CAC→CCC) yghSputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (10/46);  new base G (40/0);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 7.58e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133153‑3133451
                                                                                                                                                     |                                                                                                                                                     
aaTGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                                                                                                      <  1:1041126/151‑1 (MQ=255)
   gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgggtgg                                                                                                                                                   >  2:175971/1‑149 (MQ=255)
   gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCg                                                                                                                                                   >  1:1117640/1‑151 (MQ=255)
     aGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGcc                                                                                                                                                 >  1:984392/1‑151 (MQ=255)
      gTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGccc                                                                                                                                                >  1:31554/1‑151 (MQ=255)
         gggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCAt                                                                                                                                              <  1:1459749/150‑1 (MQ=255)
         gggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCAt                                                                                                                                              <  1:1371558/150‑1 (MQ=255)
         gggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCAt                                                                                                                                              <  1:52710/150‑1 (MQ=255)
         gggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATc                                                                                                                                             >  2:4618/1‑151 (MQ=255)
           gCTGCGGCGCGGTTGTACTATTCACAGTTTGAGTTTGTAGTGCTTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGCATTGCGCGCAGTGACGCTTGTAGAATTTCATCGGCGGGGGGGCGACCAttc                                                                                                                                            >  2:14209/1‑148 (MQ=255)
           gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTccc                                                                                                                                           >  1:942974/1‑150 (MQ=255)
             tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGt                                                                                                                                         >  1:544746/1‑151 (MQ=255)
                ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCa                                                                                                                                      <  1:1203770/151‑1 (MQ=255)
                ggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTAGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCa                                                                                                                                      <  1:1082113/151‑1 (MQ=255)
                    cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGAGTCGCGCCTCCAGATcttgg                                                                                                                                   >  2:387059/1‑146 (MQ=255)
                     ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTTaaag                                                                                                                                   >  2:1222384/1‑147 (MQ=255)
                       ttGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCt                                                                                                                               >  2:254623/1‑151 (MQ=255)
                         gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCttt                                                                                                                             >  2:270278/1‑151 (MQ=255)
                           aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgc                                                                                                                           <  1:1017952/151‑1 (MQ=255)
                           aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTg                                                                                                                            <  2:569934/150‑1 (MQ=255)
                           aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTg                                                                                                                            <  2:1429747/150‑1 (MQ=255)
                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGTTCTTTCTCCCCTCCAGTTTAAGGaaaatttaccc                                                                                                                        >  1:372664/1‑141 (MQ=255)
                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTGCGcc                                                                                                                        >  1:1433027/1‑151 (MQ=255)
                                 aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGATCTTTGCGGCAtt                                                                                                                     >  2:544344/1‑151 (MQ=255)
                                  aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                    <  2:1146219/151‑1 (MQ=255)
                                     gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGaag                                                                                                                 <  2:544383/151‑2 (MQ=255)
                                      tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCCAGGATCTTTTCGCCATTGAATg                                                                                                                >  2:33749/1‑151 (MQ=255)
                                      tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCCTTGTATg                                                                                                                >  1:392497/1‑151 (MQ=255)
                                          agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTTAAGGATCTTTTCGCCCTTGAATGTc                                                                                                              >  1:1122660/1‑149 (MQ=255)
                                             ttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTg                                                                                                         >  2:643673/1‑151 (MQ=255)
                                              tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTg                                                                                                         <  2:552944/150‑1 (MQ=255)
                                              tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgc                                                                                                        >  2:1466368/1‑151 (MQ=255)
                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTgcg                                                                                                       >  2:469738/1‑151 (MQ=255)
                                                   tGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCGCCAGGTCAAGGCTGTTTGCGCCCTTGAATGGCCATTgcggcgg                                                                                                   >  2:1367371/1‑150 (MQ=255)
                                                    gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTa                                                                                                  <  2:821025/151‑1 (MQ=255)
                                                    gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAGGAACTTTGCGCCATTGAATGGCAGTTGCGGCGt                                                                                                   >  1:1074895/1‑150 (MQ=255)
                                                     cATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                                                 <  2:1439593/151‑1 (MQ=255)
                                                        ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                              <  1:1295270/151‑1 (MQ=255)
                                                         cATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGc                                                                                             <  2:513406/151‑1 (MQ=255)
                                                         cATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCCCTTTAAGGATCTTTTCGCCATTTTAAGTCAGTTGCGGCGCAAcagg                                                                                             >  2:247760/1‑148 (MQ=255)
                                                          aTTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCCGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCg                                                                                            >  1:1152283/1‑151 (MQ=255)
                                                             cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAt                                                                                          <  1:459008/150‑1 (MQ=255)
                                                                 aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttttt                                                                                      <  2:14952/150‑1 (MQ=255)
                                                                 aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttttt                                                                                      <  1:1288476/150‑1 (MQ=255)
                                                                   atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTTGCCCATCCAGGTCAAGGATCTTTGCGCCATTGAATGTTAGATGTGGGGTAACGGCGATTTTTTcc                                                                                   >  1:1129542/1‑151 (MQ=255)
                                                                     catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTg                                                                                 <  1:508063/151‑1 (MQ=255)
                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACg                                                                            <  1:756355/151‑1 (MQ=255)
                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACg                                                                            <  1:535815/151‑1 (MQ=255)
                                                                           cTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGATCAAGGAACTTTGCGCCCTTGAATGTCAGTTGGGGCGGAAAGGCGGTTTTTTTCCGTAAAGc                                                                           >  2:235447/1‑151 (MQ=255)
                                                                                 aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAg                                                                     <  2:830049/151‑1 (MQ=255)
                                                                                 aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGGTCTTTGCGCCATTGAATGTCAGTTGGGGCGGAACGGCGAGTTTTTTCTTTAACGCTGcccg                                                                     >  2:1166149/1‑149 (MQ=255)
                                                                                    ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGtt                                                                   <  2:991364/150‑1 (MQ=255)
                                                                                     tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtg                                                                  <  1:206786/150‑1 (MQ=255)
                                                                                     tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGCGCCATTGAATGTCCGTTGCGGGGTAACGGCGATTTTTTTCTGTAACGCTGCCAGTtgg                                                                 >  1:915863/1‑151 (MQ=255)
                                                                                      ggTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtggg                                                                >  1:682395/1‑150 (MQ=255)
                                                                                            cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATATTTGCGCCCTTGAATGTCCGTTGCGGCGGAACGGCGGTTTTTTTCTGTAAAGCTGCCCGTtggggggc                                                            >  2:392082/1‑147 (MQ=255)
                                                                                              ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCCTTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCAGTTGGTGGTCAggg                                                        >  1:1179633/1‑150 (MQ=255)
                                                                                                aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGttt                                                      <  1:910287/151‑1 (MQ=255)
                                                                                                aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGttt                                                      <  2:1122059/151‑1 (MQ=255)
                                                                                                   gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAGGTCAAGGATCTTTTGGCCCTTTAATGTCGATTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCAGTTTGGGGGCAGGTTTa                                                     >  1:1360791/1‑149 (MQ=255)
                                                                                                      cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAAAGGCGATTTTTTTCTGTAAAGCTGCCAGCTGGGGGGCAGGGTTACgcgc                                                >  1:641437/1‑151 (MQ=255)
                                                                                                       ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCa                                               <  1:1332062/151‑1 (MQ=255)
                                                                                                       ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCAGTTGGGGGGCAGGTTTACgcgc                                                >  1:245042/1‑150 (MQ=255)
                                                                                                           tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTCAAAGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGTACGGCGATTTTTTTCTGTAACGCCGCCAGTTTGTGGGGCGGTTTACGCGCaaaac                                           >  2:316101/1‑149 (MQ=255)
                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCCGTTGCGGCGTAACGGCGATTTTTTTCCGTAACGGTGCCCGTTGGTGGTCAGGGTTACGCGCGAAagc                                          >  1:993270/1‑151 (MQ=255)
                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAAAGCTGCCAGTTGGTGGGCAGGTTTACGCGCAAAcgc                                          >  1:1385964/1‑151 (MQ=255)
                                                                                                             aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGGGGTCAGGGTTACGCGCAAAcgcg                                         >  1:469251/1‑151 (MQ=255)
                                                                                                             aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCATCCATGTCAAGGATCTTGTCGCCCTTTAATTTGAGACGGGGCGGAACGGGGCTTTTTTCCTGTAACACTGCCAGTTTGGGGTCAGGTGTCTGCGCAAAcgcg                                         >  1:900311/1‑151 (MQ=255)
                                                                                                              aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTGCGCCATTTAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGTGGTCAGGTTTTCGCGCAAACGCgg                                        >  2:1180883/1‑151 (MQ=255)
                                                                                                               ttGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCCAGGAACTTTGCGCCATTGAATGTCAGTTGCGGGGTAACGGGGATTTTTTTCTGTTACGCTGCCCGTTGGTGGGCAGGTTTACGCGCAAACGCggg                                       >  2:1388174/1‑150 (MQ=255)
                                                                                                                tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTGCGCCATTGAATGTCAGGTGCGGCGGAACGGCGATTTTTTCCTGTAAAGCCGCCCGTTGGTGGGCAGGGTTTCGCGCAAAAGCGGtt                                      >  2:60976/1‑151 (MQ=255)
                                                                                                                tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCCGTTGCGGCGGAACGGCGATTTTTTTCTGTAACGCCGCCAGTTGGTGGGCAGGGTTTCGCGCAAACGCgggt                                      >  2:190177/1‑149 (MQ=255)
                                                                                                                   gcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGGGGTCAGGGTTACGCGCAAACGCGGGTTGt                                   >  1:395462/1‑151 (MQ=255)
                                                                                                                   gcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGt                                   <  1:1176021/151‑1 (MQ=255)
                                                                                                                     gTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCTGCCCGTTGGGGGGCAGGGTTACGCGCAAACGCGGgttgttt                                 >  1:1178557/1‑150 (MQ=255)
                                                                                                                       aGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtca                                <  2:942966/150‑1 (MQ=255)
                                                                                                                       aGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtca                                <  1:1222377/150‑1 (MQ=255)
                                                                                                                       aGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtca                                <  2:544745/150‑1 (MQ=255)
                                                                                                                        gTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtcatc                              <  1:175968/151‑1 (MQ=255)
                                                                                                                           aCGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAAt                           <  2:1151846/151‑1 (MQ=255)
                                                                                                                            cGCTTGCAGAATTTCATCGGCGGGGTGTTGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATg                          <  2:711254/151‑1 (MQ=255)
                                                                                                                               ttGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTTAAGGATCTTTTCGCCCTTGAATGTCAGTTGGGGCGGAACGGCGATTTTTTTCTGTTACGCTGCCCGTTTGTGGTCCGGTTTTTGCGGAAAAGCGGGTTGTTTCTTAATGccc                       >  2:59541/1‑150 (MQ=255)
                                                                                                                                tGCATAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGACGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGcc                        <  1:973571/149‑1 (MQ=255)
                                                                                                                                tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGcc                        <  1:1001525/149‑1 (MQ=255)
                                                                                                                                tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCa                       <  1:749221/150‑1 (MQ=255)
                                                                                                                                tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCa                       <  1:591976/150‑1 (MQ=255)
                                                                                                                                tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCa                       <  2:790754/150‑1 (MQ=255)
                                                                                                                                tGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCa                       <  1:525381/150‑1 (MQ=255)
                                                                                                                                 gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                     <  1:14209/151‑1 (MQ=255)
                                                                                                                                 gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                     <  1:1448376/151‑1 (MQ=255)
                                                                                                                                 gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAg                     <  1:33748/151‑1 (MQ=255)
                                                                                                                                 gCAGAATTTCATCGGCGGGCTGTCGCCCGTCCAGCTCAAGGATCTTTCCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCCGTTTGTTCATCAATGCCAAg                     <  1:878712/151‑1 (MQ=255)
                                                                                                                                   aGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGGTTGTTTATCAATGCCCAGAc                   >  1:401705/1‑151 (MQ=255)
                                                                                                                                    gAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACg                  <  1:1375156/151‑1 (MQ=255)
                                                                                                                                      aTTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTTTTCATCAATGCCAAGACGa                 >  1:750827/1‑150 (MQ=255)
                                                                                                                                             cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca          <  2:560040/150‑1 (MQ=255)
                                                                                                                                              ggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAAt        >  1:861760/1‑151 (MQ=255)
                                                                                                                                               gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACACTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAAt        >  2:837323/1‑150 (MQ=255)
                                                                                                                                                   ggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg   >  1:435304/1‑151 (MQ=255)
                                                                                                                                                    gtgtCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGc  <  1:469734/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133153‑3133451

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