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

AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGC  >  NC_000913/3133153‑3133451
                                                                                                                                                     |                                                                                                                                                         
aaTGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                                                                                                          <  2:301947/151‑1 (MQ=255)
      gTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGcc                                                                                                                                                     >  1:605456/1‑150 (MQ=255)
          ggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                                                                                                <  2:819753/151‑1 (MQ=255)
           gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCa                                                                                                                                               <  1:1092703/151‑1 (MQ=255)
           gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGAGGCGCCCAtac                                                                                                                                                >  1:1038767/1‑148 (MQ=255)
            cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTcc                                                                                                                                                >  2:974982/1‑149 (MQ=255)
              gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGtt                                                                                                                                            <  1:42963/151‑1 (MQ=255)
              gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGt                                                                                                                                             <  2:64520/150‑1 (MQ=255)
              gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGtt                                                                                                                                            >  1:297715/1‑151 (MQ=255)
              gcggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGGCCATCCCGtt                                                                                                                                            >  2:475611/1‑151 (MQ=255)
                  cgcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCCTCCCGt                                                                                                                                             >  1:742605/1‑146 (MQ=255)
                          tACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCttt                                                                                                                                >  2:301731/1‑150 (MQ=255)
                           aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGTTTAAGGAACttt                                                                                                                                >  2:1212817/1‑149 (MQ=255)
                            cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGCCCCTCCAGGTTAAggggact                                                                                                                                  >  1:229318/1‑142 (MQ=255)
                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgcg                                                                                                                             <  1:973891/149‑1 (MQ=255)
                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGcc                                                                                                                          >  1:359428/1‑151 (MQ=255)
                                 aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCAt                                                                                                                        <  1:667066/150‑1 (MQ=255)
                                  aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                      <  1:76324/151‑1 (MQ=255)
                                  aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTg                                                                                                                      <  2:1262047/151‑1 (MQ=255)
                                  aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTTCAAGATCTTTTCGCCCtt                                                                                                                       >  2:391276/1‑150 (MQ=255)
                                   cAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGa                                                                                                                     <  2:1264662/151‑1 (MQ=255)
                                     gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAAt                                                                                                                   <  2:844704/151‑1 (MQ=255)
                                     gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAAt                                                                                                                   <  1:1403101/151‑1 (MQ=255)
                                      tttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAAt                                                                                                                   <  2:1057189/150‑1 (MQ=255)
                                         gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTc                                                                                                                <  1:872133/150‑1 (MQ=255)
                                           gtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTTAAGGATCTTTTCGCCCttttatggca                                                                                                               >  2:167249/1‑141 (MQ=255)
                                           gtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAAGAACTTTTCGCCCTTTAATGTCAg                                                                                                             >  1:188398/1‑150 (MQ=255)
                                              tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAAGGTCTTTTCGCCATTGAATGTccggt                                                                                                            >  2:1347182/1‑145 (MQ=255)
                                              tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCTTCCCGGTCAAGGAGGTTTTGCGCCTTTGATGGCAGtt                                                                                                           >  2:20431/1‑149 (MQ=255)
                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcg                                                                                                        <  2:923102/151‑1 (MQ=255)
                                                 aGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcgg                                                                                                       <  1:947730/150‑1 (MQ=255)
                                                  gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcggcg                                                                                                     <  1:390959/151‑1 (MQ=255)
                                                   tGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTTCGCCATTGGATGTCAGTTgcggcg                                                                                                     >  1:1202979/1‑150 (MQ=255)
                                                   agCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGAGCTTTGCGCCATTGAATGTCAGGTgcggcg                                                                                                     >  2:753631/2‑150 (MQ=255)
                                                    gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTa                                                                                                   <  1:269982/151‑1 (MQ=255)
                                                     cATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                                                  <  2:854119/151‑1 (MQ=255)
                                                      aTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGCTCAAAGATCTTTTGGCCATTTAATGTTAGTTGCGGCGTaa                                                                                                  >  1:696419/1‑150 (MQ=255)
                                                        ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCAGGTCCAAGATCTTTGCGCCCATGGAATTCTGTTTCGGCGTAACgg                                                                                               >  1:1124556/1‑151 (MQ=255)
                                                            tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCCGTTCAAGGGACTTTTCGCCCTTGAATGTCAGTTGCGGGGTGACGGCg                                                                                            >  2:8970/1‑149 (MQ=255)
                                                            tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAAGGTCTTTTTGCCCTTTAATGTAAGGTgcggcggtaaggcggt                                                                                           >  2:805743/1‑141 (MQ=255)
                                                            tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGGTCAAGGATCTTTTCGCCCTTTGATGGCAGGTGCGGCGTAACgg                                                                                               >  2:18982/1‑147 (MQ=255)
                                                             cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAt                                                                                          <  1:496028/150‑1 (MQ=255)
                                                                 aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttttt                                                                                      <  1:969638/150‑1 (MQ=255)
                                                                 aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGTGCCCCACCACTTCAAGGGTTTTTGCGCCATTTAAAGTCCGTTGAGGCGTAAAGGCGAgttttt                                                                                     >  2:1425486/1‑151 (MQ=255)
                                                                 aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGGCCACCCCTATCAACGATCTTTTCGCCATTTAATGTCAGGTGCGGCGTAACGGCGAtttttt                                                                                     >  2:466969/1‑151 (MQ=255)
                                                                  aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTGTGCCCTTGAATGTCCGTTGCGGCGTAAAGGGGGcttttt                                                                                     >  2:83693/1‑150 (MQ=255)
                                                                  aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGTCCCTCCAGGTCAAGGGGCTTTTTGCCATTGAATGTGAAGTGGTGCGTAAAGGGGAgtttttt                                                                                    >  2:461579/1‑150 (MQ=255)
                                                                     catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTg                                                                                 <  2:1218853/151‑1 (MQ=255)
                                                                      atcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGAACTTTGTGCCATTTAATGTCAGTTGTGGCGGAACGGCGGTTTTTTTCTGt                                                                                >  2:749083/1‑151 (MQ=255)
                                                                            tATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTTAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGc                                                                           >  1:1379422/1‑150 (MQ=255)
                                                                            tATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGGCCATCCAGGTCAAGGAGCTTTGCGCCATTGAATGTGAGGTGGGGGGGAAAGGGGAGTTTTTTCTGTaaa                                                                             >  2:937918/1‑147 (MQ=255)
                                                                             aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTg                                                                         <  1:258244/151‑1 (MQ=255)
                                                                              tCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCCCGTCAAAGATCTTTGCGCCATTTAATGTTAGTTTCGGCGGGACGGCGAGTTTTTTCTTTAACGc                                                                           >  2:1308447/1‑148 (MQ=255)
                                                                                  ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGCTCAAAGATCTTTTCGCCCTTTAATGTCAGTTGCGGCGTAACGGCGGATTTTTTCTTTTAAcccgcc                                                                        >  2:1456360/1‑143 (MQ=255)
                                                                                    ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGGTCTTTTCGCCCTTGAATGTCAGTTTCGGCGGAACGGCGGTTTTTTTCTTTTAAGCTGCCCGtt                                                                   >  1:1084411/1‑150 (MQ=255)
                                                                                         gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtggtgg                                                              <  2:799099/150‑1 (MQ=255)
                                                                                         gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAAGATCTTTGCGCCATTTAATGTCAGGTGCGGCGGAACGGCGGTTTTTTTCTGTTAAGCTGGCCGTtgggggg                                                             >  2:1162625/1‑150 (MQ=255)
                                                                                          aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGGTGCGGCGTAAAGGCGGTTTTTTTCTGTAAAGCTGCCAGTtggggg                                                              >  1:75022/1‑149 (MQ=255)
                                                                                           tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGAGCCCAACCAGTTCCAAGGGCTTTGTGCCATTTAATGTCGAGTGGGGGGGGAAGGCGGTTTTTTTCTTTAACGCTGCCCGGtggtgg                                                              >  2:1123580/1‑148 (MQ=255)
                                                                                            cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCa                                                           <  2:64683/150‑1 (MQ=255)
                                                                                            cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAg                                                          <  2:5131/151‑1 (MQ=255)
                                                                                            cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAg                                                          <  1:694536/151‑1 (MQ=255)
                                                                                               gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGAGTTTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGtt                                                       <  1:463758/151‑1 (MQ=255)
                                                                                                  aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAc                                                    <  1:505371/151‑1 (MQ=255)
                                                                                                     gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGGTCAAGGGACTTTTCGCCCTTTGATGTCAGGTGGGGCGTAACGGGGGTTTTTTTCTGTAACGCTGCCCGGTGGTGGTGAGGgttttgcg                                                 >  2:381455/1‑146 (MQ=255)
                                                                                                      cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACgcgc                                                <  1:693051/151‑1 (MQ=255)
                                                                                                       ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGGATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGGGGGCAGGGTTTcgc                                                  >  1:1377288/1‑148 (MQ=255)
                                                                                                        gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaa                                              <  2:1187320/151‑1 (MQ=255)
                                                                                                        gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCaa                                              <  2:1368458/151‑1 (MQ=255)
                                                                                                         cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTTAATGTCAGTTGCGGCGGAACGGCGGTTTTTTTCTGTAAAGCCGCCAGTTGGGGGTCAGGGTTTCGCGCaa                                              >  1:1011675/1‑150 (MQ=255)
                                                                                                           tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGTAAAGGCGGTTTTTTTCTGTAAAGCTGCCAGTTGGGGGGCAGGGTTACGCGCaaa                                             >  1:1051627/1‑149 (MQ=255)
                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAAcgc                                          <  2:359379/151‑1 (MQ=255)
                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGGTCAAGGGACTTTGCGCCATTGAATGTCAGGTGCGGCGGAACGGCGGTTTTTTTCTTTAAAGCTGCCAGGTGGGGGGCAGGGTTAAGCGCAAAcgc                                          >  2:375642/1‑151 (MQ=255)
                                                                                                                tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGt                                       <  1:1440578/150‑1 (MQ=255)
                                                                                                                tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGt                                       <  2:882123/150‑1 (MQ=255)
                                                                                                                 gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCCGTTTCAGGATATTTTCGCCCTTTAATTTCAGATGCGGCGGAAAGGCGATTTTTTTCTGTAACGCTGCCCCGTGGTGGGCAGGGTTTCGCGCAAAAGCgggt                                      >  2:88290/1‑148 (MQ=255)
                                                                                                                 gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTGTAAAGCTGCCAGTTTGGGGGCAGGGTTACGCGCAAACGCGGgtt                                     >  1:526519/1‑151 (MQ=255)
                                                                                                                    cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGtt                                  <  2:297682/151‑1 (MQ=255)
                                                                                                                    cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGtt                                  <  1:475666/151‑1 (MQ=255)
                                                                                                                        gTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGGTCAAGGAACTTTTCGGCAATTAATGTCAGTTGCGGCGGAACGGGGGTTTTTTTCTGTAACACCGCCCGTTGGGGGGCAGGGTTTAGCGCCAAAACGGGTTGTtcat                               >  2:1100805/1‑150 (MQ=255)
                                                                                                                         tGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGAGCTTTGCGCCATTGGATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTGTAACGCTGCCCGTTGGTGGGCAGGGTTACGCGCAAACGCGGTTTGTtcct                               >  2:503746/1‑147 (MQ=255)
                                                                                                                             gCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGc                         <  1:1347289/151‑1 (MQ=255)
                                                                                                                                tGCAGAATTTCATCGGCGGGGGGTACACCCTCCACATCAAGGAAATTTTCGCCGTTGAACGTCTGATTCGGAGGGAAGGGGATTTTTTTCTCTAACGCAGCCAGTTGGTGGTTAAGGTTTCCCGCAAAAGCGGGTTTATCTTCaaatcc                        >  1:648019/1‑145 (MQ=255)
                                                                                                                                  cAGAATTTCATCGGCGGGGTGATGGCCATCCAATTCCAGGAGCCTTTTGCCCTTTAATTTAAGTTGCGGCGGAACAGCGCTTTTTTTCCGTAAAACTGCCAGCTGGTGGGGAGGGTTAAGCGGAAACGCGGGTTGTtgatccaaaa                         >  2:1220756/1‑141 (MQ=255)
                                                                                                                                    gAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTTAAGGGTCTTTTCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTGTAAAGCTGCCCGTTGGGGGTCCGGGTTACGCGCCAACCCGGGTTGTTCATCAATGCCaaaa                    >  1:191930/1‑147 (MQ=255)
                                                                                                                                       tttCATCGGCGGGGTGTTGCCCCTCCCGTTTAAGGAACTTTTCGCCATTGGATGTCAGTTGCGGCGGAACGGCGATTTTTTCCTGTAACGATGCCCGTTGGTGGGTAGGTTTACGCGCAAAAGCGGTTTGTTCCTCAAAGCCCAGACGaa                >  2:345788/1‑150 (MQ=255)
                                                                                                                                          cATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaa            <  2:1038667/151‑1 (MQ=255)
                                                                                                                                          cATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATca             <  2:574880/150‑1 (MQ=255)
                                                                                                                                             cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca          >  1:104216/1‑150 (MQ=255)
                                                                                                                                             cggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaaca          >  2:153366/1‑150 (MQ=255)
                                                                                                                                                  gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACgg    <  2:75014/151‑1 (MQ=255)
                                                                                                                                                    gtgtCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGc  <  1:301766/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 < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: