Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,480 T→C 14.4% K109E (AAA→GAA)  PP_1366 ← transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,560,4800TC14.4% 865.9 / 109.6 329K109E (AAA→GAA) PP_1366transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base T (160/119);  new base C (22/25);  total (183/146)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 346 total aligned reads displayed.

AGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTGGAT  >  NC_002947/1560331‑1560629
                                                                                                                                                      |                                                                                                                                                      
aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                                                                                                                        <  2:4302055/150‑1 (MQ=255)
aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCGTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                                                                                                                        >  2:7097575/1‑150 (MQ=255)
  gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                                                                                                                                      >  2:5006072/1‑150 (MQ=255)
        tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCCACGTGGCTTTCCActct                                                                                                                                                 >  2:6101603/1‑150 (MQ=255)
          ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTc                                                                                                                                             >  2:1269176/1‑150 (MQ=255)
          ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTGGTTGTGCGAAGCAAGAGGGCGTATGCGGAGAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTGTc                                                                                                                                             >  2:3652017/1‑150 (MQ=255)
           tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                                                                                                                            <  1:5006072/150‑1 (MQ=255)
           tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCa                                                                                                                                            <  2:2755104/150‑1 (MQ=255)
              aTCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACGCTTTCAggg                                                                                                                                         >  1:5082731/1‑150 (MQ=255)
                   cTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtg                                                                                                                                    <  2:5823148/150‑1 (MQ=255)
                     gATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCTCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                                                                                                                                  <  2:6403559/150‑1 (MQ=255)
                         ttcgcacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCTCAAGTGGCTTTCCACTCTGTCAGGGTCTTGTGGTTg                                                                                                                              >  2:4082167/8‑150 (MQ=38)
                           aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                                                                                                                            >  1:5520843/1‑150 (MQ=255)
                                aaCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTTCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt                                                                                                                       >  1:4607844/1‑150 (MQ=255)
                                 aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGGGGTTGCCGCCttt                                                                                                                      >  2:6869943/1‑150 (MQ=255)
                                  cAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTCCAGGGTCTTGTGGTTGCCGCCTTTg                                                                                                                     >  2:1781340/1‑150 (MQ=255)
                                  cAAAAACGCCGGCAATGTGCCGGCGATTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGGTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg                                                                                                                     <  1:3155044/150‑1 (MQ=255)
                                   aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg                                                                                                                    >  1:2139522/1‑150 (MQ=255)
                                   aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGAGGCTTTCCACTCTTTCAGGGTCGTGTGGTAGCCGCCTTTgg                                                                                                                    >  2:7127152/1‑150 (MQ=255)
                                       aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCCCTATTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTc                                                                                                                >  2:6311579/1‑150 (MQ=255)
                                       aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc                                                                                                                >  1:1605379/1‑150 (MQ=255)
                                       aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc                                                                                                                >  2:391427/1‑150 (MQ=255)
                                            ggCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGGTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCa                                                                                                           >  2:1950772/1‑150 (MQ=255)
                                               aaTGTGCCGGCGTTTTTGTTTTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCGGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACtt                                                                                                        >  1:1572099/1‑150 (MQ=255)
                                                      cGGCGTTTTTGTTTTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCTAAGAGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTACACCGtt                                                                                                 >  2:4413171/1‑150 (MQ=255)
                                                      cGGCGTTTTTGTTGTGCGAAGCAATAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGtt                                                                                                 >  1:385497/1‑150 (MQ=255)
                                                      cGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGtt                                                                                                 >  2:1035908/1‑150 (MQ=255)
                                                      cGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTCGATCACTTCACCGtt                                                                                                 <  2:2520181/150‑1 (MQ=255)
                                                      cGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTCGATCACTTCACCGtt                                                                                                 <  2:285531/150‑1 (MQ=255)
                                                        gCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCTCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTtgt                                                                                               <  1:1765070/150‑1 (MQ=255)
                                                        gCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTtgt                                                                                               <  1:1316516/150‑1 (MQ=255)
                                                        gCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTACAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCATTCAGGGTCTTGTGGTTGCCGCCGGTGGTTTCGATCACATCACCGTtgt                                                                                               >  2:1672510/1‑150 (MQ=255)
                                                         cGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTtgtg                                                                                              <  2:1633997/150‑1 (MQ=255)
                                                          gTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGc                                                                                             >  1:4237937/1‑150 (MQ=255)
                                                           tttttGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCg                                                                                            >  1:6066006/1‑150 (MQ=255)
                                                             tttGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCGGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCggg                                                                                          <  2:5929843/150‑1 (MQ=255)
                                                                gttgtGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTc                                                                                       >  1:4335097/1‑150 (MQ=255)
                                                                 ttgtGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGTCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCt                                                                                      >  2:5799058/1‑150 (MQ=255)
                                                                  tgtgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTTCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCtt                                                                                     <  2:3606448/150‑1 (MQ=255)
                                                                  tgtgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCtt                                                                                     <  1:3808940/150‑1 (MQ=255)
                                                                  tgtgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTAGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCtt                                                                                     <  2:7322905/150‑1 (MQ=255)
                                                                    tgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCTGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGt                                                                                   >  1:114109/1‑150 (MQ=255)
                                                                    tgCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTTCACTCTTTCAGGGTGTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCAGCGTTGTGCGGGTTCTTGt                                                                                   >  2:5343892/1‑150 (MQ=255)
                                                                      cGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAt                                                                                 <  2:1605379/150‑1 (MQ=255)
                                                                      cGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTAGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAt                                                                                 >  1:6367073/1‑150 (MQ=255)
                                                                       gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTTTAtt                                                                                >  1:5753420/1‑150 (MQ=255)
                                                                       gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                                >  2:3302188/1‑150 (MQ=255)
                                                                       gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                                <  2:4349300/150‑1 (MQ=255)
                                                                       gAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTAtt                                                                                >  1:2324477/1‑150 (MQ=255)
                                                                         aGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCGGCTTACAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGc                                                                              >  1:2731479/1‑150 (MQ=255)
                                                                         aGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGc                                                                              <  2:6606226/150‑1 (MQ=255)
                                                                         aGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTCCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGc                                                                              <  2:776963/150‑1 (MQ=255)
                                                                         aGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGc                                                                              <  1:6763946/150‑1 (MQ=255)
                                                                          gCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCt                                                                             <  2:2634384/150‑1 (MQ=255)
                                                                           cAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTTTATTGCtt                                                                            <  1:4057140/150‑1 (MQ=255)
                                                                            aaGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTg                                                                           >  1:4132501/1‑150 (MQ=255)
                                                                              gagGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGAc                                                                         <  2:3433702/150‑1 (MQ=255)
                                                                                 ggCGTATGCGGATAATCAGTCCAGCAGGGGTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGTCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCtt                                                                      <  2:6520416/150‑1 (MQ=255)
                                                                                   cGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgc                                                                    <  1:5725845/150‑1 (MQ=255)
                                                                                   cGTATGCGGATAATCAGGCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgc                                                                    <  2:2241819/150‑1 (MQ=255)
                                                                                    gTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTAAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCtt                                                                      <  1:1004321/147‑1 (MQ=255)
                                                                                    gTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTAAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCtt                                                                      >  2:1004321/1‑147 (MQ=255)
                                                                                        gCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgcgtgcg                                                               >  1:824191/1‑150 (MQ=255)
                                                                                         cGGATAATCAGTCCAGCAGGGTTGCCGAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTgcgc                                                              >  1:3353774/1‑150 (MQ=255)
                                                                                         cGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTgcgc                                                              <  2:4134142/150‑1 (MQ=255)
                                                                                             tAATCAGTCCAGCAGGGTTGCCCTGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCtt                                                          <  1:186519/150‑1 (MQ=255)
                                                                                              aaTCAGTCCAGCAGGGTTGCCCAGCTTTCACCCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCCGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTg                                                         <  2:462912/150‑1 (MQ=255)
                                                                                               aTCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTAGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTgg                                                        <  1:2412040/150‑1 (MQ=255)
                                                                                               aTCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTgg                                                        >  1:1695204/1‑150 (MQ=255)
                                                                                                  aGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTCGa                                                                                                              <  1:7204661/93‑1 (MQ=255)
                                                                                                   gTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTTCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGAc                                                    >  2:4568229/1‑150 (MQ=255)
                                                                                                            ggTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCtt                                           <  2:1537001/150‑1 (MQ=255)
                                                                                                              ttGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCCGCTTGGTGACGGTAGTCTTGAc                                        <  2:2170781/150‑1 (MQ=255)
                                                                                                                 cccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCTACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCg                                      <  2:5074968/150‑1 (MQ=255)
                                                                                                                 cccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTg                                                         <  2:3809219/131‑1 (MQ=255)
                                                                                                                 cccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCg                                      >  1:3222791/1‑150 (MQ=255)
                                                                                                                 cccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTGGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCg                                      >  1:6961373/1‑150 (MQ=255)
                                                                                                                  ccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTACAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCg                                      >  2:4523266/1‑150 (MQ=255)
                                                                                                                  ccAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGa                                     >  1:4836240/1‑150 (MQ=255)
                                                                                                                   cAGCTTTCAACCATGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTGCGTGCGCGCTTGGTGACGGTAGTCtt                                           >  1:5873440/1‑142 (MQ=255)
                                                                                                                   cAGCTTTCAACCACGTCGCCGCCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGa                                     >  1:6752509/1‑150 (MQ=255)
                                                                                                                   cAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTAGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGaa                                    >  2:3113266/1‑150 (MQ=255)
                                                                                                                    aGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCCTGCGCGCTTGGTGACGGTAGTCTTGACCGAAc                                   <  1:2682160/150‑1 (MQ=255)
                                                                                                                     gCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAAcc                                  <  1:7100217/150‑1 (MQ=255)
                                                                                                                       tttCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTATTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCAc                                >  2:6725108/1‑150 (MQ=255)
                                                                                                                         tCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGa                              <  2:5032860/150‑1 (MQ=255)
                                                                                                                          cAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAg                             <  2:1712405/150‑1 (MQ=255)
                                                                                                                                cGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCt                       >  1:7397528/1‑150 (MQ=255)
                                                                                                                                    gccgccCCAAGTGGCTTTCCACTCTTCCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCt                   >  2:5474149/1‑150 (MQ=255)
                                                                                                                                     ccgccCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTc                  >  2:474686/1‑150 (MQ=255)
                                                                                                                                       gccCCAAGTGGCTTACCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAg                <  2:4786492/150‑1 (MQ=255)
                                                                                                                                         ccccAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTGACCGAACCACGAGGTGCCTTGCTCAGtt              <  1:6583136/150‑1 (MQ=255)
                                                                                                                                             aaGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCg           >  1:2895590/1‑150 (MQ=255)
                                                                                                                                               gTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgc                                                                    >  1:6199014/1‑91 (MQ=255)
                                                                                                                                               gTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTtgc                                                                    <  2:6199014/91‑1 (MQ=255)
                                                                                                                                                 ggCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGActc       <  1:5448024/150‑1 (MQ=255)
                                                                                                                                                   cTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTg     <  2:1131945/150‑1 (MQ=255)
                                                                                                                                                   cTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCCCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTg     >  1:3608575/1‑150 (MQ=255)
                                                                                                                                                    tttCCACTCTTTCAGGGTCTTGTGGTTGACGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTgg    >  1:881944/1‑150 (MQ=255)
                                                                                                                                                    tttCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTgg    >  1:5279800/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCACCGTTGTGCGGGTTCTTGTATTGCTTGACCTTGCGTGCGCGCTTGGTGACGGTAGTCTTGACCGAACCACGAGGTGCCTTGCTCAGTTTCGACTCTGGAT  >  NC_002947/1560331‑1560629

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: