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

GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652092‑1652573
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                                                                                                                                                                                                           <  1:123296/251‑1 (MQ=255)
gtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CAGTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTag                                                                                                                                                                                                                                           <  2:354925/251‑2 (MQ=255)
  gtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                                                                                                                                                                                          <  1:372082/250‑1 (MQ=255)
                        gCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCCTTCATCTTATTAAGCCCGCAATCGACGCCCTCCGGAGTCAGATGG‑CATTAAAGAGACTGAG‑CAAATAGCATATCATTTATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAAACAAAATGc                                                                                                                                                                                                                 >  1:530429/1‑251 (MQ=255)
                             ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCt                                                                                                                                                                                                             <  1:495348/250‑1 (MQ=255)
                                  gCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGCACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAAAGGGGCCTTCTCCTTTTTTAACACTCAAGCCGCGCGCTACGGAGTCAGATGA‑CATTTAAGAGCCATAGACAAA‑AGCATATAATTTATGTTCACTTCAGAGAGGATTGCAGAA‑CAAcaaaaacaaaaagaaccatttttt                                                                                                                                                                                                        >  2:723/1‑239 (MQ=255)
                                         gCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAAACAAAAGGGTCATTCATTTTATTAAGCCAGCACTAGACGCTCTCCCGAGTCAAATGA‑CATTAAAGAGACTGAG‑CACAAAGCATATCATTTATGTTCACCTCAAAGAGGATGGCAGAAACA‑‑AAAAAAAAAAATGGACACTTGttttttaaac                                                                                                                                                                                                >  2:2477/1‑246 (MQ=255)
                                                  tGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAAATGGTCATTCATCTTATTAAGCCAGCAATCCACGCTCTCCGGAGTAAGATGAACTTTAA‑GAGACCGAG‑CAAAAAGCAGATCATTTATGTTTACCTCAAAGAGAATGGGCGAA‑AA‑GAAAAACAAAAATGCAACCTTGTTTTTCAACCTTAAGGGTc                                                                                                                                                                                       >  2:190719/1‑251 (MQ=255)
                                                                     aaaaTGGAACAATGAGGAGCAGGAGAAATTTAGCGAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACAAATCGGGTCATTCATCTTAGTAAGCCAGCAATCGACACTCTCCGGAGTCAGATGA‑CATGAACGAGACTGAG‑CAAAGAGCATATCACTGAGGTTCACCTCATAGAGTAGGGCAGAA‑CA‑GAAAAACAAAAAAGCACCTTTGTTGTGCAACCTGAAGTGGCAGCGAGAGTAAAAAATTAt                                                                                                                                                                    <  2:315132/251‑1 (MQ=255)
                                                                                  gatgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTa                                                                                                                                                       <  1:318181/251‑1 (MQ=255)
                                                                                     gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAGACAATCGGGTCATTCATCTTATTAAGCCAGCCATCGACACACTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATACTTGGTGTTCACCTCAGAGGATATGGCAGAA‑CA‑AGAAAACAAAAAAGCCACTTTGGTTTTGAACCTGAAATGTAAAAGAGAGTAAAAAATTATGGTGAACAtttttccc                                                                                                                                                    >  1:362064/1‑247 (MQ=255)
                                                                                             gaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGGCATTCATCTTTTTAAGCCAGCAATCGACGCTCTCCCGAGTCAGATGA‑CATTAACGAGACTGAG‑CCAAAAGCCTATCATTTATGTTTACCCCAGAGAGTATGGCAGAA‑CA‑AAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGGTGACCATTTTACCGTTTACTc                                                                                                                                             >  1:188439/1‑250 (MQ=255)
                                                                                                      cAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAAACGGGGCATTCATCTTATTAAGCCAGCAATCGGCGCTCTCCGGCGGCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATAATTGGTGTTCACCTCCGAGAGGATGGCAGAA‑CA‑AAAAAACAAAAATGCCCCTTTTTTTTTCACCCTCAAGGGTCAGCAGGAGGAAAAAAATTTTGTGACCATTTTTCCGTTTACTCaaaaaagca                                                                                                                                    >  2:143969/1‑243 (MQ=255)
                                                                                                                                            aCCAATCGGGTCAGTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCa                                                                                             <  1:140004/251‑1 (MQ=255)
                                                                                                                                             ccAATCGGGTCATGCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAATATCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTTACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCAt                                                                                            <  2:530429/251‑1 (MQ=255)
                                                                                                                                                       cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                                                                                                                                                                            >  1:365219/1‑161 (MQ=255)
                                                                                                                                                       cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                                                                                                                                                                            <  2:365219/161‑1 (MQ=255)
                                                                                                                                                          tcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCatat                                                                                <  1:3130/250‑1 (MQ=255)
                                                                                                                                                                 attaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCgag                                                                        <  1:143969/251‑1 (MQ=255)
                                                                                                                                                                       ccAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAc                                                                  <  1:320889/251‑1 (MQ=255)
                                                                                                                                                                           cAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTAt                                                              <  1:114197/251‑1 (MQ=255)
                                                                                                                                                                             aTCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGCAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGGGCCGGCCTCCGCCATCGAAAATCATATCAGTCgggcaatacttaggc                                                            >  2:493343/1‑237 (MQ=255)
                                                                                                                                                                                           gAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTg                                               >  1:255244/1‑250 (MQ=255)
                                                                                                                                                                                                     a‑CATTAACGATACTGAG‑CAAAGCGCATATCATTGATGTTCACCTCAGAGAGTATGGCCGAA‑CA‑GCAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGTCTCAATAAGGCAGATGTGAGACGCAGCGATAAATAGTGCCGGCCTCAGCCATCGAACATCCTCTCAGTCGAGCAATACTTATCACTGCTGGTcactaccagctgttgc                                    >  2:455789/1‑236 (MQ=255)
                                                                                                                                                                                                         ttAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCGGCTGGTCGCTGACAGCTGGTGCTa                                  >  1:289663/1‑250 (MQ=255)
                                                                                                                                                                                                           aaCGAGACTGAG‑CAAAGAGCATATCAATGAGGTTCACCTAAGAGAGGATGGCAAAA‑CA‑GAAAAACAAAAATGCACCTTTTTTTTTAAACCCGAAGTGTCAGCGAGATTAAAAAATTATGGTGACCATTTTACCGTTGACCCAATAAGGCAGATTTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTAGGCCTGCTGTTCGCTGACAGCTGGTGCTAAc                                <  2:218227/250‑1 (MQ=255)
                                                                                                                                                                                                               agaCTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGc                           >  1:175868/1‑251 (MQ=255)
                                                                                                                                                                                                                  cTGAG‑CAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                        >  1:463719/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ttGATGTTAACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTACGCCATCGAAAATCATATCAATCGAGACATAATTATGCCTGATGGACGCTGACAGCTGGTGCTAACACGGCATTTATAGCAAACACGATGGGc     >  2:121661/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aggTTCACCTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGGGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCGGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  <  2:475284/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGA‑CATTAACGAGACTGAG‑CAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAA‑CA‑GAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652092‑1652573

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