Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,251,378	G→A	intergenic (+904/‑1128)	malM → / → ubiC	maltose regulon periplasmic protein/chorismate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,251,378	0	G	A	93.3%	36.8 / ‑3.3	15	intergenic (+904/‑1128)	malM/ubiC	maltose regulon periplasmic protein/chorismate lyase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (6/8);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GATAATGTCATATTTAGATGTCGACGTTTACACAATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGATCTCGATGGTTATATCAATAACATTATAAAAATGATTGATACATTGCCAGATAATGCAATGATATTGAAATCCGTTCTGGCCGTAAAACTGGTGATGCAATTAAAAATAC  >  NC_000913/4251156‑4251624                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       gATAATGTCATATTTAGATGTCGACGTTTACACAATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAAGTCAAACTTCGATAATGGTGATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAAAAAGATCATTTATTAAATCGAATTTGAc                                                                                                                                                                                                                            <  1:76503/251‑1 (MQ=255)       gTCATATTTAGATGTCGACGTTTACACAATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAACAGAGTTTCTGTCTGCCAAAAAAATATCAGTTATTAAATCGAATTGGACAGCTg                                                                                                                                                                                                                       <  1:4840/250‑1 (MQ=255)                                                                    tCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCACATTGGTCCTAGCAATATTCGCGACTGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTg                                                                                                                                                         >  1:209110/1‑251 (MQ=255)                                                                                                    gACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTCTCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAATATCATTTATTAAATCGAATTTGAc                                                                                                                                                                                                                            >  1:33395/1‑152 (MQ=255)                                                                                                        tctcGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGGTATGTTTGATTATGTTCGAATGTCGACAGGTAATTTTAAAGATTGCATTACAGAACGAGTTGAATTAATTATGGATTAGTCAGAGATATTTTGGa                                                                                                                       <  1:77887/249‑1 (MQ=255)                                                                                                            gCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAAAAATATCATTTATTAAANCGANNNNNAc                                                                                                                                                                                                                            >  2:129947/1‑144 (MQ=255)                                                                                                                               aGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAAtt                                                                                                                                                                                   >  1:8788/1‑166 (MQ=255)                                                                                                                                         gcgAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGAtc                                                                                                              >  1:49846/1‑225 (MQ=255)                                                                                                                                                 cAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCTCACTTCTCCATAGCAGATTTTCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTGTGTTCGAATGTCGACGGGGTATTTTAAAGATTTCATTAGATAACAATTTGAATTAACTTTTGATTATTCAGTTAGATTTTGGAATGAATAGCTCGATGGTTATATCAATACCATTATAAAAAtgat                                                                            <  1:89667/251‑1 (MQ=255)                                                                                                                                                     ttAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCACAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTAGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGATCTCGATGGTTATATCAATAACATTATAAAAAtgattga                                                                         >  1:117674/1‑250 (MQ=255)                                                                                                                                                                          aTNGTTNTNNTNNCGCNTCGCACTTCACCAGAGCAGAGNNTCTGTCTGCNAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTNCGAATGTCGNNAGNGAATGNTNNNGNNNGCNTTACAGAACAATTGGAATTNACNNTTGATTNTTCAGATATATTTTGGAATGAAGAtc                                                                                                              <  2:49846/192‑1 (MQ=255)                                                                                                                                                                                                  tNNCCAGAGNNGAGTTTCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAAtt                                                                                                                                                                                   <  2:8788/99‑1 (MQ=255)                                                                                                                                                                                                                 ttCTGTCTGCCAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGATCTCGATGGTTATATCAATAACATTATAAAAATGATTGATACATTGCCAGATAATGCAATGATATTGAAATCCGTTCTGGCCGTAAAACTGGTGATGCa             <  1:3260/250‑1 (MQ=255)                                                                                                                                                                                                                     gtctgCCAAAAAAATATCATTTATTAAATCGAATTTGCCAGCTGTTATGTTTGATTATGTGCGACTGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGATCTCGATGGTTATATCAATAACATTATAAAAATGATTGATACATTGCCAGATAATGCAATGATATTGAAATCCGTTCTGGCCGTAAAACTGGTGATGCAATTa         >  1:67426/1‑250 (MQ=255)                                                                                                                                                                                                                           cAAAAAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGATCGCGATGGTTATATCAATAACATTATAAAAATGATTGATACATTGCCAGATAATGCAATGATATTGAAATCCGTTCTGGCCTTAAAACTGGTGATGCAATTAAAAATan  <  1:129947/251‑2 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       GATAATGTCATATTTAGATGTCGACGTTTACACAATGTGATTTTTATCAAAGCGAGTGGTGAATGTGTCGATTTCAGCAAAAATATTCTTGATACAGTTGACTTCTCGCAGAGTCAACTTGGTCATAGTAATTTTCGCGAATGTCAGATTAGAAATTCAAACTTCGATAATTGTTATCTTTACGCTTCGCACTTCACCAGAGCAGAGTTTCTGTCTGCCAAAGAAATATCATTTATTAAATCGAATTTGACAGCTGTTATGTTTGATTATGTGCGAATGTCGACAGGGAATTTTAAAGATTGCATTACAGAACAATTGGAATTAACTATTGATTATTCAGATATATTTTGGAATGAAGATCTCGATGGTTATATCAATAACATTATAAAAATGATTGATACATTGCCAGATAATGCAATGATATTGAAATCCGTTCTGGCCGTAAAACTGGTGATGCAATTAAAAATAC  >  NC_000913/4251156‑4251624

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