Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,193,873 G→T 70.3% 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,8730GT70.3% 11.8 / 13.0 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 (3/2);  new base T (6/6);  total (9/8)
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.52e-01

ATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTC  >  CP009273/1193742‑1194016
                                                                                                                                   |                                                                                                                                               
aTTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAg                                                                                                                               >  2:1209856/1‑150 (MQ=255)
  tatataCTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAgaa                                                                                                                             >  1:3266219/1‑150 (MQ=255)
     ataCTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAgaaatg                                                                                                                          >  2:1075531/1‑150 (MQ=255)
               gTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaa                                                                                                                <  1:1931537/150‑1 (MQ=255)
                  ggTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaata                                                                                                             >  1:474700/1‑150 (MQ=255)
                                                   gACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatgatgaaatga                                                                            <  2:1626586/150‑1 (MQ=255)
                                                    aCCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatgatgaaatgat                                                                           >  1:3782693/1‑150 (MQ=255)
                                                    aCCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatgatgaaatgat                                                                           >  1:3782697/1‑150 (MQ=255)
                                                                                gATCAGTAGCTCGTATCGTGGTGGAGTAGCGCTACAAATTTGTCGCTGGGTTAGAATCTGTATTTGATAGATATTAAATAATGAAATAATGATATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTcc                                               <  1:3361210/150‑1 (MQ=38)
                                                                                gATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTcc                                               <  1:2837900/150‑1 (MQ=255)
                                                                                gATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTcc                                               <  1:2837964/150‑1 (MQ=255)
                                                                                gATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTcc                                               <  1:3361188/150‑1 (MQ=255)
                                                                                  tCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGATTTTCCAg                                             >  1:4607361/1‑150 (MQ=255)
                                                                                  tCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAg                                             >  1:4607429/1‑150 (MQ=255)
                                                                                               tCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATgg                                <  2:4611015/150‑1 (MQ=255)
                                                                                               tCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATgg                                <  2:4610982/150‑1 (MQ=255)
                                                                                                                             cTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTc  >  1:2272690/1‑150 (MQ=255)
                                                                                                                                   |                                                                                                                                               
ATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTC  >  CP009273/1193742‑1194016

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