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% 101.1 / ‑5.0 28intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (13/14);  total (14/14)
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.97e-01

GGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACT‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652089‑1652575
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGAGTACTTTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCagaga                                                                                                                                                                                                                                                <  1:113883/250‑1 (MQ=255)
                 aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGaaa                                                                                                                                                                                                                               >  1:198555/1‑250 (MQ=255)
                                   gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAGTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAGAAGGCACCTTTGt                                                                                                                                                                                                            <  2:126543/251‑1 (MQ=255)
                                    aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGATCAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAAC‑AATCGGGTCATTCATCTTATTAAACCACCAATCGCCC‑CGCTCCGGAGTCAGAAGACATTTACGAGACTGAGCCAAGAGGATATCATTGATGTTC‑CTCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAAATAACCTTTGtt                                                                                                                                                                                                           >  2:132083/1‑251 (MQ=255)
                                             cACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                                                                                                  <  1:60524/251‑1 (MQ=255)
                                                    aTGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AAACGGGTCATTTATCTTATTAAACCAGCAATCGACG‑CTCTCCGGAGTCAGAAGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtgt                                                                                                                                                                                           >  1:812/1‑251 (MQ=255)
                                                      gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                                                                                                                                                                                          <  2:154696/250‑1 (MQ=255)
                                                             cATAGAGCTGAAAAATGGCACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCAGAAACCCAAACGG‑ACATTTATATTATTAAGACAGCAATAGCCGACTTTCCG‑AGACCGATGACATTAACGATACTTAGCAAAAAGCATATCATTTATGTTTACC‑TCAGCGAGTATGGACGAACAGAAAAACAATAATACACCTTTGTTTTTTAACCTTAAgtgacaagcag                                                                                                                                                                                    >  2:197473/1‑241 (MQ=255)
                                                                aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCAGC‑TCAGAGAGTATGGCATTACATAAATACAGAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑Ta                                                                                                                                                                               <  1:26308/250‑1 (MQ=255)
                                                                                  aaTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCAtt                                                                                                                                                            <  1:41404/251‑1 (MQ=255)
                                                                                   atgatgATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCAttt                                                                                                                                                           <  1:130941/251‑1 (MQ=255)
                                                                                         atCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AAAAGGGGCATTCCTCTTTTTAAACCCGCAATCCACG‑CCCTCCCGAGTCAGATGACATTTACGAGACTTAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAATATGGCAGAACAGAAAAACAAAAATGCCCCTTTTTTTTTTTACCTGAAGTGTCAGCGAGAGGAAAAAAATTAGGGT‑ACCATTTTTCCGt                                                                                                                                                     >  2:147429/1‑251 (MQ=255)
                                                                                             ggAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGAc                                                                                                                                                 <  2:160638/251‑1 (MQ=255)
                                                                                                gaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTc                                                                                                                                               >  1:51272/1‑250 (MQ=255)
                                                                                                  aaTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AAACGGGTCATTCCTCTTATTAAGCCAGCAATCGACG‑CCCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATCTCCTTGATGTTCACC‑TCCGAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGCCCATTTTACCGTTGACTCaa                                                                                                                                             >  1:137095/1‑250 (MQ=255)
                                                                                                                                               aCC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCa                                                                                               <  2:198555/251‑1 (MQ=255)
                                                                                                                                                    aTCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAACATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCAATCg                                                                                            >  2:140255/1‑250 (MQ=255)
                                                                                                                                                            attcatCTTATTAAGCCAGCAATCGCCG‑CTCTCCGGAGTCAGATGAAATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                                                                                                                                                                             >  1:251311/1‑130 (MQ=255)
                                                                                                                                                            attcatCTTATTAAGCCAGCAAGCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGTCAGAACAGAAAAACAAAAATGCACCTTTg                                                                                                                                                                                                             <  2:251311/130‑1 (MQ=255)
                                                                                                                                                              tcatCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGGTGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAc                                                                                                      >  2:233645/1‑230 (MQ=255)
                                                                                                                                                              tcatCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGGTGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCAGTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAc                                                                                                      <  1:233645/230‑1 (MQ=255)
                                                                                                                                                                atCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                                                 >  2:76772/1‑249 (MQ=255)
                                                                                                                                                                             aGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACAAGACTGAGAACAGAGAATATCATTGATGTTCACC‑TCAGCGAGTCTGCCCGAACAGAAAAACAAAAATGCACCTTTGTTTTTCACCCTGAAGTGTAAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGAAGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACtt                                                                  >  1:167236/1‑251 (MQ=255)
                                                                                                                                                                                       g‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                                        <  1:70051/251‑1 (MQ=255)
                                                                                                                                                                                               ggAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTg                                                 <  1:213375/251‑1 (MQ=255)
                                                                                                                                                                                                                           cAAAGAGCATATCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGCG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                       >  1:112295/1‑249 (MQ=255)
                                                                                                                                                                                                                                     atCATTGATGTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc               >  1:60670/1‑247 (MQ=255)
                                                                                                                                                                                                                                              gTTCACC‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTGGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCTGCATTTATAGGAAACCAGATGTGGCATGc  <  2:112295/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
GGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACC‑AATCGGGTCATTCATCTTATTAAGCCAGCAATCGACG‑CTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACT‑TCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAG‑TAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652089‑1652575

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