Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,193,873 G→T 100% intergenic (‑180/‑278) ymfE ← / → lit e14 prophage; putative inner membrane protein/e14 prophage; cell death peptidase, inhibitor of T4 late gene expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,193,8730GT82.2% 25.3 / 5.1 17intergenic (‑180/‑278)ymfE/lite14 prophage; putative inner membrane protein/e14 prophage; cell death peptidase, inhibitor of T4 late gene expression
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (7/7);  total (8/9)
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.04e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTAATTAAGTGTCATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTA  >  CP009273/1193728‑1194009
                                                                                                                                                 |                                                                                                                                        
tttAATTAAGTGTCATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGaa                                                                                                                                      <  2:4201640/150‑1 (MQ=255)
                     aCTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAgaaatgaa                                                                                                                 <  2:4412489/150‑1 (MQ=255)
                       ttttCCGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGAACTAGCTACGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATAGATAgaaatgaaat                                                                                                               >  2:2196069/1‑150 (MQ=255)
                       ttttCAGGAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAgaaatgaaat                                                                                                               >  2:2695481/1‑150 (MQ=255)
                                  ttagttaCAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaataat                                                                                                    <  2:1213413/150‑1 (MQ=255)
                                    agttaCAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaataatga                                                                                                  >  1:983640/1‑150 (MQ=255)
                                    agttaCAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaataatga                                                                                                  >  1:983631/1‑150 (MQ=255)
                                              ttttGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaa                                                                                        <  1:1453698/150‑1 (MQ=255)
                                              ttttGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaa                                                                                        <  1:238760/150‑1 (MQ=255)
                                                  gTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgat                                                                                    >  1:1618270/1‑150 (MQ=255)
                                                        gTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatg                                                                              >  1:1409338/1‑150 (MQ=255)
                                                        gTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGAGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatg                                                                              >  1:1409402/1‑150 (MQ=255)
                                                                                        tatTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACaga                                              <  1:1227501/150‑1 (MQ=255)
                                                                                                                 ggaggagTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTc                     <  1:3006981/150‑1 (MQ=255)
                                                                                                                               aCAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTa       >  1:4330112/1‑150 (MQ=255)
                                                                                                                                    tttGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTa  <  1:3882376/150‑1 (MQ=255)
                                                                                                                                    tttGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTa  <  1:2845175/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
TTTAATTAAGTGTCATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTA  >  CP009273/1193728‑1194009

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