Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,652,331 T→C intergenic (+2346/+397) ydfD → / ← ynfP lysis protein/protein YnfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3310TC92.0% 65.0 / ‑5.2 25intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/6);  minor base G (1/1);  total (18/7)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652087‑1652575
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGCAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGAC‑AATCCTCTTTTTAAGCCCGCACTC‑CACGCTCTCCGGCGGCAGATGACAATAACGAGACTGAGC‑AAAAAGCATATCATTGATGTTCACCTCag                                                                                                                                                                                                                                                  >  1:59244/1‑249 (MQ=255)
             gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCCATC‑GACGCTCTCCGGGGGCCGATGACATTAACGAGACTGAGC‑AAAAAGCAAATCATTGATGTTCACCCCCGAGAGTATGGgagaa                                                                                                                                                                                                                                    >  1:117469/1‑250 (MQ=255)
             gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTTAGATGACATTAACGAGACGGATA‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcaga                                                                                                                                                                                                                                     <  1:50614/249‑1 (MQ=255)
               ttATACTGGAGTTCGCCCTGGTGAGCTGCGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAGAC‑AAACAGGTA‑ATTCATCTTATTAAGCCAGCAAAA‑AACGCCACCCAGGGTCCGATGACAATAACGAGACTGCGC‑CAACAGCATATCATTGACGTTCACCTCATAGAGTAACGcagaaca                                                                                                                                                                                                                                  >  2:243108/1‑250 (MQ=255)
                  taCTGGAGTTCGGCCTGGTGAGCTGTGGGCACTTGGATGGGAGGACATAGATCTGCAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAAAC‑AAACGGGTCCATT‑TTATTATTAATCCCGCACAC‑GACCCCCTCCGGGGTCAGATGCAATTAACGCGACAGAGA‑CAAGAGCATATTCTTGATGTTCACCTCAGAGAGGATGGcagaacaaaa                                                                                                                                                                                                                               >  2:200596/1‑247 (MQ=255)
                   aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCATAACAGaaa                                                                                                                                                                                                                              <  2:45629/250‑1 (MQ=255)
                         gTTCGCCCTGGTGAGCTGGGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGCTGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCAAT‑GGGGC‑CTTCATCTTATTAAACCAACAAACTGA‑GCTTTCCGGAGGAAGAAGACATTTAAGAGACTG‑GCGAAAGAGAAAATCCTTTATGTTCACCTCCGCAAGTATGGCAGAACATGAAAACaaa                                                                                                                                                                                                                        >  2:85919/1‑250 (MQ=255)
                         gTTCGCCATGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACC‑AATCGGGGC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGGCATTAACGAAAACGAGC‑AAAAAGCATATCATTGATGTTCACCTCCGAGAGTATGGCAGAACAGAaaaaaaaa                                                                                                                                                                                                                        >  1:275036/1‑250 (MQ=255)
                            cGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                                                                                                                                                                     <  2:69791/250‑1 (MQ=255)
                              cccTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAAGC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                                                                                                                                                                                                  <  2:117469/251‑1 (MQ=255)
                                            gtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaa                                                                                                                                                                                                    <  1:2648/251‑1 (MQ=255)
                                                  ttGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGG‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAAAGCACCCTTTTTTTTTAACCTGa                                                                                                                                                                                               >  1:866/1‑250 (MQ=255)
                                                                       tGAAAAATGGCACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGGC‑ATTCATCTTATTAAGCCAGCAAAC‑GACGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAAGGTCAGCGAGAGTAAAAAAt                                                                                                                                                                         >  1:97633/1‑251 (MQ=255)
                                                                                             cAGGAGAAATTTAGCAAAAGCCCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAACCAACAATC‑GACGCTCTCCCGCGTCAGATGACATTAACGAGACTGAGC‑AAAAAGCATCTCATTGATGTTCACCTCAGAGAGTATGGCCGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCCGAAATGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcccgtt                                                                                                                                                   >  2:299902/1‑247 (MQ=255)
                                                                                                               aGACCGTTTCACGGTACCAAAAACCCAGGCGGGAACC‑AATCTGGTC‑AGTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATGAACGAGACTGAGC‑ATAGAGCATATCATTGATGTTCACGTCAGAGAGTATGTGAGAACAGAAAAACAAAAATGCGCCTTTTTTTTTCATCCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTTAACAGTTTACCTTAGTCTCAATAAGGCAGAtgt                                                                                                                                 <  2:346921/251‑1 (MQ=255)
                                                                                                                         aCGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGc                                                                                                                       <  2:866/251‑1 (MQ=255)
                                                                                                                                                         gggTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaaa                                                                                         >  2:107856/1‑250 (MQ=255)
                                                                                                                                                                   atCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                                                 >  1:19809/1‑249 (MQ=255)
                                                                                                                                                                                           gCTCTCCGGAGTCAGATGACCTTCACGAGACTGAGC‑AAAGAGCATATCCTTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACATTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAACAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGCCATACTTATGCctgctg                                                        >  2:296785/1‑251 (MQ=255)
                                                                                                                                                                                                ccGGAGTCAGATGACATTCACGAGACTGAGC‑AAAGAGCATATCCTTGATGTTCACCTCAGAGAGTATGGAAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGc                                                   >  1:60767/1‑251 (MQ=255)
                                                                                                                                                                                                  ggAGACAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGAGAGAGAAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGCTCAACCGTCCCGGACTACACCATCGACAATCATATCAGTCGAGACATACTTATGCATGCCGGTCGCTg                                                 >  2:254166/1‑251 (MQ=255)
                                                                                                                                                                                                                 tAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                                  >  1:13982/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaCGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAcc                                 >  1:188147/1‑251 (MQ=255)
                                                                                                                                                                                                                                    agCATATCATTGATGTTCACCTCCGAGAGTCTGGCAGAAAAGAAAAACAAAAATGCACCTTTGTTTTTCACCCTGCAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATACGGCCGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAc                >  2:295159/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ttCACGTCCGAGAGTATGGCAGACCAGCAAAACAAACATGCACCTTTGTTTTTCAACCTGAAGAGTCAGCGAGAGTAAAAAATCATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAAAGTGCAGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGc  >  2:23003/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTC‑ATTCATCTTATTAAGCCAGCAATC‑GACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGC‑AAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652087‑1652575

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