Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,316,710 G→C 13.9% intergenic (‑223/+55) PP_1144 ← / ← rapA membrane protein/RNA polymerase‑associated protein RapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,316,7100GC13.9% 1334.4 / 66.3 462intergenic (‑223/+55)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base G (160/231);  new base C (59/12);  total (219/243)
Fisher's exact test for biased strand distribution p-value = 2.28e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 472 total aligned reads displayed.

AGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGA  >  NC_002947/1316561‑1316859
                                                                                                                                                      |                                                                                                                                                     
aggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAAGGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCg                                                                                                                                                       <  1:1085206/150‑1 (MQ=255)
  gaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCATACACAGCGGTGGATGCCACCCGc                                                                                                                                                     >  2:2622008/1‑150 (MQ=255)
     tACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGCTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTACTGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCttt                                                                                                                                                  <  2:6095223/150‑1 (MQ=255)
     tACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCttt                                                                                                                                                  <  2:8048330/150‑1 (MQ=255)
       cTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGc                                                                                                                                                <  2:1346099/150‑1 (MQ=255)
        ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGGAGCAAGCCCAGCGGGGGTTGCCCCCggggtttgc                                                                                                                                               >  2:5917829/1‑143 (MQ=255)
        ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGGAGCAAACACAGCGGTGGATGCCACCCGCGTTGcc                                                                                                                                               >  1:3543149/1‑150 (MQ=255)
        ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCCAGCAACCACCGCGGTGGGTGCCCCCGGCGGTGcc                                                                                                                                               >  2:3758613/1‑150 (MQ=255)
        ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCAAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGACGCCACCGGCTTTGcc                                                                                                                                               >  2:3517481/1‑150 (MQ=255)
         tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                              <  1:7722890/150‑1 (MQ=255)
         tAAGCCTTCCTGACTAATGGCACAATGAGTTCTACGATTGGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                              <  2:7411295/150‑1 (MQ=255)
            gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtg                                                                                                                                           <  1:4024563/150‑1 (MQ=255)
            gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACACCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtg                                                                                                                                           <  1:890095/150‑1 (MQ=255)
             ccttccttACTAATGGCACAATGATTTCTACGATTTGGTTATCACTCCATTCGCTTTGAGTAATAATTTTGGGTCCCTAATGCCTTCCTGTGGTAGCGGGCGAGCCCGCGAAGCATCCACTGCGTTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                                          <  1:1007003/150‑1 (MQ=38)
             ccTTCCTGACTAATGGCACAATGATTTCTCCGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                                          <  1:8601889/150‑1 (MQ=255)
               ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATTACTCCCATCGCAGTTAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCTCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTc                                                                                                                                        <  2:8744893/150‑1 (MQ=38)
               ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAATTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtggca                                                                                                                                        <  2:4838462/150‑4 (MQ=255)
                 ccTGAATAATGGCACAATTATTTCTACTTTTTTTTTATAACTCCAATCGCAGTGAGTAATAATTTGGGGGCCATAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTcgc                                                                                                                                      <  1:4544145/150‑1 (MQ=37)
                   tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCTGGGAGGCATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGACGCCGCCCATGCGGTGGAGGCCCCCGGCTGTGGCGGGGGTGGCgg                                                                                                                                    >  2:5581238/1‑150 (MQ=255)
                   tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                    <  2:6393402/150‑1 (MQ=255)
                   tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCCACCAACCACAGCCGGGGATGCCACCCGCGTTGCCGGGGTTCGCgg                                                                                                                                    >  2:3205675/1‑150 (MQ=255)
                     aCTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGc                                                                                                                                  <  2:2721992/150‑1 (MQ=255)
                         aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgc                                                                                                                              >  2:4287918‑M2/1‑148 (MQ=255)
                         aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgc                                                                                                                              >  2:7282154‑M2/1‑148 (MQ=255)
                           ggCACAATGATTTCTACGATTGGGTTATAACTCCACGCGCAGTGCGGAATAAATTTGGGGCCCTACGGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgccc                                                                                                                            >  1:1736134‑M2/1‑146 (MQ=255)
                              acaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAAACAAAGCGGTGGATTCCACCCGCATTGCCCGTGTTCGCCGGCCGGccccca                                                                                                                         >  1:8815625/1‑147 (MQ=255)
                                  tGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccca                                                                                                                     >  1:1319074‑M2/1‑139 (MQ=255)
                                  tGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCCCGCATCAACCACAGCGGTGGTTGGCACCCGCGTTTGCGGGGTTTGCGGGGGGGGCCGCGCCca                                                                                                                     >  2:2264931/1‑150 (MQ=255)
                                  tGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCCAAGCAAACACAGCCGTGGATGCCACCCGCGTTTCCGGTGTTGGCGGGCGTgccccctcccc                                                                                                                     <  2:5478845‑M2/150‑12 (MQ=255)
                                         tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGGGGAGCATCCACTGCGGGGGGTGCCCCCCGCGTTTCCGGTGGTCGCGGGGGTGGCCGCACCCACAAaggg                                                                                                              >  2:1925573/1‑150 (MQ=255)
                                         tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCCTCCACTGCGGGGGGAGCCACCCGCATTTCCGGGGTTGGCGGGGGGGGCCCCACCCACAAgggg                                                                                                              >  1:2136236/1‑150 (MQ=255)
                                         tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACAGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgggg                                                                                                              >  2:7943119/1‑150 (MQ=255)
                                         tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAGCAACCACAGCGGTGGGTTCCACCCGCGTTGGCCGTGTTGGCGGGCGGGCCCCCACCCACaaacgg                                                                                                              >  1:6412587/1‑146 (MQ=255)
                                            gATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGgcgc                                                                                                           <  2:1566927/150‑1 (MQ=255)
                                              ttgggTTATAACGCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                                         <  2:8662755/147‑1 (MQ=255)
                                               tgggTTAAAACTCCAATCGCAGTGAGTAATAATTTTGGGGCCCTAATGCCTTCCTGGGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCggg                                                                                                        <  1:6069889/148‑1 (MQ=255)
                                                 ggTTAGACCGCCAATCGCAGGGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGGGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGc                                                                                                      <  1:6974024/150‑1 (MQ=37)
                                                  gTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCCCCCACCAAACACAGCGGTGGGTGCCCACCCCCTTTACCGTGTTGGCGGGCGGGCCCCCGACCCCAAAGGGCGCCGGGcc                                                                                                     >  1:526832/1‑150 (MQ=255)
                                                  gTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGAGGGAGCGGGCGAGCCCCCCAACCAAACACCGCGGTGGTTGCCCCCGGCGTTGGCGGTGTTCGCGGGCGTGCCCGCACCCACAAGGGGCGCgggggc                                                                                                     >  2:7038573/1‑148 (MQ=255)
                                                     ataACTCCAAGCGCAGTGCGGAATAACTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTgc                                                                                                  <  2:7715098/150‑1 (MQ=255)
                                                          tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCATCCACAGCGGGGGATTGCACCCGCGGTTCCGGTGTTCGCGGGCGTGCCCGCGCCCCCCAGGGGCGCGGGGCCAGCGCTag                                                                                             >  2:6658031/1‑150 (MQ=255)
                                                           ccACTCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCGTTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagc                                                                                            <  1:2027307/150‑1 (MQ=255)
                                                            cAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCGGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGAGCGCAGGCCCCgcgcccctt                                                                                           <  2:8313335/150‑6 (MQ=37)
                                                             cctCGCTGTGCGGAATAATTTGGGGGCCCTATTGCCTTCGTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                                          <  1:1785343/148‑1 (MQ=25)
                                                                 gCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAACAAACACAGCCGTGGATGCCACCCGCTTTTCCCGTGTTCGCGGGCGTGCCCCCAACCACAAAGGGCGCGGGGCCAGCGCAAgccaccca                                                                                      >  1:6587587/1‑150 (MQ=255)
                                                                  ccgTGAGTAATAAATTTGGGGCCCTAAGGCCTGCCTTTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccac                                                                                     <  1:577472/148‑1 (MQ=255)
                                                                  ccgTGAGGAATAAATTTGGGGCCCCAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccac                                                                                     <  1:4220798/148‑1 (MQ=255)
                                                                   aGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccacc                                                                                    <  1:686459/150‑1 (MQ=255)
                                                                          aaTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCAcgc                                                                             <  2:7836365/150‑1 (MQ=255)
                                                                          aaTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAACCAAACACTGCGGTGGATTCCACCCGCATTTCCCGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCAGCGCAAGCCAGCCACCACCAcgc                                                                             >  2:5705854/1‑150 (MQ=255)
                                                                              attttGGGGCCCTATGGCGTTCGTGGGGGCGCGGGCGACCCCGCGAAGCATCCACGGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                         <  1:5914358/148‑1 (MQ=37)
                                                                              aaTTGGGGGCCCTAATGCCTTCCTGGGGGCGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                         <  1:3751506/150‑1 (MQ=255)
                                                                               aTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATg                                                                        <  2:6511214/150‑1 (MQ=255)
                                                                                         cTACTGCCTTCCTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGgc                                                              <  1:2518548/150‑1 (MQ=255)
                                                                                          tAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACAGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCCCCCGCACGCCAATGGCTTCCAGgcg                                                             >  2:7173529/1‑150 (MQ=255)
                                                                                              gCCTTCCTGTGGGAGCGGGCGAGGCCGGGAAGCATACAAAGCGGTGGATTGCACCCGCGTTTGCGGTGTTGGCGGGCGTGGCCGCGCCCCCCACGGGCGCGGGGGCGGCGCGAGCCCGCCCCCACCACCCCCATGGCGTCCCGGCGCacg                                                         >  2:8207295/1‑148 (MQ=255)
                                                                                              gCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAgg                                                         <  2:3646113/150‑1 (MQ=255)
                                                                                              gCCTTCCTGGGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAgg                                                         <  1:2697592/150‑1 (MQ=255)
                                                                                                cTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggc                                                       <  2:4416531/150‑1 (MQ=255)
                                                                                                   ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                                    <  1:570160/150‑1 (MQ=255)
                                                                                                    cTGTGGGCGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGcc                                                   <  2:2759941/150‑1 (MQ=255)
                                                                                                        gggAGCGGGCGAGCCCCCCAAGCAACCACTGCGGTGGATGCCACCCGCGTTGCCGGTGTTCGCCGGCGTgcccgctcccacaggaaggcaatagggccccaaatgtattcctcacggcgattggagttataaccaaatcgtagaaatcat                                               <  2:384789‑M2/150‑82 (MQ=255)
                                                                                                          gAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTcc                                             <  2:1684257/150‑1 (MQ=255)
                                                                                                             cGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTACAGGCGCAGGGCGGCCTTCTCCAGc                                          <  2:4166340/150‑1 (MQ=255)
                                                                                                              gggCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCa                                         <  2:7795051/150‑1 (MQ=255)
                                                                                                                gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgcgcccacaggaaggcattagggccccaaaattattactaacagcgattggggttataaccaaaaaggagaaataattgttccat                                       <  1:3219919‑M2/150‑90 (MQ=255)
                                                                                                                gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccat                                       <  1:3619589‑M2/150‑90 (MQ=255)
                                                                                                                  gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccatta                                     <  1:4932798‑M2/150‑92 (MQ=255)
                                                                                                                  gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc                                     >  1:2599290/1‑150 (MQ=255)
                                                                                                                  gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc                                     <  2:1817656/150‑1 (MQ=255)
                                                                                                                  gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGGGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCCGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAGTGc                                     >  2:1084828/1‑150 (MQ=255)
                                                                                                                     cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                  <  1:8024208/150‑1 (MQ=255)
                                                                                                                     cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                  <  2:2540041/150‑1 (MQ=255)
                                                                                                                     cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                  <  2:2881610/150‑1 (MQ=255)
                                                                                                                      ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAgg                                 >  2:2451273/1‑150 (MQ=255)
                                                                                                                      ccGCGAAACATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgctcccacaggaaggcagtagggccccaaatttattcctcacggcgattggagttataaacaaatcgtagaaatcattgtgccattagtca                                 <  1:2325135‑M2/150‑96 (MQ=255)
                                                                                                                       cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCCGCCCCCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc                                >  1:4852694/1‑150 (MQ=255)
                                                                                                                       ccggtaGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc                                >  2:5222169/6‑150 (MQ=255)
                                                                                                                        gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgctcccaccggaaggcagtagggccccaaatttattcctcactgcgattggaggtataaccaactcgtagaaatccttgtgccagtagtcagg                               <  2:499591‑M2/150‑98 (MQ=255)
                                                                                                                        gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc                               >  1:1500721/1‑150 (MQ=255)
                                                                                                                        gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc                               >  2:6188230/1‑150 (MQ=255)
                                                                                                                        gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc                               >  1:7372909/1‑150 (MQ=255)
                                                                                                                         cgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaag                                                                                                              >  2:3564052‑M2/1‑52 (MQ=255)
                                                                                                                           tggCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATc                            >  2:5044772/3‑150 (MQ=255)
                                                                                                                               aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg                        >  2:2075685/1‑150 (MQ=255)
                                                                                                                               aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg                        >  1:5037692/1‑150 (MQ=255)
                                                                                                                                  cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCgcgc                     >  1:7534358/1‑150 (MQ=255)
                                                                                                                                  cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCgcgc                     >  2:3399681/1‑150 (MQ=255)
                                                                                                                                  cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCgcgc                     >  2:1452294/1‑150 (MQ=255)
                                                                                                                                  cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCgcgc                     >  2:3024830/1‑150 (MQ=255)
                                                                                                                                     tGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGcgccgc                  >  1:6652987/1‑150 (MQ=255)
                                                                                                                                        ggTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTg               <  1:8719520/150‑1 (MQ=255)
                                                                                                                                           ggATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTgcgc            <  2:4559096/150‑1 (MQ=255)
                                                                                                                                            gATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCa           <  1:3024830/150‑1 (MQ=255)
                                                                                                                                             aTGCCACCGGCTGTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAg          <  2:6888900/150‑1 (MQ=255)
                                                                                                                                              tGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAgg         <  1:8150766/150‑1 (MQ=255)
                                                                                                                                              tGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAgg         <  2:3891110/150‑1 (MQ=255)
                                                                                                                                                  cACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCg      >  1:7150021/1‑150 (MQ=255)
                                                                                                                                                   aCCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGt     >  2:1181739/1‑150 (MQ=255)
                                                                                                                                                      ggCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGa  <  1:4300127/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGA  >  NC_002947/1316561‑1316859

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