Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,768,441 (G)6→7 coding (245/3057 nt) ydiJ ← putative FAD‑linked oxidoreductase YdiJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,768,4411.G87.0% 57.9 / 4.5 24coding (245/3057 nt)ydiJputative FAD‑linked oxidoreductase YdiJ
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (6/14);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCGGCGCAAAAAAGTAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGG‑‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCACCAGTTGAACGACGCCGGGTGCCTGGGAAATCTGTGG  >  NC_000913/1768193‑1768679
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
gACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCGGCGCAAAAAAGTAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCggg‑ggggt                                                                                                                                                                                                                                                  <  2:162371/251‑3 (MQ=255)
        gCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCGGCGCAAAAAAGTAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCa                                                                                                                                                                                                                                       <  1:782/251‑1 (MQ=255)
                              cGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCGGCGCAAAAAAGTAGCAGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTgcgc                                                                                                                                                                                                                 >  1:181257/1‑251 (MQ=255)
                                                   tCGAAAGTTCCGGCGCAAAAAAGTAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCTCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGGGT‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACg                                                                                                                                                                                                        <  1:101348/240‑1 (MQ=255)
                                                   tCGAAAGTTCCGGCGCAAAAAAGAAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGGTGAGCGCCTGGCCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGG‑‑GGGAAAGATCAGCGATGAATAGCGTTCCTGCGCgggcagaag                                                                                                                                                                                                        >  2:101348/1‑233 (MQ=255)
                                                                  cAAAAAAGTAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATTTGGCGGGACATATGAACAATAATCCCCTGGTTGATCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCATCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGa                                                                                                                                                                             <  1:247410/251‑1 (MQ=255)
                                                                             ccGAAAGGTTTCAGGTACTGATTGAGTTTATCTGTTATCCCCCCGGCCTCAACGCGCACTCAACCCTCTTCAGGGTTAATTTCGATGATGCGTTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTTACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACAGCTGCGGGTGAACGCTGAAATa                                                                                                                                                                  <  2:250619/251‑1 (MQ=255)
                                                                                  cGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTTGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGGTGGG‑‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGGGATCAGCGCCACATCTGCGGTTGAACGCGGAAATaccacc                                                                                                                                                             <  1:69908/251‑1 (MQ=255)
                                                                                      ttCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCAt                                                                                                                                                         >  1:124879/1‑251 (MQ=255)
                                                                                        cAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGG‑‑GG‑AAAGATCAGCGATGAATAGCGGTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCg                                                                                                                                                       >  1:109295/1‑251 (MQ=255)
                                                                                             aCTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGga                                                                                                                                                   <  2:124879/250‑1 (MQ=255)
                                                                                                        ggatgttttATCACCCCGGCCTCACCGCGCACCCAGCCCTCTTCAGTGAGAATTTCGATGATGCGGTTCATATGGCGGGAGATGTCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGATCGCGGAAATACCACCGCATCGGGGAGGATTTGGTaa                                                                                                                                        <  2:312888/245‑1 (MQ=255)
                                                                                                                                                                   gatgCGGTTCATATGGCGGGACATAGCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACTCGGAAATACCACCGCATCGGGGATAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCTTATCGCCGGTa                                                                            <  2:247072/251‑1 (MQ=255)
                                                                                                                                                                           tCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAA‑TACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCttg                                                                    <  1:119669/250‑1 (MQ=255)
                                                                                                                                                                            cATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCttgt                                                                   <  1:181643/251‑1 (MQ=255)
                                                                                                                                                                                         atatCAACAATACTCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCtt                                                      >  2:107254/1‑251 (MQ=255)
                                                                                                                                                                                         atatCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCtt                                                      <  2:109295/251‑1 (MQ=255)
                                                                                                                                                                                               acaaTAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCaaaa                                                <  1:366171/251‑1 (MQ=255)
                                                                                                                                                                                                         ccTGGTTGAGCGCCTGACCGTTAGTGCCGGCGGCGCAGCACCGGG‑GGGA‑T‑‑AAAGATCAGAGATGAAAAGCGTTCCTGCGCGGCAAGCCGGCCGATCAGCGCCACATATGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGCAGTTGGAACATCCTGTTGTCGGTCGACATTGTCAGACGCTCGGCATCACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCAc                                      >  1:223778/1‑251 (MQ=255)
                                                                                                                                                                                                          cTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCAcc                                     >  1:12330/1‑251 (MQ=255)
                                                                                                                                                                                                                gAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCACCAGTTGa                               <  1:239075/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccggtgccgCCGCCGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAAGTAAGCACCAGTTGAACGACGCCGGGTGCCTggg            <  1:27767/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gccgccgccGCGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGATCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCACCAGTTGAACGACGCCGGGTGCCTGGGAAATc       >  2:30717/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ccgccgcGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCACCAGTTGAACGACGCCGGGTGCCTGGGAAATCtgtg   <  1:32941/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cgccgcGGG‑GGGG‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCACCAGTTGAACGACGCCGGGTGCCTGGGAAATCTGTgg  >  1:93608/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
GACCGGATGCATCGGTATTGATCATCCCACCGAGCGTTGCCCGGTTGCTGGTCGAAAGTTCCGGCGCAAAAAAGTAGCCGAACGGTTTCAGGTACTGATTGAGTTGATCTTTTATCACCCCGGCCTCAACGCGCACCCAGCCCTCTTCAGGGTTAATTTCGATGATGCGGTTCATATGGCGGGACATATCAACAATAATCCCCTGGTTGAGCGCCTGACCGTTAGTGCCGGTGCCGCCGCCGCGGG‑GGG‑‑T‑‑AAAGATCAGCGATGAATAGCGTTCCTGCGCGGCAAGACGGGCGATCAGCGCCACATCTGCGGTTGAACGCGGAAATACCACCGCATCGGGGAGAAGTTGGTAAATACTGTTGTCGGTCGACATTGTCAGACGATCGGCATAACTTGTCGCCGTATCGCCGGTAAAACCTTGTTGCTCCAGCTCTTGCAAAAAATTAAGCACCAGTTGAACGACGCCGGGTGCCTGGGAAATCTGTGG  >  NC_000913/1768193‑1768679

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑