Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 615,604 G→C 7.5% intergenic (‑119/+47) xseB ← / ← ribB exodeoxyribonuclease 7 small subunit/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947615,6040GC 7.5% 161.7 / 4.6 53intergenic (‑119/+47)xseB/ribBexodeoxyribonuclease 7 small subunit/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base C (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

ATCGAGGGAGGCTTTTTTGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAATTGTACTGCCGTGCATATACCGCCA  >  NC_002947/615459‑615751
                                                                                                                                                 |                                                                                                                                                   
aTCGAGGGAGGCTTTTTTGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGTAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGaaa                                                                                                                                                  <  2:76599/149‑1 (MQ=255)
   gagggaggCTTTTTTGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAgg                                                                                                                                                <  8:252600/148‑1 (MQ=255)
    agggaggCTTTTTTGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGATGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGc                                                                                                                                              >  1:25943/1‑149 (MQ=255)
          gCTTTTTTGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCg                                                                                                                                        >  3:240365/1‑149 (MQ=255)
            ttttttGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGca                                                                                                                                      >  1:2914/1‑149 (MQ=255)
            ttttttGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGca                                                                                                                                      >  6:212940/1‑149 (MQ=255)
              ttttGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGca                                                                                                                                      <  8:1969/147‑1 (MQ=255)
                   cGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGa                                                                                                                                <  5:212940/148‑1 (MQ=255)
                    gggCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGAcc                                                                                                                              >  6:236697/1‑149 (MQ=255)
                       ccATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGc                                                                                                                            >  7:73923/1‑148 (MQ=255)
                         aTGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGCGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCaaa                                                                                                                         <  2:2914/149‑1 (MQ=255)
                                         tGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTTTCGCGACACAGGGCCGCt                                                                                                                                          >  5:109382/1‑116 (MQ=38)
                                         tGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTTTCGCGACACAGGGCCGCt                                                                                                                                          <  6:109382/116‑1 (MQ=38)
                                            cATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGAt                                                                                                      >  1:119610/1‑149 (MQ=255)
                                                 aGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGc                                                                                                 >  1:81434/1‑149 (MQ=255)
                                                           gACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGtt                                                                                                              <  6:20393/126‑1 (MQ=255)
                                                           gACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGtt                                                                                                              >  5:20393/1‑126 (MQ=255)
                                                               cTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTccc                                                                                   <  5:219197/149‑1 (MQ=255)
                                                                tAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGACAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCg                                                                                  >  8:83128/1‑149 (MQ=255)
                                                                  gCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGtt                                                                                >  6:124939/1‑149 (MQ=255)
                                                                       gcggggctggGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGc                                                                           <  2:147341/149‑1 (MQ=255)
                                                                             cggggGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGcc                                                                     <  5:124939/149‑1 (MQ=255)
                                                                                 gggAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCAc                                                                 <  8:73923/149‑1 (MQ=255)
                                                                                 gggAGCAAACTTCCGACAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCAc                                                                 <  7:83128/149‑1 (MQ=255)
                                                                                     gCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTccc                                                                                                                  <  4:54653/96‑1 (MQ=255)
                                                                                     gCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTccc                                                                                                                  >  3:54653/1‑96 (MQ=255)
                                                                                       aaaCTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTAcc                                                           >  3:208675/1‑149 (MQ=255)
                                                                                       aaaCTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTAcc                                                           >  5:232664/1‑149 (MQ=255)
                                                                                           ttCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCgcgc                                                       >  2:189881/1‑149 (MQ=255)
                                                                                           ttCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCgcgc                                                       >  5:233939/1‑149 (MQ=255)
                                                                                                   tcTGCATCAGCCGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGctccg                                               >  2:53195/2‑148 (MQ=255)
                                                                                                      gCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCgcgc                                                       <  3:129785/138‑1 (MQ=255)
                                                                                                      gCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGa                                            >  1:125677/1‑149 (MQ=255)
                                                                                                      gCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGCCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCgcgc                                                       >  4:129785/1‑138 (MQ=255)
                                                                                                      gCATCAGCCGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGc                                                   >  8:232506/1‑142 (MQ=255)
                                                                                                      gCATCAGCCGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGc                                                   <  7:232506/142‑1 (MQ=255)
                                                                                                       cATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGAt                                           <  5:236697/149‑1 (MQ=255)
                                                                                                       cATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGAt                                           <  2:25943/149‑1 (MQ=255)
                                                                                                           ggCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGt                                        >  1:179934/1‑148 (MQ=255)
                                                                                                           ggCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGt                                        <  2:119610/148‑1 (MQ=255)
                                                                                                            gCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCAtt                                                                        <  6:161944/115‑1 (MQ=255)
                                                                                                            gCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCAtt                                                                        >  5:161944/1‑115 (MQ=255)
                                                                                                                   cAGACCGATGGGGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACAGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCgg                               >  1:205265/1‑149 (MQ=255)
                                                                                                                    aGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGc                              <  6:233939/149‑1 (MQ=255)
                                                                                                                        cGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGctcct                                               >  3:176670/1‑128 (MQ=255)
                                                                                                                        cGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGctcct                                               <  4:176670/128‑1 (MQ=255)
                                                                                                                         gATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGGCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAAtt                         <  1:189881/149‑1 (MQ=255)
                                                                                                                         gATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAAtt                         <  2:125677/149‑1 (MQ=255)
                                                                                                                             gaggagCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAATTGTa                      <  1:53195/148‑1 (MQ=255)
                                                                                                                                gagCGGCCCTTTCGCGACACAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCAc                                                                 <  4:187472/102‑1 (MQ=38)
                                                                                                                                gagCGGCCCTTTCGCGACACAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCAc                                                                 >  3:187472/1‑102 (MQ=38)
                                                                                                                                gagCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAATTGTACTGc                  >  3:203044/1‑149 (MQ=255)
                                                                                                                                  gCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCAGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAATTGTACTGCCg                >  8:53825/1‑149 (MQ=255)
                                                                                                                                                cgAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAATTGTACTGCCGTGCATATACCGCCa  >  6:125478/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATCGAGGGAGGCTTTTTTGCGGGCCATGGCGACTGTCTCGGTGGCATTCAGAACGGCGCGACACTAGCAGAGCGGGGCGGGGGGAGCAAACTTCCGGCAGCTGCATCGGCTGATGCAGACCGATGGAGGAGCGGCCCTGTGTCGCGAAAGGGCCGCTCGCACAAGGACCGCAAAGCGTCCCGGTTCGGTGGATCAGGCCGGCTCTGCCTCCCGTTCCAGCATTGCCTCACGGTACCGCGCCAGCTCCTCGATGGTCAGCACCGGCAAATTGTACTGCCGTGCATATACCGCCA  >  NC_002947/615459‑615751

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