Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,316,690 C→A 7.0% intergenic (‑203/+75) 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,6900CA 7.0% 1338.3 / 14.7 453intergenic (‑203/+75)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base C (145/256);  new base A (46/3);  total (194/259)
Fisher's exact test for biased strand distribution p-value = 1.12e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 459 total aligned reads displayed.

TGTCCAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGC  >  NC_002947/1316541‑1316838
                                                                                                                                                     |                                                                                                                                                     
tGTCCAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATc                                                                                                                                                       <  1:167538/150‑1 (MQ=255)
tGTCCAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATc                                                                                                                                                       <  2:8500945/150‑1 (MQ=255)
  tCCAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCa                                                                                                                                                     <  2:8337270/150‑1 (MQ=255)
  tCCAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCa                                                                                                                                                     >  2:6274690/1‑150 (MQ=255)
   ccAGAGGCATTTTAGTAAGGAGTTCTTAAGCCTTCCTGACTATTGTCACATTGATTTCTACGTTTTGTTTATAACTCCAATCGCAGTGAGGAATAATTTTGGGGCCCTAAGGCCTTCCTGTGGGTGCGGGCGAGCCCGCGAAGCATCCAc                                                                                                                                                    <  1:6256598/150‑1 (MQ=255)
    cAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACt                                                                                                                                                   <  1:6051987/150‑1 (MQ=255)
     agagGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTg                                                                                                                                                  <  2:2017797/150‑1 (MQ=255)
        ggCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTATTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAATTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCgg                                                                                                                                               <  1:3181516/150‑1 (MQ=255)
         gCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTCCGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGt                                                                                                                                              <  1:1753164/150‑1 (MQ=255)
          cAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTg                                                                                                                                             <  1:3349143/150‑1 (MQ=255)
            gTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGa                                                                                                                                           <  2:2635828/150‑1 (MQ=255)
            gTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGa                                                                                                                                           <  2:2713213/150‑1 (MQ=255)
            gTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGa                                                                                                                                           <  2:7356720/150‑1 (MQ=255)
               taGGAAGGAGTACTTAAGCCTTCCTGACTAATGGAACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATGc                                                                                                                                        <  1:284905/150‑1 (MQ=255)
               taGGAAGGAGTACTTAAGCCATCCTGACTAATGGCACAATGAGTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGCAATAAATGTGGGGCCCTCATGCCTACCTGTGGGAGCGGGCGAGCCCCCGCAGCACCCCCGGCGGTggcggc                                                                                                                                        >  2:5218312/1‑146 (MQ=255)
                 ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTTATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATGctt                                                                                                                                      <  2:7005659/150‑3 (MQ=255)
                 ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                      >  2:5984094/1‑150 (MQ=255)
                  gaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGCAGCAACCACAGCGGGGGATGCCAc                                                                                                                                     >  1:7239244/1‑150 (MQ=255)
                  gaaggaGTACTTAAGCCTTCCTAACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGCAGCATCCACGGCGGGGGGTGCCAc                                                                                                                                     >  2:3870083/1‑150 (MQ=255)
                            ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGGAGCCTCCCCGGCGGGGGATGCCACCGGCGTTGcc                                                                                                                          >  2:7871121/1‑150 (MQ=255)
                            ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCCAGCAACCACCGCGGTGGGTGCCCCCGGCGGTGcc                                                                                                                          >  2:3758613/1‑150 (MQ=255)
                            ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCAAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGACGCCACCGGCTTTGcc                                                                                                                          >  2:3517481/1‑150 (MQ=255)
                            ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTGGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGcc                                                                                                                          <  1:2674592/150‑1 (MQ=255)
                             tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                         <  2:8804717/150‑1 (MQ=255)
                             tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                         <  2:6111207/150‑1 (MQ=255)
                             tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                         <  2:3412195/150‑1 (MQ=255)
                             tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                         <  1:8052387/150‑1 (MQ=255)
                              aaGCCTTCCTGACTAATTGAAAAATGATTTCTACTATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCgg                                                                                                                        <  1:4126456/150‑1 (MQ=255)
                              aaGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGAAGCATGCACTTCGGTGGTTTGCACCGGCTTTGGCgg                                                                                                                        >  2:6733404/1‑150 (MQ=255)
                               aGCCTTCCTGACTAATGGCACAATGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                       <  1:6818957/150‑1 (MQ=255)
                                gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCATCCACGGCGGTGGCTGCCACCGGCGTTGCCgggg                                                                                                                      >  2:5325995/1‑148 (MQ=255)
                                 ccttccttACTAATGGCACAATGATTTCTACGATTTGGTTATCACTCCATTCGCTTTGAGTAATAATTTTGGGTCCCTAATGCCTTCCTGTGGTAGCGGGCGAGCCCGCGAAGCATCCACTGCGTTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                     <  1:1007003/150‑1 (MQ=38)
                                 ccTTCCTGACTAATGGCACAATGATTTCTCCGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                     <  1:8601889/150‑1 (MQ=255)
                                 ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGCGGCCTCCCCGGCGGGGGATGGCCCCCGCGTGGCCGgtgt                                                                                                                     >  2:8257014/1‑150 (MQ=255)
                                 ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                     <  2:1113794/150‑1 (MQ=255)
                                   ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTc                                                                                                                   <  2:7230750/150‑1 (MQ=255)
                                      cTGACTAATGGCACAATTTTTTCTACTATTTGGTTATAACTCCAATCGCAGTGAGGAATAATTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTcgcg                                                                                                                <  1:6629128/150‑1 (MQ=255)
                                      cTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAAACACAGCGGGGGATGCCACCGGCGTTGCCGGTGTTcgcg                                                                                                                >  2:8150766/1‑150 (MQ=255)
                                      cTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAACATACACAGAGGTGGAATCCAACCCCGTTGCCCGTGTTagcg                                                                                                                >  2:189067/1‑150 (MQ=255)
                                       tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGGAGCAACCACAGGCGGGGGTGGCACCCGCGTTGGCGGTGTTCGCgg                                                                                                               >  1:5703780/1‑150 (MQ=255)
                                       tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCCACCAACCACAGCCGGGGATGCCACCCGCGTTGCCGGGGTTCGCgg                                                                                                               >  2:3205675/1‑150 (MQ=255)
                                         aatAATGGCAATTTGATTTCTAAGATTTGGTTATAACTCCACTCGCTTTGAGTAATAAATTTGGGGCCCTATTGCCTTCCTTTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGc                                                                                                             <  1:2191354/148‑1 (MQ=25)
                                           tAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACAGCGGGGGGTGCCACCGGCGTTGCCGGTGTTAGCGGGCGt                                                                                                           >  1:6561972/1‑150 (MQ=255)
                                             aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCATTGCCGGTGTTAGCGGGCGTgc                                                                                                         <  2:4300127‑M2/150‑3 (MQ=255)
                                             aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAACCCCAGCGGGGGGAGCCACCCGCGGTTGCCGGGTTGGCGGGGGTgc                                                                                                         <  2:453262‑M2/150‑3 (MQ=255)
                                                  acaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAAACAAAGCGGTGGATTCCACCCGCATTGCCCGTGTTCGCCGGCCGGccccca                                                                                                    >  1:8815625/1‑147 (MQ=255)
                                                   caATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAAAAACAGCGGTGGATTCCAACCGCGTTTGCGGTGTTAGCGGGCGTgcccccac                                                                                                   <  2:1782402‑M2/150‑9 (MQ=255)
                                                    aaTGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCATCCACTGCGGTGGCTGCCACCGGCGTTGCCGGGGTTCGCGGGCGGgcccgcgcc                                                                                                  <  2:5640963‑M2/150‑10 (MQ=255)
                                                     aTGATTTCTACGATTTGTTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccc                                                                                                 >  1:3545798‑M2/1‑140 (MQ=255)
                                                      tGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGGGGATGCCACCGGCGTTGCCGGTGTTCGCGGGCGTgcccgctcccc                                                                                                <  2:1387435‑M2/150‑12 (MQ=255)
                                                        aTTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAACCACAGCGGGGGGTGCCCCCCGCGTTTCCGGGGTTAGCGGGCGGgcccccacccaca                                                                                              <  1:1447521‑M2/150‑14 (MQ=255)
                                                         tttCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACaa                                                                                             <  2:7660231/150‑1 (MQ=255)
                                                             tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCATACACGGCGGGGGGTGCCACCCGCGTTGCCCGGGGTTGCGGGCGTGGCCGCACCCCCCAgggg                                                                                         >  1:2584223/1‑150 (MQ=255)
                                                             tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgggg                                                                                         >  2:7966598/1‑150 (MQ=255)
                                                                gATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGgcgc                                                                                      <  2:1566927/150‑1 (MQ=255)
                                                                gATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGgcgc                                                                                      <  1:4987184/150‑1 (MQ=255)
                                                                  tttGGTTATACCTCCAATCGCAGTGTTGAATAATTTTGGGGGCCTACTGCGTTGGTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                    <  1:4803488/150‑1 (MQ=25)
                                                                  tttGGTCATAACTCCACTCGCCGTGAGGAATAAATTTGGGGCCCTCCGGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                    <  1:1916524/150‑1 (MQ=37)
                                                                  ttgggTTATAACGCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                    <  2:8662755/147‑1 (MQ=255)
                                                                     ggTTATACCTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGc                                                                                 <  2:7759930/150‑1 (MQ=255)
                                                                     ggTTAGACCGCCAATCGCAGGGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGGGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGc                                                                                 <  1:6974024/150‑1 (MQ=37)
                                                                         atcaccccAATCGCAGTGAGGAATAAATTGGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTgc                                                                             <  1:4241073/144‑1 (MQ=255)
                                                                         ataACTCCAATCGCCGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGGGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTgc                                                                             <  1:8211331/150‑1 (MQ=255)
                                                                            aCTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCt                                                                          <  2:232875/150‑1 (MQ=255)
                                                                             cTCCACTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTa                                                                         <  2:2529795/150‑1 (MQ=255)
                                                                              tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGAAGCATACACTGCGGTGGATGGCGCCGGGTTTTGGGGTGTTGGCGGGGGTGGCCGCTCCCCAAAAGGGGGCGGGGGCAGCGCTag                                                                        >  2:7086689/1‑150 (MQ=255)
                                                                              tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCCCCCAACAAACAAAGCGGTGGTTGGCCCCCGCCGTTCCGGGGTTGGGGGGCGGGGCCCCTCCCACCAAGGGGGCGGGGGCGGCGCGag                                                                        >  2:1001572/1‑150 (MQ=255)
                                                                                cAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCCCGGAGCCAAAACTGCGGTGGGTGCCACCCGCGTTTGCGGTGTTGGCGGGCGGGCCCCCGCCCACCACGGGGGCGGGGGCGGCGCTagcc                                                                      >  2:3019132/1‑150 (MQ=255)
                                                                                  aTCGCTTTGTGGAATAATTGGGGGGCCCTAAGGCCTGCGGGTGGGCGGGGGCGCCCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagccag                                                                    <  1:4452290/150‑1 (MQ=25)
                                                                                  aTCGCTGTTAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGAGCGCAGGCCCCgcgccccttgt                                                                    <  1:2494341/150‑8 (MQ=25)
                                                                                  aTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGAGCGCAGGCCCCgcgccccttgt                                                                    <  2:2912865/150‑8 (MQ=25)
                                                                                      cAGTGAGGAATAAATTTGGGGCCCTACTGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccac                                                                <  1:577769/150‑1 (MQ=255)
                                                                                                  aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGGGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                    <  2:5796879/150‑1 (MQ=255)
                                                                                                     ttGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGc                                                 <  2:6582133/150‑1 (MQ=255)
                                                                                                       ggggCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCATCCCCAGCGGTGGCTGCCACCCGCTTTGCCGGGGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCGAGCCAGCCACCAGCCCCCGCAAGGCtt                                               >  2:7813858/1‑150 (MQ=37)
                                                                                                                     ttCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                 <  1:3854479/150‑1 (MQ=255)
                                                                                                                        cTGTGGGCGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGcc                              <  2:2759941/150‑1 (MQ=255)
                                                                                                                        cTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGcc                              <  1:371413/150‑1 (MQ=255)
                                                                                                                             ggAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTtctc                         <  2:1212705/150‑1 (MQ=255)
                                                                                                                              gAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTcc                        <  2:1684257/150‑1 (MQ=255)
                                                                                                                                gCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACCAGGGGCGCGGGGCCTGCGCGAGCCAGCCCCCAGCACCCGACTGGCTTCCGGGCGCAGGGCGGCCTTCTcccc                      >  2:643998/1‑148 (MQ=37)
                                                                                                                                   ggCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccccaggaaggccttagggcccccaaggtattcctcacggcgattggagtgataacccaatcggagaaatcattgtgcca                   <  2:4855137‑M2/150‑89 (MQ=255)
                                                                                                                                    gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagg                                                                                 <  1:529434‑M2/87‑27 (MQ=255)
                                                                                                                                        gCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaaaagtagaaatcattgtgccattagt              <  1:755938‑M2/150‑94 (MQ=255)
                                                                                                                                         cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg             <  2:2881610/150‑1 (MQ=255)
                                                                                                                                          ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAgg            >  2:4191183/1‑150 (MQ=255)
                                                                                                                                           ccggtaGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc           >  2:2568303/6‑150 (MQ=255)
                                                                                                                                            gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgctcccaccggaaggcagtagggccccaaatttattcctcactgcgattggaggtataaccaactcgtagaaatccttgtgccagtagtcagg          <  2:499591‑M2/150‑98 (MQ=255)
                                                                                                                                            gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc          >  1:7942062/1‑150 (MQ=255)
                                                                                                                                            gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc          >  1:7372909/1‑150 (MQ=255)
                                                                                                                                            gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc          >  2:5531587/1‑150 (MQ=255)
                                                                                                                                             cgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaag                                                                                         >  2:3564052‑M2/1‑52 (MQ=255)
                                                                                                                                              gAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCAt        <  2:3933708/150‑1 (MQ=255)
                                                                                                                                               tggCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATc       >  2:5044772/3‑150 (MQ=255)
                                                                                                                                                 gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCtt     <  1:5053279/150‑1 (MQ=255)
                                                                                                                                                   aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg   >  2:2075685/1‑150 (MQ=255)
                                                                                                                                                   aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg   >  1:2068691/1‑150 (MQ=255)
                                                                                                                                                    tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaaggggcgcggggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcaggaaggcttaag  <  2:4715054‑M2/150‑106 (MQ=255)
                                                                                                                                                    tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc  >  2:3670911/1‑150 (MQ=255)
                                                                                                                                                    tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc  >  1:1088880/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCCAGAGGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGC  >  NC_002947/1316541‑1316838

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: