Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,065 T→G 19.7% intergenic (‑123/+31) PP_4560 ← / ← csbD ribonuclease BN/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,0650TG19.7% 155.0 / 32.1 81intergenic (‑123/+31)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base T (39/24);  new base G (8/7);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAC  >  NC_002947/5181920‑5182211
                                                                                                                                                   |                                                                                                                                                               
ccGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAc                                                                                                                                                              >  5:332512/1‑149 (MQ=255)
  gcgcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGAGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                            >  1:194494/1‑149 (MQ=255)
     cAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                                         <  6:332512/149‑1 (MQ=255)
       ggTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                                       <  4:74546/149‑1 (MQ=255)
       ggTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                                       <  3:115327/149‑1 (MQ=255)
       ggTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGGCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                                       <  2:263903/149‑1 (MQ=255)
        gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                                         >  4:121007/1‑146 (MQ=255)
        gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                                         >  3:67753/1‑146 (MQ=255)
         tCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaaga                                                                                                                                                     <  3:328955/149‑1 (MQ=255)
            gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                            <  2:82670/139‑1 (MQ=255)
            gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                            <  2:77148/139‑1 (MQ=255)
            gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                            >  1:77148/1‑139 (MQ=255)
            gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                            >  1:82670/1‑139 (MQ=255)
             ggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                                                 <  1:287830/149‑1 (MQ=255)
                      tAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGc                                                                                                                                        >  8:86918/1‑148 (MQ=255)
                          gTCTCCTTTCGCCGCTTGATCTTAACACCGCAACGTATAATACACCGGCCCCGCGCTGTGCTGCTCCCGGCCGTGGTTTAACAATTTGCTTACCAAAAAGGGCCGGCCTCTTCGCGGGTCCACCCTAGAAGTGGCCGGCACCGCCtcaa                                                                                                                                   <  2:9745/149‑1 (MQ=255)
                                tttCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCgg                                                                                                                                              >  6:218634/1‑132 (MQ=255)
                                tttCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCgg                                                                                                                                              <  5:218634/132‑1 (MQ=255)
                                            aTCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAACAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCAcc                                                                                                                  <  5:345190/148‑1 (MQ=255)
                                                ttACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATc                                                                                                             >  3:25062/1‑149 (MQ=255)
                                                    acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTc                                                                                                         >  5:214011/1‑149 (MQ=255)
                                                    acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCAGGCACCGCCTCAAGTCAAGGCTTCTTCACCCCATCCtt                                                                                                          >  5:208858/1‑148 (MQ=255)
                                                      aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAAATTTGCTTACGAATAAGTGCCGGCCTATTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGATTCTTCACCGCATCCTTc                                                                                                         >  1:103558/1‑148 (MQ=255)
                                                      aGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCAACCTTca                                                                                                        >  2:171003/1‑148 (MQ=255)
                                                             tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGTTTGAAGAATTTGCTTACTAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCCGGCACCGCCTCAAGTCAAGGCTTCTTCAc                                                                                                                   <  1:225781/131‑1 (MQ=255)
                                                             tCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCACGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                                                >  3:258511/1‑149 (MQ=255)
                                                               tAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                                              >  6:23330/1‑149 (MQ=255)
                                                                    gACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑Cct                                                                               >  3:283038/1‑149 (MQ=255)
                                                                       cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                                    >  4:72113/1‑88 (MQ=37)
                                                                       cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                                    <  3:72113/88‑1 (MQ=37)
                                                                       cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑Cctgct                                                                            >  7:124131/1‑149 (MQ=255)
                                                                               aGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTACTGcgcctcg                                                                  >  3:66555/1‑149 (MQ=255)
                                                                                   gtgCTGCTCCCGGCCGTGGTTGATGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCtt                                                              >  4:50948/1‑149 (MQ=255)
                                                                                    tgctgcTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCttt                                                             >  7:104055/1‑149 (MQ=255)
                                                                                     gctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGcgcctc                                                                   >  1:39173/1‑142 (MQ=255)
                                                                                     gctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGcgcctc                                                                   <  2:39173/142‑1 (MQ=255)
                                                                                      ctgctCCCGGCCGTGGTTGATGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTg                                                            <  3:50948/148‑1 (MQ=255)
                                                                                      ctgctCCCGGCCGTGGTTCAATATTTTGCTTACGAATACGGGCCGGCCTCTTCGCGGGTGCCCCCGAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGa                                                           <  4:67753/149‑1 (MQ=255)
                                                                                       tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAg                                                          >  8:203652/1‑149 (MQ=255)
                                                                                          tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                               >  3:319652/1‑73 (MQ=255)
                                                                                          tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                               <  4:319652/73‑1 (MQ=255)
                                                                                          tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTc                                                       >  7:103334/1‑149 (MQ=255)
                                                                                            ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCt                                                     >  7:209023/1‑149 (MQ=255)
                                                                                              ggCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTgc                                                   <  7:86918/149‑1 (MQ=255)
                                                                                                 cGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCCTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcc                                                >  2:283126/1‑149 (MQ=255)
                                                                                                        gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTGCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCATCGCCTTTGAGCTCCTGAgcctttcctt                                         >  7:129235/1‑149 (MQ=255)
                                                                                                        gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgcctt                                         >  7:134404/1‑149 (MQ=255)
                                                                                                          aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgcctt                                         <  4:283038/148‑1 (MQ=255)
                                                                                                            acaaTTTCCTAACGAATACGGGCCGGCCTCTCCGCGGGTCCACCCGACAAGCGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCg                                       <  8:55574/147‑1 (MQ=255)
                                                                                                                       aCGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGc                              <  1:105196/147‑1 (MQ=255)
                                                                                                                          aaTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                         <  3:121007/149‑1 (MQ=255)
                                                                                                                          aaTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                         <  6:214011/149‑1 (MQ=255)
                                                                                                                           atcctggCCGGCCTCTTCGCGGGTTCACCGTATAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                         <  6:208858/142‑1 (MQ=255)
                                                                                                                           aTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg                        >  6:150598/1‑149 (MQ=255)
                                                                                                                              aGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTg                     <  3:281633/149‑1 (MQ=255)
                                                                                                                                tGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc                   <  2:103558/149‑1 (MQ=255)
                                                                                                                                 gCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                  >  3:280100/1‑149 (MQ=255)
                                                                                                                                 gCCGGCCTCTGCGCGGGTCCCCCCCACATGTGCCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                  <  1:171003/149‑1 (MQ=255)
                                                                                                                                  ccGGCCTCTTCGCGGGTTCACCCGAGAATTGGCCGGCACCGCCTCAACTCAATGCTTCTTCACCGCATCCTTCACATTTACCCTTA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt                    <  1:283126/147‑1 (MQ=255)
                                                                                                                                  ccGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCgg                 >  3:366756/1‑149 (MQ=255)
                                                                                                                                  ccGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                  <  4:25062/148‑1 (MQ=255)
                                                                                                                                   cGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                      <  1:280321/143‑1 (MQ=255)
                                                                                                                                   cGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                      >  2:280321/1‑143 (MQ=255)
                                                                                                                                   cGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt                >  7:143919/1‑149 (MQ=255)
                                                                                                                                    ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt                >  7:345754/1‑148 (MQ=255)
                                                                                                                                    ggCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CATGCTGCGCATCGCCTTTGAGCTCCTGCGCCTTGACTTCGGCCCGCAGCTTTTCTTTGTCTGt                >  8:256212/1‑148 (MQ=255)
                                                                                                                                     gCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCGTTTCGTTGTAGGTCa              >  1:360316/1‑149 (MQ=39)
                                                                                                                                        tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTgcg                                                  >  1:320823/1‑110 (MQ=255)
                                                                                                                                        tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTgcg                                                  <  2:320823/110‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt           >  7:66302/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGGTTCACCCGAGAAGAGGCCGCCACCGCCACAAGTCAAGGCTTCTTCACCGCCTCCTTCACATTTACCATGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCTCATCGCATTAGAACTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt           >  2:186581/1‑149 (MQ=255)
                                                                                                                                         ctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTgc                                                   >  7:104402/1‑108 (MQ=38)
                                                                                                                                         ctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTgc                                                   <  8:104402/108‑1 (MQ=38)
                                                                                                                                          ttCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg         >  6:360708/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg         >  7:203200/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg         >  4:128061/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg         >  1:6700/1‑149 (MQ=255)
                                                                                                                                            cgcgGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc        <  8:103334/148‑1 (MQ=255)
                                                                                                                                             gcgGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAAGGTGAAGGGGATGATGAGCACACCACTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCgg                 >  1:122949/1‑149 (MQ=255)
                                                                                                                                                 gTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAacac  <  5:150598/149‑1 (MQ=255)
                                                                                                                                                 gTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAacac  >  4:145459/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                               
CCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGA‑‑‑‑‑‑‑‑‑‑CCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAC  >  NC_002947/5181920‑5182211

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: