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

ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑TTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGAT  >  NC_000913/1652087‑1652566
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCagag                                                                                                                                                                                                                                        <  2:47863/251‑1 (MQ=255)
     gtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTAt                                                                                                                                                                                                                                   <  1:374911/251‑1 (MQ=255)
            tGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACACTCTCCGGAGTCAGATGAAATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGc‑aga                                                                                                                                                                                                                            >  1:289049/1‑250 (MQ=255)
                           tCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCCATCAGGTCATTCATATTATTAAACCAGCAATCGACGCCCTCCGGAGGAAGATGCAATTAACGAGCCTGAGCCAAGAGGATAACATTGATGTTCACCTCCAG‑GAGTATGGCCAGAA‑AGAAAAACaaaaa                                                                                                                                                                                                             >  2:193642/1‑250 (MQ=255)
                                            gtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCGCCTTTGTTTTTCa                                                                                                                                                                                            <  1:117663/250‑1 (MQ=255)
                                                         ggaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCTTCTTATTAAACCAGCCATCGCCGCCCTCCGGGGACAGATGAAATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CCCAGAGAGTATGGC‑CGAACAGAAAAACAAAAAATCACCTTTGTTTTTCAACAAGAAGTGTc                                                                                                                                                                                >  2:303693/1‑249 (MQ=255)
                                                                          aaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AAACCATCGGGTCATTCCTCTTATTAAGCAAGCAAAAGAAGCTCTCCGGAGTCAGATGACATTAACGAGACTGGGCAAAAAGGATTTCATTGGTGTTCAC‑CTCAGAGAGAATGGA‑GGAACAAAAAAAAAAAAATTCACATTTGTTTTTTAACCTGAAGGGTCAGGGGGAGTAAAAAATTAt                                                                                                                                                             >  2:227041/1‑251 (MQ=255)
                                                                                                              aaGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGAtg                                                                                                                         <  1:303693/251‑1 (MQ=255)
                                                                                                                  ccGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTggg                                                                                                                     <  1:203389/251‑1 (MQ=255)
                                                                                                                    gTTTCACGGTACCAAAAACACAGGCGGGGAACCA‑TCTGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATTATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTTTTTTTCAACCTGAAGTGTCATCGAGGGTAAAAAAATATTGTGACCATTTTAACGTTGTCTCAATAATGCAGATGTgggg                                                                                                                    >  2:106261/1‑249 (MQ=255)
                                                                                                                                cAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaaa                                                                                                       <  1:106261/251‑1 (MQ=255)
                                                                                                                                                         gggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaaaa                                                                               <  1:129729/251‑1 (MQ=255)
                                                                                                                                                                tcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTGACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCatat                                                                         <  1:227041/250‑1 (MQ=255)
                                                                                                                                                                 catCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                                        <  1:372919/250‑1 (MQ=255)
                                                                                                                                                                     ttattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTGTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTc                                                                    <  2:45697/250‑1 (MQ=255)
                                                                                                                                                                                  aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                                      >  1:273732/1‑251 (MQ=255)
                                                                                                                                                                                               cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCt                                         <  1:234164/251‑1 (MQ=255)
                                                                                                                                                                                                                aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGAC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAAACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAcc                        >  1:370626/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTGAGCAAAGAGCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                 >  1:304241/1‑251 (MQ=255)
                                                                                                                                                                                                                                gagCATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCaa         >  1:127620/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cATATCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCa     <  1:193642/251‑1 (MQ=255)
                                                                                                                                                                                                                                      atCATTGATGTTCAC‑CTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGAt  <  2:21325/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCAC‑TTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGAT  >  NC_000913/1652087‑1652566

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