Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,313,636 G→C 5.3% R341G (CGC→GGC)  rapA ← RNA polymerase‑associated protein RapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,313,6360GC 5.3% 232.0 / 3.4 75R341G (CGC→GGC) rapARNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

TCGCCGTCGCTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAGGCCGTATTCGGC  >  NC_002947/1313488‑1313783
                                                                                                                                                    |                                                                                                                                                   
tCGCCGTCGCTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGcgcg                                                                                                                                                     <  1:91813/149‑1 (MQ=255)
  gCCGTCGCTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCgg                                                                                                                                                    <  1:193855/148‑1 (MQ=255)
       cGCTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCACAGCATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGACGAGCAGTTCCTGCACCGCCTCGGCTACAGGGCCGAACTGCACGCTTTCGGCGCGGAAAGCg                                                                                                                                              >  2:386268/1‑149 (MQ=255)
         cTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGc                                                                                                                                            >  5:382092/1‑149 (MQ=255)
         cTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCCGCGCCGAAAGCgag                                                                                                                                            <  1:156590/149‑3 (MQ=255)
         cTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCCGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGc                                                                                                                                            >  7:220332/1‑149 (MQ=255)
          tGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCCAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGcc                                                                                                                                           >  1:344700/1‑149 (MQ=255)
                cTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCg                                                                                                                                     >  1:276145/1‑149 (MQ=255)
                   ccAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTg                                                                                                                                   >  5:20484/1‑148 (MQ=255)
                     agcagGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGaa                                                                                                                                 >  6:129272/1‑148 (MQ=255)
                          ggCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGt                                                                                                                           <  7:210402/149‑1 (MQ=255)
                           gCTTCGCCGTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGtt                                                                                                                          <  6:219383/149‑1 (MQ=255)
                            cTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTg                                                                                                                         <  3:205918/149‑1 (MQ=255)
                                  cTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTc                                                                                                                    >  2:91840/1‑148 (MQ=255)
                                     cGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGa                                                                                                                  <  8:301190/147‑1 (MQ=255)
                                     cGGCACCCAGGAAGCCCAGGAGGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGa                                                                                                                  <  6:70938/147‑1 (MQ=255)
                                      ggCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAg                                                                                                                 >  3:46663/1‑147 (MQ=255)
                                        cACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGGGAAAACGGTTGGGGTCGAGTAgg                                                                                                             >  1:156722/1‑149 (MQ=255)
                                         aCCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGtt                                                                                                                          >  4:315920/1‑135 (MQ=255)
                                         aCCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGtt                                                                                                                          <  3:315920/135‑1 (MQ=255)
                                         aCCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGgc                                                                                                            >  4:85818/1‑149 (MQ=255)
                                         aCCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGgc                                                                                                            >  4:151351/1‑149 (MQ=255)
                                          cccAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGa                                                                                                                  <  4:304243/142‑1 (MQ=255)
                                          cccAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGa                                                                                                                  >  3:304243/1‑142 (MQ=255)
                                            cAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGgcg                                                                                                           >  2:281129/1‑147 (MQ=255)
                                                aaGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCTCCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTTGGGTCGAGTAGGCGCAggcg                                                                                                     <  8:219390/149‑1 (MQ=255)
                                                   cccAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTTGGCGCAggcgggc                                                                                                  <  8:317682/149‑1 (MQ=255)
                                                    ccAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAggcgggc                                                                                                  >  7:49761/1‑148 (MQ=255)
                                                          tggtggCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTgg                                                                                           <  5:129272/149‑1 (MQ=255)
                                                             tggCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTg                                                                                        >  8:170935/1‑149 (MQ=255)
                                                                cGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGATAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTcc                                                                                     <  1:386268/149‑1 (MQ=255)
                                                                 gTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCt                                                                                    <  6:382092/149‑1 (MQ=255)
                                                                      ccTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTggg                                                                                                                       >  5:108955/1‑109 (MQ=255)
                                                                      ccTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTggg                                                                                                                       <  6:108955/109‑1 (MQ=255)
                                                                       cTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGcc                                                                                >  6:299599/1‑147 (MQ=255)
                                                                       cTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCa                                                                              <  7:290949/149‑1 (MQ=255)
                                                                             cGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCCGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAgctgct                                                                        <  8:220332/149‑1 (MQ=255)
                                                                                cccaaggCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGGCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCgg                                                                     <  6:313317/144‑1 (MQ=255)
                                                                                cAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCACCGCTACCGGGCGGTACTCCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCCCAGGCGGGCGAAGTCGCTGTCCTGGCCCAGCAGCTCgg                                                                     >  4:360660/1‑149 (MQ=255)
                                                                                 aGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCggg                                                                    >  4:167597/1‑149 (MQ=255)
                                                                                  gACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCgggg                                                                   <  1:281129/149‑1 (MQ=255)
                                                                                   aCGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCgg                                                                     >  3:295599/1‑146 (MQ=255)
                                                                                   aCGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCgg                                                                     <  4:295599/146‑1 (MQ=255)
                                                                                     ggCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTg                                                                <  2:276145/149‑1 (MQ=255)
                                                                                        ccTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGGGGTAGTGCTCGGTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAgc                                                                            <  4:213286/134‑1 (MQ=255)
                                                                                        ccTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCg                                                             >  5:824/1‑149 (MQ=255)
                                                                                        ccTCGTCGAGCAGTTCCTGCACCGCCTCGCCAACCGGGAGGTAGTCCTCGCTTTCGGCGCGGAAACCGGCCACGACGTCAAAACGGTTGGGGTCCAGAAGGCGCAGGCGCGCGAACTGGCTGTCCAGTACCAgc                                                                            >  3:213286/1‑134 (MQ=255)
                                                                                        ccTCGTCGAGCAGGTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCg                                                             <  4:46663/149‑1 (MQ=255)
                                                                                            gtcgAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCTGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCa                                                         <  2:344700/149‑1 (MQ=255)
                                                                                              cgAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTc                                                                                                                    <  3:167653/88‑1 (MQ=255)
                                                                                              cgAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTc                                                                                                                    >  4:167653/1‑88 (MQ=255)
                                                                                                  cAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCaa                                                     >  5:78014/1‑147 (MQ=255)
                                                                                                  cAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACa                                                   <  1:128606/149‑1 (MQ=255)
                                                                                                       ccTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACg                                              >  3:171606/1‑149 (MQ=255)
                                                                                                         tGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGcc                                            >  1:296766/1‑149 (MQ=255)
                                                                                                         tGCACCGCCTCGGCTACCGGGCGGTAGTGCTAGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGcc                                            >  1:191752/1‑149 (MQ=255)
                                                                                                          gCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGccc                                           >  1:317508/1‑149 (MQ=255)
                                                                                                           cACCGCCTCGGCTACCGGGCTGTAGTGCTCGCTTTCGGCGCGGAAAGCTGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGTTCAGCAACAGTACGccc                                           <  2:191752/148‑1 (MQ=255)
                                                                                                           cACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTgg                                                               <  7:275298/128‑1 (MQ=255)
                                                                                                           cACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTgg                                                               >  8:275298/1‑128 (MQ=255)
                                                                                                            aCCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCggggg                                                                  <  8:125270/124‑1 (MQ=255)
                                                                                                            aCCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCggggg                                                                  >  7:125270/1‑124 (MQ=255)
                                                                                                             ccGCCTCGGCTACCTGGCGGTAGTGCTCGCTTTCGGCGCGGAATGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGa                                        <  3:360660/149‑1 (MQ=255)
                                                                                                             ccGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGa                                        <  3:167597/149‑1 (MQ=255)
                                                                                                                      cTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTgg                               <  4:171606/149‑1 (MQ=255)
                                                                                                                      cTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTgg                               <  5:299599/149‑1 (MQ=255)
                                                                                                                            ggCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGt                         >  6:200321/1‑149 (MQ=255)
                                                                                                                                 tAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTc                    <  3:151351/149‑1 (MQ=255)
                                                                                                                                  aGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTc                    <  2:156722/148‑1 (MQ=255)
                                                                                                                                   gTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCg                   <  8:49761/148‑1 (MQ=255)
                                                                                                                                    tGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGAc                 <  2:317508/149‑1 (MQ=255)
                                                                                                                                       tCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAg              >  2:251397/1‑149 (MQ=255)
                                                                                                                                       ccGCTTTCCGCGCCGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAg              >  8:284349/2‑149 (MQ=255)
                                                                                                                                        cGCTTTCCGCGCCGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACa                                                   <  4:43974/111‑1 (MQ=255)
                                                                                                                                        cGCTTTCCGCGCCGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACa                                                   >  3:43974/1‑111 (MQ=255)
                                                                                                                                         gCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTg                                <  1:161409/129‑1 (MQ=255)
                                                                                                                                          cTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAGGcc           <  1:251397/149‑1 (MQ=255)
                                                                                                                                                 cgcgGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGGCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAGGCCGTATTCg    <  2:49920/149‑1 (MQ=255)
                                                                                                                                                   cgGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAGGCCGTATTCGgc  <  6:78014/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGCCGTCGCTGACCGCTGCCAGCAGGGCTTCGCCTTCGGCACCCAGGAAGCCCAGGATGGTGGCGTGGGCCTTGGGCGACAGACGGCCCTCGTCGAGCAGTTCCTGCACCGCCTCGGCTACCGGGCGGTAGTGCTCGCTTTCGGCGCGGAAAGCGGCCAGGTCGTGAAAACGGTTGGGGTCGAGTAGGCGCAGGCGGGCGAAGTGGCTGTCCTGGCCCAGCTGCTCGGGGGTGGCGGTCAGCAACAGTACGCCCGGAATCACTTGGGCCAGTTGCTCGACCAGGCCGTATTCGGC  >  NC_002947/1313488‑1313783

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