Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,311,766 T→A 6.1% intergenic (‑228/+50) 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,311,7660TA 6.1% 328.5 / 8.8 109intergenic (‑228/+50)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base T (52/48);  new base A (2/5);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 114 total aligned reads displayed.

TACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGATTTC  >  NC_002947/1311617‑1311914
                                                                                                                                                      |                                                                                                                                                    
tACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCttt                                                                                                                                                      <  1:3910755/150‑1 (MQ=255)
  cTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGc                                                                                                                                                    <  1:382729/150‑1 (MQ=255)
   ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCAACCACAGCCGGGGGTGGCCCCCGCTTTGcc                                                                                                                                                   >  1:2636565/1‑150 (MQ=255)
   ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATTCCACCGGCGTTTcc                                                                                                                                                   >  2:4806724/1‑150 (MQ=255)
   ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCCTCCCCAGCGGGGGATGCCCCCCGCGTTGcc                                                                                                                                                   >  2:3130097/1‑150 (MQ=255)
    tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                                  <  1:3876513/150‑1 (MQ=255)
      aGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAATTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                                                <  1:2006119/150‑1 (MQ=255)
      aGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                                                <  2:3080794/150‑1 (MQ=255)
       gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATACCTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtg                                                                                                                                               <  1:1362005/150‑1 (MQ=255)
       gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtg                                                                                                                                               <  1:1361999/150‑1 (MQ=255)
       gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCAAATCGCCGTGAGGAATAATTTTGGGGCCCTCCTGCCTTCCTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtg                                                                                                                                               <  1:2456883/150‑1 (MQ=255)
        ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCCAACCATACACAGCGGTGGGTGCCACCCGCGGTGGCGgtgt                                                                                                                                              >  2:3412657/1‑150 (MQ=255)
          ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCCGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtggca                                                                                                                                            <  1:1983038/150‑4 (MQ=255)
           tCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCACCAACCACGGCGGGGGTTGCCACCCGCGGTGCCGGTGTTcg                                                                                                                                           >  2:522601/1‑150 (MQ=255)
            ccTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtggcatc                                                                                                                                          <  2:3175099/150‑6 (MQ=255)
             cTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTcgcg                                                                                                                                         >  2:2664898/1‑150 (MQ=255)
              tGACTATTGGCACAATGATTTCTACTATTTGTTTATAACGCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                        <  1:237143/150‑1 (MQ=255)
              tGACTAATGGCACAATTATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                        <  1:2428590/150‑1 (MQ=255)
                aCTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAACAATCCACATAGGTGGTATGCAACCGCGTTTCCGGGGTTAGCgggg                                                                                                                                      >  2:59314/1‑149 (MQ=255)
                      ggCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgccc                                                                                                                                >  1:1679751‑M2/1‑146 (MQ=255)
                       gCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCCCAGCCACCCCAGCGGTGGGTGGCCCCCGCGTTGCCGGTGTTAGCGGGCGGgcccg                                                                                                                               <  2:2749867‑M2/150‑6 (MQ=255)
                           caggattTCTACGAGTGGGTTCTAACGCCAATCGCAGTGAGGAATAACTGTGGGGCCCTAATGCCGTCCTGGGGGAGCGGGCGAGCCCGCGACGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGGGTTCGCGGGCGTgcccgctcc                                                                                                                           >  2:4717758‑M2/4‑141 (MQ=255)
                             tGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACAGGGGGGGGAGGCCCCCGCGTTTCCGGGGTTCGCGGGCGGGCCCCCGcccc                                                                                                                         >  2:3684185/1‑149 (MQ=255)
                               aTTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAGCAAACACAGCGGGGGATGCCACCCGCTTTTCCGGTGTTCGCGGGCGGgcccgcgcccaca                                                                                                                       <  2:2307597‑M2/150‑14 (MQ=255)
                                  tCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCAACCACTGCGGTGGGTGCCCCCCGCGTGGCCCGTGTTCGCGGGCGTgcccccgcccaccagg                                                                                                                    <  1:4839563‑M2/150‑17 (MQ=255)
                                   cTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAggg                                                                                                                   <  2:4467907/150‑1 (MQ=255)
                                   cTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCATACACAGCGGTGGATGCCACCGGCGGTGCCGGGGTTCGCCGGCGGGCCCGCTCCCACCAggg                                                                                                                   >  2:4750366/1‑150 (MQ=255)
                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                                  >  1:4387704/1‑118 (MQ=255)
                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                                  <  2:4387704/118‑1 (MQ=255)
                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAAACAATGCGGGGGGTGCCACCCGCTTTTCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgggg                                                                                                                  >  2:3289467/1‑150 (MQ=38)
                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGCAGCAACCACAGCGGGGGATGCCCCCCGCGTTGCCCGGGTTCGCCGGCGGGCCCCCTCCCCCAAcggg                                                                                                                  >  2:3586603/1‑150 (MQ=255)
                                         tttGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                                             <  2:4763931/150‑1 (MQ=255)
                                         tttGGTTATAACTCCAATCGCAGGGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                                             <  2:1979826/150‑1 (MQ=255)
                                            ggTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAGCAAACAAAGCGGGGGATTCCACCCGCGTTGGCCGGGTTCGCGGGCGTGCCCCCTCCCACAACGGGCGCGGGGc                                                                                                          >  2:1841965/1‑150 (MQ=255)
                                               tataACCCCAATCGCCGTGAGTAATAATTTTGGGGCCCTAAGGCCTTCGTGTGGGCGCGGGCGCGCCCGCGACGCATCCACTGCGGGGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTg                                                                                                       <  1:907913/150‑1 (MQ=25)
                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCCAAGCATACACAGCGGGGGTTGCCACCCGCGTTGCCGGTGTTCGCGGGCGTGCCCCCTCCCCAAACGGGCGCGGGGCCGGCGCTag                                                                                                 >  2:3068148/1‑150 (MQ=37)
                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCGACCAACCACACCCGGGGCTGCCCCCCGCGTTGCCCGGGGTCGCCGGCCGGGCCCCCCCCCCAACGGGCGCCGGGCCAgccccac                                                                                                 >  2:1444373/1‑145 (MQ=255)
                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAACAAACCCAGCGGGGGGTTGCACCCCCTTTTCCGGGGGTTGCGGGCGGGCCCCCGCCCCCAACGGGCGCCGGGCCAgcccaag                                                                                                 >  1:432941/1‑145 (MQ=255)
                                                      ccAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCATACACAGCGGTGGATGCCACCCGCGTTGCCCGTGTTAGCGGGCGTGCCCCCTCCCACAAGGGGCGCGGGGCCAGCGCTagc                                                                                                >  2:1280811/1‑150 (MQ=255)
                                                        aaTCGCAGTGAGGAATAAATTTGGGGCCCTCATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                                              <  2:4224767/150‑1 (MQ=255)
                                                            gCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGCGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAgccagcca                                                                                          <  2:4347960/150‑1 (MQ=255)
                                                               gTGAGGAATAATTTTGGGGCCCTACGGCCTTCGTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccacca                                                                                       <  1:3220572/150‑1 (MQ=255)
                                                                       tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccc                                                                                                                          >  1:353722‑M2/1‑97 (MQ=255)
                                                                       tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccc                                                                                                                          <  2:353722‑M2/107‑11 (MQ=255)
                                                                         attttGGGGCCCTTTTGCCTTCGTGTGGGAGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                             <  1:4750366/148‑1 (MQ=255)
                                                                         aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                             <  1:2822558/150‑1 (MQ=255)
                                                                         aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                             <  2:5129936/150‑1 (MQ=255)
                                                                         aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                             <  2:2117132/150‑1 (MQ=255)
                                                                                            ttCCTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCGAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                                          <  1:2740760/150‑1 (MQ=255)
                                                                                            ttCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                                          <  2:1996645/150‑1 (MQ=255)
                                                                                            ttCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                                          <  2:3936505/150‑1 (MQ=255)
                                                                                              cgtgtgGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                                        <  1:1354479/148‑1 (MQ=255)
                                                                                              ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                                        <  1:4267959/150‑1 (MQ=255)
                                                                                              ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTGGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                                        <  1:4268126/150‑1 (MQ=255)
                                                                                                  tgGGAGCGGGCGAGCCCGCCAAGCATACACTGCGGTGGATGCCCCCGGCATTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCCGGGCCAGCCCAAGCCACCCACCAACACGCGAATTGCATACAGGCGCCGGGCGGCCGtc                                                    >  2:2235594/1‑150 (MQ=25)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                                           >  2:2952054/1‑150 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCGTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                                           >  2:4696499/1‑150 (MQ=255)
                                                                                                              aGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattag                                        <  1:3185955‑M2/150‑93 (MQ=255)
                                                                                                                cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                      >  2:3748696/1‑150 (MQ=255)
                                                                                                                cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                      >  2:1306251/1‑150 (MQ=255)
                                                                                                                 ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaaggggcgcggggcctcaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcagg                                     <  1:381620‑M2/150‑96 (MQ=255)
                                                                                                                 ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTGCCAGGCGCAGGGCGGCCTTCTCCCGCATTGCCAgg                                     >  2:49804/1‑150 (MQ=255)
                                                                                                                  cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc                                    >  2:3480590/1‑150 (MQ=255)
                                                                                                                  cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc                                    <  1:3480590/150‑1 (MQ=255)
                                                                                                                       aGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcaggaagg                               <  1:4995665‑M2/150‑102 (MQ=255)
                                                                                                                        gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCa                                                >  1:4746595/1‑132 (MQ=255)
                                                                                                                        gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCa                                                <  2:4746595/132‑1 (MQ=255)
                                                                                                                        gCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCtt                              >  1:4219828/1‑150 (MQ=255)
                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcc                                                                                               <  1:2390496/84‑1 (MQ=255)
                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcc                                                                                               >  2:2390496/1‑84 (MQ=255)
                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTc                             <  2:2561006/150‑1 (MQ=255)
                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTc                             <  2:5127644/150‑1 (MQ=255)
                                                                                                                          aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg                            >  2:1792957/1‑150 (MQ=255)
                                                                                                                            ccACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgcg                          >  1:841259/1‑150 (MQ=255)
                                                                                                                             cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                                              >  1:4359944/1‑65 (MQ=255)
                                                                                                                             cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                                              <  2:4359932/65‑1 (MQ=255)
                                                                                                                             cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                                              <  2:4359944/65‑1 (MQ=255)
                                                                                                                             cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                                              >  1:4359932/1‑65 (MQ=255)
                                                                                                                                tGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcaggaaggcttaagtac                      <  1:5377029‑M2/150‑111 (MQ=255)
                                                                                                                                  cGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCtt                    >  1:1657228/1‑150 (MQ=255)
                                                                                                                                      ggATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTgcgc                <  1:522601/150‑1 (MQ=255)
                                                                                                                                       gATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCa               >  2:556812/1‑150 (MQ=255)
                                                                                                                                       gATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCa               >  2:556431/1‑150 (MQ=255)
                                                                                                                                       gATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCa               >  1:377270/1‑150 (MQ=255)
                                                                                                                                           aacACCGGCAAAGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAggg            >  2:254821/3‑150 (MQ=255)
                                                                                                                                             cACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCg          >  2:1230430/1‑150 (MQ=255)
                                                                                                                                             cACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCg          >  2:1260560/1‑150 (MQ=255)
                                                                                                                                              aCCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGt         >  1:643484/1‑150 (MQ=255)
                                                                                                                                              aCCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGt         >  2:1680141/1‑150 (MQ=255)
                                                                                                                                              aCCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGt         >  2:582316/1‑150 (MQ=255)
                                                                                                                                              aCCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCGAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCCGGGCGGCCTTCTCCAGCATTGCCAGGCCCTCTTCGCGCCGCTTGCGCAGGGCGt         >  1:643425/1‑150 (MQ=255)
                                                                                                                                               ccGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTc        >  2:4426911/1‑150 (MQ=255)
                                                                                                                                               ccGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTc        >  1:2614968/1‑150 (MQ=255)
                                                                                                                                               ccGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTc        <  1:2609646/150‑1 (MQ=255)
                                                                                                                                                cGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCg       >  2:4001302/1‑150 (MQ=255)
                                                                                                                                                  gCTTTGCCGGTGTTCGCGGGCGTGCCCGCTACCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGAt     >  1:5222392/1‑150 (MQ=255)
                                                                                                                                                    tttGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGAttt   >  1:3786695/1‑150 (MQ=255)
                                                                                                                                                    tttGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGAttt   <  1:1306251/150‑1 (MQ=255)
                                                                                                                                                    tttGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTAGCGCCGCTTGCGCAGGGCGTCGAttt   >  1:3786089/1‑150 (MQ=255)
                                                                                                                                                     ttGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGATTTc  <  1:2749867/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGATTTC  >  NC_002947/1311617‑1311914

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