Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,843,613 T→C N408S (AAC→AGC)  gacS ← sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,6130TC98.3% 1447.1 / 4.9 415N408S (AAC→AGC) gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (184/224);  total (188/227)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 416 total aligned reads displayed.

AGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAAC  >  NC_002947/1843466‑1843762
                                                                                                                                                   |                                                                                                                                                     
aGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTg                                                                                                                                                     >  2:1050125/1‑150 (MQ=255)
aGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTg                                                                                                                                                     <  2:1807268/150‑1 (MQ=255)
  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGt                                                                                                                                                   <  2:1660374/150‑1 (MQ=255)
  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGt                                                                                                                                                   >  2:3157343/1‑150 (MQ=255)
  tgtgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGt                                                                                                                                                   >  2:2510810/1‑150 (MQ=255)
    tgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCTGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGa                                                                                                                                                 <  1:5277815/150‑1 (MQ=255)
    tgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGa                                                                                                                                                 <  1:1015227/150‑1 (MQ=255)
    tgCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGa                                                                                                                                                 <  1:4257909/150‑1 (MQ=255)
       gcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGa                                                                                                                                              >  1:232944/1‑150 (MQ=255)
            tCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGc                                                                                                                                         <  1:6307482/150‑1 (MQ=255)
            tCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGc                                                                                                                                         >  2:2691174/1‑150 (MQ=255)
            tCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGc                                                                                                                                         <  1:3255649/150‑1 (MQ=255)
            tCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGc                                                                                                                                         >  1:1610435/1‑150 (MQ=255)
                gAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGa                                                                                                                                     <  1:4379375/150‑1 (MQ=255)
                     gACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGTGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGc                                                                                                                                >  2:7070907/1‑150 (MQ=255)
                      aCAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCa                                                                                                                               >  1:5876618/1‑150 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                             >  2:1927162/1‑150 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                             >  2:480874/1‑150 (MQ=255)
                        aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAAc                                                                                                                             >  1:4957003/1‑150 (MQ=255)
                          gCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACgg                                                                                                                           >  1:7396107/1‑150 (MQ=255)
                            cGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGt                                                                                                                         >  2:5433329/1‑150 (MQ=255)
                               tGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGc                                                                                                                      <  2:2549452/150‑1 (MQ=255)
                                gCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGcc                                                                                                                     <  2:2844054/150‑1 (MQ=255)
                                    gtgtCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTAcc                                                                                                                 <  2:4479007/150‑1 (MQ=255)
                                         cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGa                                                                                                            >  1:458603/1‑150 (MQ=255)
                                         cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGa                                                                                                            <  2:3824377/150‑1 (MQ=255)
                                         cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGa                                                                                                            >  1:6930152/1‑150 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                           >  1:77828/1‑150 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                           <  1:7905904/150‑1 (MQ=255)
                                          tGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                           >  1:2395262/1‑150 (MQ=255)
                                             aCGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCa                                                                                                        <  2:371915/150‑1 (MQ=255)
                                             aCGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCa                                                                                                        <  2:393798/150‑1 (MQ=255)
                                                cTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                                                     <  2:2295342/150‑1 (MQ=255)
                                                cTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCACACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                                                     <  1:1923835/150‑1 (MQ=255)
                                                 tGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtg                                                                                                    >  2:1331352/1‑150 (MQ=255)
                                                 tGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtg                                                                                                    >  2:7155283/1‑150 (MQ=255)
                                                  gATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                   <  2:6164037/150‑1 (MQ=255)
                                                  gATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                   >  2:2285609/1‑150 (MQ=255)
                                                  gATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                   <  2:3497243/150‑1 (MQ=255)
                                                   aTGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgtg                                                                                                  >  2:6602388/1‑150 (MQ=255)
                                                       gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTcgc                                                                                              >  2:6342021/1‑150 (MQ=255)
                                                       gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTcgc                                                                                              >  1:5396120/1‑150 (MQ=255)
                                                         aaCTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                                                            >  2:6037040/1‑150 (MQ=255)
                                                          aCTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGt                                                                                           <  1:7486887/150‑1 (MQ=255)
                                                               gcgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATc                                                                                      <  2:7447599/150‑1 (MQ=255)
                                                               gcgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATc                                                                                      >  1:3907411/1‑150 (MQ=255)
                                                                cgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                     <  2:1463003/150‑1 (MQ=255)
                                                                  cAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAgg                                                                                   <  1:6578643/150‑1 (MQ=255)
                                                                  cAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAgg                                                                                   <  2:5692650/150‑1 (MQ=255)
                                                                      cTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCa                                                                               <  1:3702464/150‑1 (MQ=255)
                                                                            gCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAAtt                                                                         >  1:5053446/1‑150 (MQ=255)
                                                                              tCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCa                                                                               <  2:3384785/142‑1 (MQ=255)
                                                                              tCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCa                                                                               >  1:3384785/1‑142 (MQ=255)
                                                                              tCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTcc                                                                       >  1:4101783/1‑150 (MQ=255)
                                                                              tCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTcc                                                                       <  2:2569933/150‑1 (MQ=255)
                                                                               cGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCa                                                                      <  2:1570352/150‑1 (MQ=255)
                                                                                gTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAg                                                                     <  2:1999083/150‑1 (MQ=255)
                                                                                 tCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgc                                                                    >  2:6839689/1‑150 (MQ=255)
                                                                                   ttCCATCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctg                                                                  >  2:4302248/1‑150 (MQ=255)
                                                                                   ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctg                                                                  >  2:2041698/1‑150 (MQ=255)
                                                                                   ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctg                                                                  >  1:7447727/1‑150 (MQ=255)
                                                                                   ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctg                                                                  >  2:408926/1‑150 (MQ=255)
                                                                                   ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctg                                                                  <  2:3006790/150‑1 (MQ=255)
                                                                                   ttCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgc                                                                 >  1:2968907/1‑150 (MQ=255)
                                                                                       aGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTg                                                              <  1:3058556/150‑1 (MQ=255)
                                                                                       aGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTg                                                              >  1:7237410/1‑150 (MQ=255)
                                                                                                  cACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggcggc                                                   <  2:7316394/150‑1 (MQ=255)
                                                                                                   acgacgATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGcc                                                  <  2:1178409/150‑1 (MQ=255)
                                                                                                      acgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGa                                               <  2:3699134/150‑1 (MQ=255)
                                                                                                      acgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGa                                               <  1:7072543/150‑1 (MQ=255)
                                                                                                        gaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCATGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGc                                             <  2:641486/150‑1 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                           <  1:1741243/150‑1 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                           <  1:3794268/150‑1 (MQ=255)
                                                                                                          tGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGa                                           <  2:2183782/150‑1 (MQ=255)
                                                                                                             tGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGa                                        <  1:4770512/150‑1 (MQ=255)
                                                                                                                cTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATgg                                     <  1:2706611/150‑1 (MQ=255)
                                                                                                                cTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATgg                                     >  2:358181/1‑150 (MQ=255)
                                                                                                                 ttCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGt                                    >  2:7827339/1‑150 (MQ=255)
                                                                                                                 ttCGCGGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGt                                    >  1:6713811/1‑150 (MQ=255)
                                                                                                                   cgcgGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCa                                  <  2:7626381/150‑1 (MQ=255)
                                                                                                                   cgcgGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCa                                  <  1:5475628/150‑1 (MQ=255)
                                                                                                                     cgGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAgg                                <  1:278969/150‑1 (MQ=255)
                                                                                                                     cgGGTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAgg                                <  1:5345399/150‑1 (MQ=255)
                                                                                                                      gggTGAACTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAggg                               <  1:1214747/150‑1 (MQ=255)
                                                                                                                            aCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTgg                                                        <  1:6078782/119‑1 (MQ=255)
                                                                                                                            aCTTGATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCt                         <  1:5631192/150‑1 (MQ=255)
                                                                                                                                gATGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGAt                     <  1:4492363/150‑1 (MQ=255)
                                                                                                                                 aTGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATc                    <  2:6737149/150‑1 (MQ=255)
                                                                                                                                  tGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCa                   <  1:2051083/150‑1 (MQ=255)
                                                                                                                                  tGGCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGATGCCGGAGCGAGGATGGTCAGGGTGTCCTGTATCa                   <  2:2383583/150‑1 (MQ=255)
                                                                                                                                   ggCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCTGGAGCGAGGATGGTCAGGGTGTCCTGGATCAg                  <  1:148044/150‑1 (MQ=255)
                                                                                                                                   ggCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAg                  <  1:6448090/150‑1 (MQ=255)
                                                                                                                                   ggCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAg                  >  1:7515595/1‑150 (MQ=255)
                                                                                                                                   ggCGTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGGGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAg                  <  2:4395043/150‑1 (MQ=255)
                                                                                                                                      gTTGCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTc               <  1:2214094/150‑1 (MQ=255)
                                                                                                                                         gCTGACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCg            <  2:5062447/150‑1 (MQ=255)
                                                                                                                                            gACCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAg         <  1:6152428/150‑1 (MQ=255)
                                                                                                                                             aCCAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAgg        >  1:5946007/1‑150 (MQ=255)
                                                                                                                                               cAGGCTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGTTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGtt      >  2:5184020/1‑150 (MQ=255)
                                                                                                                                                   ttGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAAc  <  1:2413326/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAAC  >  NC_002947/1843466‑1843762

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