Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,024,413 C→T intergenic (‑65/‑58) yccF ← / → helD conserved inner membrane protein YccF/DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,024,4130CT94.6% 111.8 / ‑1.7 37intergenic (‑65/‑58)yccF/helDconserved inner membrane protein YccF/DNA helicase IV
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (13/22);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGC  >  NC_000913/1024276‑1024544
                                                                                                                                         |                                                                                                                                   
tCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTAc                                                                                                                                   >  2:581151/1‑140 (MQ=255)
   ccAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGtt                                                                                                                                <  1:192185/140‑1 (MQ=255)
   ccAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGtt                                                                                                                                <  1:662237/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:190489/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:865168/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:726674/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:490238/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:40560/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:155853/140‑1 (MQ=255)
         gccagccCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  1:281565/140‑1 (MQ=255)
                   agTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTaca                                                                                                                <  1:91725/140‑1 (MQ=255)
                    gtggtgGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          >  1:195892/1‑129 (MQ=255)
                    gtggtgGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTctct                                                                                                                          <  2:195892/129‑1 (MQ=255)
                    gtggtgGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTacac                                                                                                               <  1:536477/140‑1 (MQ=255)
                        tggCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTaca                                                                                                                <  2:697560/135‑1 (MQ=255)
                        tggCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTaca                                                                                                                >  1:697560/1‑135 (MQ=255)
                        tggCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGAt                                                                                                           >  1:6463/1‑140 (MQ=255)
                         ggCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATa                                                                                                          <  2:667525/140‑1 (MQ=255)
                                        cacaAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCgtg                                                                                           <  2:626471/140‑1 (MQ=255)
                                                  aGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACa                                                                                 >  1:571/1‑140 (MQ=255)
                                                  aGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACa                                                                                 >  1:824316/1‑140 (MQ=255)
                                                         tCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCAt                                                                          >  2:777741/1‑140 (MQ=255)
                                                             aaCGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGaa                                                                      <  2:238925/140‑1 (MQ=255)
                                                                ggTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTg                                                                   >  1:406391/1‑140 (MQ=255)
                                                                  tACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGaa                                                                 >  1:508933/1‑140 (MQ=255)
                                                                   aCGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGaaa                                                                >  1:695023/1‑140 (MQ=255)
                                                                          aaaCCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACaa                                                         <  2:824316/140‑1 (MQ=255)
                                                                              cTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCt                                                     >  1:379626/1‑140 (MQ=255)
                                                                                     tACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGaaaa                                              >  2:718624/1‑140 (MQ=255)
                                                                                       cTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTAGAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACg                                            >  2:726286/1‑140 (MQ=255)
                                                                                              tttCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGca                                     <  1:644635/140‑1 (MQ=255)
                                                                                                 cAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAg                                  <  2:6463/140‑1 (MQ=255)
                                                                                                  aGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGc                                 <  1:777741/140‑1 (MQ=255)
                                                                                                         ggcggcAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAg                                  >  2:245285/1‑132 (MQ=255)
                                                                                                         ggcggcAATTGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAg                                  <  1:245285/132‑1 (MQ=255)
                                                                                                                 ttGTAACGTTTTTTTGGGCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCggg                  <  2:379626/140‑1 (MQ=255)
                                                                                                                                 ggCTGGAGTACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGc  <  2:508933/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
TCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGTTACACTGATAAGAAATAATTTCGTGTGGATCTACAGAGTCATGGAACTGAAAGCGACAACGCTTGGAAAACGTCTGGCACAGCACCCTTACGATCGGGCGGTGATCCTCAATGC  >  NC_000913/1024276‑1024544

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