Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	310	(T)9→8	intergenic (+55/‑27)	thrL → / → thrA	thr operon leader peptide/fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	310	0	T	.	81.8%	52.7 / 11.1	22	intergenic (+55/‑27)	thrL/thrA	thr operon leader peptide/fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (10/8);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

GAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGT  >  NC_000913/75‑543                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           gAACTGGTTCCCTGCCGTGAGTAAATTAAAATTTTATTGACTGATGTAACTAAAATCTTTAACCATTATAGGCATAGCTCACAGACAGATAAAATTTACAGAGTACACAACATCCATGATACTCATTAGCACCACCATTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTCCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACg                                                                                                                                                                                                                             <  2:384032/251‑1 (MQ=255)       gTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTaac                                                                                                                                                                                                                      <  1:119893/251‑1 (MQ=255)               ccGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                                                                                                                              <  1:97297/251‑1 (MQ=255)                                                                      aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCgggt                                                                                                                                                      >  1:81191/1‑248 (MQ=255)                                                                              cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGat                                                                                                                                               >  1:270520/1‑250 (MQ=255)                                                                                   cgacagaTACAAATTACAGAGTACCCAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAACAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGCAACGAAGTAACAACCATGCGAGTGTTGAAGTACGGAGGTACATCAGTGGCA‑AATGAAGAACGGTTTCTGGGGGTTGGCGATTTTCtt                                                                                                                                         >  2:245490/2‑250 (MQ=255)                                                                                             aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGAGCTTTTTTTT‑CGACCAAAGTTAACGGGGTAACAACCATGCGAGTTTTGAAGTTCTGCGGTCCATCAGTGGCT‑AATGCAAAACGTTTTCTGCGGGTTGCCGATCTTCTGGAAAGCAATg                                                                                                                               >  1:137131/1‑251 (MQ=255)                                                                                                   cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCcaggc                                                                                                                         <  1:245490/251‑1 (MQ=255)                                                                                                       gTACACAACCTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCAAggcaggg                                                                                                                     >  1:264046/1‑251 (MQ=255)                                                                                                              acaTCCAATAAACGCATTAGCCCCCCCATTCACACCACCATCACCATTACCACAGGTAACGGTGAGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                                                              <  2:81191/251‑1 (MQ=255)                                                                                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccct                                                                                        >  1:167112/1‑250 (MQ=255)                                                                                                                                          ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGCGGCCACCGTCCTCTCTGCCCCCGCCaa                                                                                    >  1:95047/1‑250 (MQ=255)                                                                                                                                                caaccatctccatTACCAAGGGAAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGACCGCCCATGACAGGGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACAATGCGAGTGTTTAAGTTCGGCGGTACATCAGTGGCA‑AAAGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGTAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAAGCacca                                                                            <  2:108720/249‑1 (MQ=255)                                                                                                                                                 caacatcaccatTACCACAGGAAACGGTGCGAGCTAACGCGGAAAGGAAACACAGAAAAAAGCCCGCACCTGAAAGTGCGGGCTTTTAATT‑AGACCAAAGGTAAAACAGAAACAACGATGAGAGAGTTGAAGTTCGGAGGTACATCAATGGCACAAGACAGA‑CGGTTTATGCGGGGATCCGATATTCTGACAATCACAGGCAGGCAGGGGCAGGAGGCCACCGTACATTCTGGCACCGGCAAAATCaccg                                                                            >  2:375699/4‑249 (MQ=255)                                                                                                                                                          catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGCGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGGCCTCTCTGCCCCCGCCCAAATCACCAACCACCtggt                                                                  >  2:69040/1‑251 (MQ=255)                                                                                                                                                                 acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                                                           <  1:281161/251‑1 (MQ=255)                                                                                                                                                                          ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AAAGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCTCCAAAATCACCAACCAACTGGTTGCGATGTATGaaaa                                                   >  2:270049/1‑250 (MQ=255)                                                                                                                                                                                                         aaaaaGCCCCCACCTTACAGTGCGTG‑TTTGTTTTTCGACAAATGTTAATGGGGTAACAACAAGGCGAGGGTTGACGTTCGGCGTTACATCAGTGGCA‑AATGCAGAACGTTTTCGGGGGGTTGCCGATATTCTGAAAAGCAATGCCAGGCAGGGGGAGGTGGCCACCGTCCTCTCTTCCCCCGCCAAAATCACCAACCCCCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAata                   <  2:326251/251‑1 (MQ=255)                                                                                                                                                                                                                      cTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCg       >  1:122160/1‑250 (MQ=255)                                                                                                                                                                                                                        gACAGTGCGGGCTTTTTTTT‑CGACCAAAGGGAATGAGGTAACAACCATGCGAGTGTGGAAGTTCGGCGGTACACCAGTGGCA‑AATGCAGAACGTTGTCTGCGGGTTGCCGATATTCTGGAAAGCAAGGCCAGGCAGGGTCAGTTGGCACCCGTCCTCACGGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGa      <  2:336001/249‑1 (MQ=255)                                                                                                                                                                                                                           aGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCttt                           <  1:78381/225‑1 (MQ=255)                                                                                                                                                                                                                           aGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTagttcggaagggcgtcgtgtagggg  >  2:78381/1‑226 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           GAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCA‑AATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGT  >  NC_000913/75‑543

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