Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,312,039 G→A 5.4% L184L (CTC→CTT csiR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,312,0390GA 5.4% 987.0 / 28.8 312L184L (CTC→CTTcsiRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (147/147);  new base A (2/15);  total (149/163)
Fisher's exact test for biased strand distribution p-value = 2.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 317 total aligned reads displayed.

TCATCTAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCTTGTGCCGCTC  >  NC_002947/3311891‑3312188
                                                                                                                                                    |                                                                                                                                                     
  aTCTAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagc                                                                                                                                                    <  2:70760/150‑1 (MQ=255)
     tAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcc                                                                                                                                                 >  2:1242519/1‑150 (MQ=255)
     tAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcc                                                                                                                                                 >  2:2768034/1‑150 (MQ=255)
     tAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcc                                                                                                                                                 >  2:4072852/1‑150 (MQ=255)
     tAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcc                                                                                                                                                 >  1:2909178/1‑150 (MQ=255)
      aGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcca                                                                                                                                                <  1:6063496/150‑1 (MQ=255)
      aGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcca                                                                                                                                                <  2:4511762/150‑1 (MQ=255)
      aGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgcca                                                                                                                                                <  2:5365019/150‑1 (MQ=255)
       gCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTCAagcgccag                                                                                                                                               >  2:1783597/1‑150 (MQ=255)
        cGGGCACCGATGCCTTCGGCGTGCATGATCTGCACAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgccagc                                                                                                                                              <  1:6751127/150‑1 (MQ=255)
         gggCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCTGCATGCTCCTGGCGCTCAagcgccagga                                                                                                                                             <  1:6186509/150‑3 (MQ=255)
         gggCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgccagcg                                                                                                                                             >  1:3414224/1‑150 (MQ=255)
         gggCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTCAagcgccagga                                                                                                                                             <  1:799145/150‑3 (MQ=255)
         gggCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTCAagcgccagga                                                                                                                                             <  2:5395637/150‑3 (MQ=255)
         gggCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTCAagcgccagga                                                                                                                                             <  1:3428894/150‑3 (MQ=255)
          ggCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGagcgccagcgc                                                                                                                                            <  2:3585576/150‑1 (MQ=255)
               cGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCtt                                                                                                                                       >  1:466311/1‑150 (MQ=255)
                gATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGTTTTTGGCATCCCGGGCGTGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTc                                                                                                                                      <  2:335987/150‑1 (MQ=255)
                gATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTc                                                                                                                                      <  2:186601/150‑1 (MQ=255)
                    ccTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGaa                                                                                                                                  >  1:3508820/1‑150 (MQ=255)
                     cTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGaaa                                                                                                                                 >  2:5599131/1‑150 (MQ=255)
                     cTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGaaa                                                                                                                                 >  2:6482668/1‑150 (MQ=255)
                              gCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCt                                                                                                                        >  2:6673350/1‑150 (MQ=255)
                                   aTCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTc                                                                                                                   >  2:4683832/1‑150 (MQ=255)
                                    tCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCa                                                                                                                  <  2:5039967/150‑1 (MQ=255)
                                    tCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCa                                                                                                                  <  2:6319270/150‑1 (MQ=255)
                                      tGCGCATTTATGGGTACGGGTGTCATCTGGTGCGAGCGCATCCGGGCGCTGGCGTTTTTGGCATCCCGGGCGAGTATCACCTCGACCTGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGc                                                                                                                <  1:1424778/150‑1 (MQ=255)
                                           actGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCCTGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAg                                                                                                           <  1:535057/148‑1 (MQ=255)
                                           aaTGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTTCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAg                                                                                                           >  1:5884678/1‑150 (MQ=255)
                                             tGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGt                                                                                                         >  2:4457125/1‑150 (MQ=255)
                                               aTCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTgg                                                                                                       >  2:1483642/1‑150 (MQ=255)
                                                tCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGc                                                                                                      >  1:3886063/1‑150 (MQ=255)
                                                 cGGTACCGGTGTCATCTGGTGCGAGCGCATCATGGCGGTGGCGGTTTGGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCg                                                                                                     <  1:4821784/150‑1 (MQ=255)
                                                   gTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGt                                                                                                   >  2:6625592/1‑150 (MQ=255)
                                                    tACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTa                                                                                                  <  1:1503826/150‑1 (MQ=255)
                                                    tACCGGGGTCATCCGGTGCGAGCGCAGCAGGGCGCGGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCGGGCGCTTGCGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTa                                                                                                  <  1:4185986/150‑1 (MQ=255)
                                                      ccGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAgc                                                                                                >  1:2816235/1‑150 (MQ=255)
                                                           ttcatcaGGTGCGAGCGAATCATGGCGCTGGCGTTTTTGGCATCCCGGGCGAGTATAACCTCGACCAGTGCCGAATGCTCCTGGCGCTTTAGCGCAAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCg                                                                                           <  1:2275279/149‑1 (MQ=255)
                                                            tcatcaGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCgg                                                                                          >  1:125492/1‑150 (MQ=255)
                                                               tcaGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCt                                                                                       >  2:2026424/1‑150 (MQ=255)
                                                               tcaGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCt                                                                                       >  2:5854264/1‑150 (MQ=255)
                                                                   gTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTc                                                                                   >  2:149921/1‑150 (MQ=255)
                                                                   gTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCGGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTc                                                                                   <  1:2459326/150‑1 (MQ=255)
                                                                     gCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGa                                                                                 >  1:598335/1‑150 (MQ=255)
                                                                       gAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAAc                                                                               >  1:6471824/1‑150 (MQ=255)
                                                                        aGCGCACCAGGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACa                                                                              <  1:4575339/150‑1 (MQ=255)
                                                                           gcATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACaggt                                                                           <  2:4726662/150‑1 (MQ=255)
                                                                              tcatGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACaggtagg                                                                        >  2:3910951/1‑150 (MQ=255)
                                                                                   gcgcTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGgtgcgt                                                                   <  2:6442243/150‑1 (MQ=255)
                                                                                     gcTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGgtgcgtgc                                                                 >  2:1156607/1‑150 (MQ=255)
                                                                                     gcTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGgtgcgtgc                                                                 <  1:4316366/150‑1 (MQ=255)
                                                                                        gggggTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTg                                                              <  1:1790882/147‑1 (MQ=255)
                                                                                           ggTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTgcag                                                           >  1:4025189/1‑150 (MQ=255)
                                                                                             ttttGGCATCCCGGGCGAGGATCACCTCGACCTGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTgcagca                                                         <  1:5559787/150‑1 (MQ=255)
                                                                                             ttttGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTgcagca                                                         >  1:2117052/1‑150 (MQ=255)
                                                                                             ttttGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTgcagca                                                         >  2:2182154/1‑150 (MQ=255)
                                                                                              tttGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGcagcag                                                        <  1:6374882/150‑1 (MQ=255)
                                                                                              ttaggCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCAAGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGcagcag                                                        >  2:4416112/4‑150 (MQ=255)
                                                                                                  gCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGc                                                    >  2:490809/1‑150 (MQ=255)
                                                                                                   cATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCt                                                   >  1:6827402/1‑150 (MQ=255)
                                                                                                    aTCCCGGGCGAGTATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCtt                                                  <  2:1690131/150‑1 (MQ=255)
                                                                                                     tCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTg                                                 >  1:5155747/1‑150 (MQ=255)
                                                                                                     tCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTg                                                 <  1:2694320/150‑1 (MQ=255)
                                                                                                       ccGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCg                                                                                  <  1:1730674/115‑1 (MQ=255)
                                                                                                       ccGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGa                                               <  1:3293576/150‑1 (MQ=255)
                                                                                                        cGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAg                                              <  2:6215818/150‑1 (MQ=255)
                                                                                                        cGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAg                                              <  2:4162899/150‑1 (MQ=255)
                                                                                                               ggATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGc                                       >  1:6590616/1‑150 (MQ=255)
                                                                                                                gATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGTAGCAGGTGCTTGGAGCCGCAGcc                                      >  1:4911852/1‑150 (MQ=255)
                                                                                                                  tCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCg                                    >  2:4832632/1‑150 (MQ=255)
                                                                                                                    aCCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAg                                  <  1:4283413/150‑1 (MQ=255)
                                                                                                                    aCCTCGACCAGTGCCGCAGGCGCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAg                                  <  2:2884748/150‑1 (MQ=255)
                                                                                                                     ccTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAgg                                 <  2:5568670/150‑1 (MQ=255)
                                                                                                                     ccTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAgg                                 >  1:3760733/1‑150 (MQ=255)
                                                                                                                     ccTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAgg                                 <  2:812062/150‑1 (MQ=255)
                                                                                                                      cTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGc                                >  1:4504472/1‑150 (MQ=255)
                                                                                                                      cTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGc                                >  2:2499644/1‑150 (MQ=255)
                                                                                                                       tCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCa                               >  2:2399268/1‑150 (MQ=255)
                                                                                                                         gACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACag                                                                             >  1:5281724/1‑102 (MQ=255)
                                                                                                                         gACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAAt                             <  2:4221199/150‑1 (MQ=255)
                                                                                                                         gACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAAt                             >  1:543389/1‑150 (MQ=255)
                                                                                                                          aCCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATg                            <  1:4072852/150‑1 (MQ=255)
                                                                                                                             aGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCa                         <  1:4862924/150‑1 (MQ=255)
                                                                                                                             aGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCa                         <  2:6380004/150‑1 (MQ=255)
                                                                                                                              gTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAg                        >  2:5484521/1‑150 (MQ=255)
                                                                                                                                 ccGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTAt                     >  2:1984939/1‑150 (MQ=255)
                                                                                                                                 ccGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTAt                     <  2:5156095/150‑1 (MQ=255)
                                                                                                                                  cGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCcac                                                                                                             <  2:2267308/61‑1 (MQ=255)
                                                                                                                                  cGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACaggta                                                                          >  2:1085437/1‑96 (MQ=255)
                                                                                                                                  cGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACaggta                                                                          <  1:1085437/96‑1 (MQ=255)
                                                                                                                                   gCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATgg                   <  1:4983970/150‑1 (MQ=255)
                                                                                                                                    cATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGa                  >  1:1324623/1‑150 (MQ=255)
                                                                                                                                        gcgcTGGCGCTCAAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGc              >  2:3774746/4‑150 (MQ=255)
                                                                                                                                         tCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGcc             >  2:4377672/1‑150 (MQ=255)
                                                                                                                                           cTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCtt           >  2:5931279/1‑150 (MQ=255)
                                                                                                                                           cTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCtt           <  1:717248/150‑1 (MQ=255)
                                                                                                                                            tGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAgg                                                       >  2:6756744/1‑105 (MQ=255)
                                                                                                                                            tGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCTtg          <  2:3177684/150‑1 (MQ=255)
                                                                                                                                             ggCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCTtgt         <  2:4104578/150‑1 (MQ=255)
                                                                                                                                                    gAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCTTGTGCCGCTc  >  2:2693001/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCATCTAGCGGGCACCGATGCCTTCGGCGTGCATGATCTGCGCAATGATCGGTACCGGTGTCATCAGGTGCGAGCGCATCATGGCGCTGGCGGTTTTGGCATCCCGGGCGAGGATCACCTCGACCAGTGCCGCATGCTCCTGGCGCTTGAGCGCCAGCGCCTCTTCGGAAAACACCGTCTGCGTCAGCCACAGGTGGCGGTAGCGCTCGGCCTGGTCGAACAGGTAGGTGCGTGCCTGCAGCAGGTGCTTGGAGCCGCAGCCCGAGGCAATGGCAGTATGGAACGCCTTGTGCCGCTC  >  NC_002947/3311891‑3312188

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