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,3310TC96.4% 79.4 / ‑5.8 28intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/13);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGC  >  NC_000913/1652092‑1652557
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    
gtttttttGTTTATACTGGAGTGCGCCCCGGTGAGATGTGGTCACGTGGATTGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGTATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                                                                                                                                                                                           <  2:521493/244‑1 (MQ=255)
       tGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTTATTAAGCCAGCAATCGAAGCTCTCCGGAGGCAGATGACATTAAAGAGACTGAGCAAAGAGAAAATCATTGTTGTTCACCTCCGAGAGGATGGcaga                                                                                                                                                                                                                     >  1:11584/1‑250 (MQ=255)
                 ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                                                                                                                                                                                            <  1:208548/249‑1 (MQ=255)
                       cGCCCTTGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGTAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                                                                                                                                                     <  1:518725/250‑1 (MQ=255)
                             ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                                                                                                                                                            <  1:424977/251‑1 (MQ=255)
                                                   gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAGAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAAGAAACAGGGGCATACATCTTATTTAGCCAGCAAACACCACCATCCGGCGTCAGGATGACATAAAGAGACTTAGCAAAGAACATATCATTTGTGATACACTCCACGAGTAAGGCAGAAAAGAAAAACAAAAAAGCCACATTTTTTTTCAACCttcatggta                                                                                                                                                                       >  2:567364/1‑242 (MQ=255)
                                                                       aaTGGAACAATGATGATCATGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCCATTCTCTTAATAAGCAAACAATCGGCGCCCTCCGGGGGAAGATGACATTAACGGGACTGGGGGAAACGGATATCAATGGTGTTCACCTCAGGGAGGATGGGAGAAAAGAAAAACAAAAATGGCCCTTTGTTTTTTAACCCGAAGTGTCAGAGAGGGGAAAAAAtttatg                                                                                                                                                  >  1:306307/1‑247 (MQ=255)
                                                                                        cAGGATAAATTTAGCAAAAGACCTGTTCACTGTACCAAAAACACAGGTGGGAACCATTCTGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGCATGGCAGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCTGAAGTTGCAGCGAGAGTACAAAATTATGGTGACCAGGTTACCGTTg                                                                                                                                 <  2:591149/251‑1 (MQ=255)
                                                                                        cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCAACGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAAAGCATATCATTGATGTTAAACTCCGAGAATAAGGCAGAAAAGCAAAACAAAAATGACACTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAAAATTATGGGGACCATTTTACCggtg                                                                                                                                 >  1:657799/1‑248 (MQ=255)
                                                                                                aTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACt                                                                                                                              <  2:417371/246‑1 (MQ=255)
                                                                                                aTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACt                                                                                                                              >  1:417371/1‑246 (MQ=255)
                                                                                                            aCCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                                                                                                                                                                        <  2:494911/176‑1 (MQ=255)
                                                                                                                   cACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCGCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCGCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga                                                                                                                                                                <  2:438759/193‑1 (MQ=255)
                                                                                                                   cACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAAAAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga                                                                                                                                                                >  1:438759/1‑193 (MQ=255)
                                                                                                                                                  cGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGAATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGCGACCTTTTTACCGTTGACTAAATAAGGCAGATGTGGGACTAAGCGATAAAACGTGCCGGACTCCGACATCGaa                                                                        >  2:15053/1‑250 (MQ=255)
                                                                                                                                                    ggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAAt                                                                     >  1:146623/1‑251 (MQ=255)
                                                                                                                                                     gTCATTCATCTTAATAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTTCCGGACTCCGCCATCGAAAATc                                                                    <  2:84725/251‑1 (MQ=255)
                                                                                                                                                         attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCTGAGAGTAt                                                                                                                                                                                                                           <  2:378190/99‑1 (MQ=255)
                                                                                                                                                         attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                                                                                                                                                                                           >  1:378190/1‑99 (MQ=255)
                                                                                                                                                                 tattaAGCCTGCCATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTTAGCAAAGAGCATTACATTGATGTTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCGCCTTTGTGTTTCACCCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGTCCATTTTACCGTTGACTCAATAAGGCAGATGTGGTACGCAGCGATAAAACGTGCCGGACTCCGCTATCTAAAATCATATCAGt                                                            <  2:601091/248‑1 (MQ=255)
                                                                                                                                                                 tattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAGTCCGCCATCGAAAATCATATAAGTCg                                                          >  2:460601/1‑250 (MQ=255)
                                                                                                                                                                             aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATCTCAGTAGAGAACTACTTATg                                             >  2:658967/1‑251 (MQ=255)
                                                                                                                                                                             aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                             >  1:264907/1‑251 (MQ=255)
                                                                                                                                                                                                  gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACa                                                     <  1:583614/222‑1 (MQ=255)
                                                                                                                                                                                                                     gAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAttt       <  1:658967/250‑1 (MQ=255)
                                                                                                                                                                                                                        cAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat     >  1:593737/1‑249 (MQ=255)
                                                                                                                                                                                                                         aaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAACTCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc  >  1:231691/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaGAGAATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATAAGGGACGAAGCGATAAAACGTGCCGGACTCCGCCCTCACAAATCCTATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGCCACCTGGTGCTAACCCGGCCGTTATAGc  >  2:221643/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    
GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGC  >  NC_000913/1652092‑1652557

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