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

GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGC  >  NC_000913/1652092‑1652557
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    
gtttttttGTTTATACTGGAGTGCGCCCCGGTGAGATGTGGTCACGTGGATTGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑GTATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTAt                                                                                                                                                                                                                           <  2:521493/244‑1 (MQ=255)
       tGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGG‑GCATTCATCTTATTAAGCCAGCAATCGAAGCTCTCCGGAGGCAG‑ATGACATTAAAGAGACTGAGCAAAGAGAAAATCATTGTTG‑TTCACCTCCGAGAGGATGGcaga                                                                                                                                                                                                                     >  1:11584/1‑250 (MQ=255)
                 ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                                                                                                                                                                                            <  1:208548/249‑1 (MQ=255)
                       cGCCCTTGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGTAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                                                                                                                                                     <  1:518725/250‑1 (MQ=255)
                             ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑Ctt                                                                                                                                                                                            <  1:424977/251‑1 (MQ=255)
                                                   gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAGAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAAGAAACAGGG‑GCATACATCTTATTTAGCCAGCAAACACCACCATCCGGCGTCAGGATGACAT‑AAAGAGACTTAGCAAAGAACATATCATTTGTGATACACTCCA‑CGAGTAAGGCAGAAAAGAAAAACAAAAAAGCCAC‑ATTT‑TTTTTCAACCttcatggta                                                                                                                                                                       >  2:567364/1‑242 (MQ=255)
                                                                       aaTGGAACAATGATGATCATGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCCATTC‑TCTTAATAAGCAAACAATCGGCGCCCTCCGGGGGAAG‑ATGACATTAACGGGACTGGGGGAAACGGATATCAATGGTG‑TTCACCTCAGGGAGGATGGGAGAAAAGAAAAACAAAAATGG‑CC‑CTTTGTTTTTTAACCCGAAGTGTCAGAGAGGGGAAAAAAtttatg                                                                                                                                                  >  1:306307/1‑247 (MQ=255)
                                                                                        cAGGATAAATTTAGCAAAAGACCTGTTCACTGTACCAAAAACACAGGTGGGAACCATTCTGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGCATGGCAGAACAGAAAAACAAAAATGC‑AC‑CGTTGTTTTTCAACCTGAAGTTGCAGCGAGAGTACAAAATTATGGTGACCAGGTTACCGTTg                                                                                                                                 <  2:591149/251‑1 (MQ=255)
                                                                                        cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCAACGGG‑GCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAAAGAGACTGAGCAAAAAGCATATCATTGATG‑TTAAACTCCGAGAATAAGGCAGAAAAGCAAAACAAAAATG‑‑ACACTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAAAATTATGGGGACCATTTTACCggtg                                                                                                                                 >  1:657799/1‑248 (MQ=255)
                                                                                                aTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACt                                                                                                                              <  2:417371/246‑1 (MQ=255)
                                                                                                aTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑GCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACt                                                                                                                              >  1:417371/1‑246 (MQ=255)
                                                                                                            aCCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGtt                                                                                                                                                                                        <  2:494911/176‑1 (MQ=255)
                                                                                                                   cACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCGCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCGCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCgaga                                                                                                                                                                <  2:438759/193‑1 (MQ=255)
                                                                                                                   cACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑GCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGGATGGCAGAACAGAAAAAAAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCgaga                                                                                                                                                                >  1:438759/1‑193 (MQ=255)
                                                                                                                                                  cGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGAATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGCGACCTTTTTACCGTTGACTAAATAAGGCAGATGTGGGACTAAGCGATAAAACGTGCCGGACTCCGACATCGaa                                                                        >  2:15053/1‑250 (MQ=255)
                                                                                                                                                    gg‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAAt                                                                     >  1:146623/1‑251 (MQ=255)
                                                                                                                                                     g‑TCATTCATCTTAATAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTTCCGGACTCCGCCATCGAAAATc                                                                    <  2:84725/251‑1 (MQ=255)
                                                                                                                                                         attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCTGAGAGTAt                                                                                                                                                                                                                           <  2:378190/99‑1 (MQ=255)
                                                                                                                                                         attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTAt                                                                                                                                                                                                                           >  1:378190/1‑99 (MQ=255)
                                                                                                                                                                 tattaAGCCTGCCATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTTAGCAAAGAGCATTACATTGATG‑TTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑GC‑CTTTGTGTTTCACCCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGTCCATTTTACCGTTGACTCAATAAGGCAGATGTGGTACGCAGCGATAAAACGTGCCGGACTCCGCTATCTAAAATCATATCAGt                                                            <  2:601091/248‑1 (MQ=255)
                                                                                                                                                                 tattaAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAGTCCGCCATCGAAAATCATATAAGTCg                                                          >  2:460601/1‑250 (MQ=255)
                                                                                                                                                                             aaTCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATCTCAGTAGAGAACTACTTATg                                             >  2:658967/1‑251 (MQ=255)
                                                                                                                                                                             aaTCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                             >  1:264907/1‑251 (MQ=255)
                                                                                                                                                                                                  g‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACa                                                     <  1:583614/222‑1 (MQ=255)
                                                                                                                                                                                                                     gAGCAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAttt       <  1:658967/250‑1 (MQ=255)
                                                                                                                                                                                                                        cAAAGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat     >  1:593737/1‑249 (MQ=255)
                                                                                                                                                                                                                         aaaGAGCATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAACTCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc  >  1:231691/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaGAGAATATCATTGATG‑TTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATAAGGGACGAAGCGATAAAACGTGCCGGACTCCGCCCTCACAAATCCTATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGCCACCTGGTGCTAACCCGGCCGTTATAGc  >  2:221643/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    
GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGG‑TCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAG‑ATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATG‑TTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGC‑AC‑CTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGC  >  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: ‑