Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	467,244	G→A	T23T (ACC→ACT)	ybaE ←	uncharacterized protein YbaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	467,244	0	G	A	94.1%	36.7 / NA	18	T23T (ACC→ACT)	ybaE	uncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (11/5);  minor base C (1/0);  total (12/6)
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.93e-01

GCTCACCTGGGGCCAGTTGCGCCAGCTCCAGCACATCTTGCTGCTTTCCGGTTTCCAGTGCCTGTTCCATCATCGCATTGCGTAGCGATTCCGGCGTGACCAGAAAGCGTAATTGTCCGCGCTTTCCGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACGGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGC  >  NC_000913/467003‑467444                                                                                                                                                                                                                                                  |                                                                                                                                                                                                         gCTCACCTGGGGCCAGTTGCGCCAGCTCCAGCACATCTTGCTGCTTTCCGGTTTCCAGTGCCTGTTCCATCATCGCATTGCGTAGCGATTCCGGCGTGACCAGAAAGCGTAATTGTCCGCGCTTTCCGCGTCCGGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTGGCggn                                                                                                                                                                                                 <  1:125592/251‑2 (MQ=255)              cAGTTGCGCCAGCTCCCGCAAATCTTGCTGCTTTCCGGTTTCCAGTGCCTGTTCCATCATCGCCTTGCGTAGCGATTCCGGCGTGACCAGAAAGCGTAATTGTCCGCGCTTTCCGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCTTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGa                                                                                                                                                                                      >  1:118132/1‑249 (MQ=255)                                                                     atcatcGCATTGCGTAGCGATTCAGGCGTGACCAGAAAACGTAATTGTCCGCGCTTTCAGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGATTGCGGCTTTCCGGCGGACGGTTGCCACAGCCGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCg                                                                                                                             >  1:207733/1‑251 (MQ=255)                                                                                       gATTCCGGCGTGACCAGAAAGCGTAATTGTCCGCGCTTTCCGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGGTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACt                                                                                                           >  1:24598/1‑251 (MQ=255)                                                                                                  naCCAGAAAGCGTAATTGTCCGCGCTTTCCGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCcgg                                                                                                                                                                                          >  1:89949/2‑161 (MQ=255)                                                                                                                              ccGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACCGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGAGATACTCAAt                                                                    >  1:474784/1‑251 (MQ=255)                                                                                                                                                 gCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGCTCGCTGCAAAACCAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAAACGCATATAAACCCGGAACAATATTCTTTAACTATTCACTATTACtt                                                                                             >  1:322722/1‑207 (MQ=255)                                                                                                                                                            cATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTgg                                           <  1:218308/246‑1 (MQ=255)                                                                                                                                                                      cctGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGc                              <  1:4604/249‑1 (MQ=255)                                                                                                                                                                                          gTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTAttt                                                                                                               >  1:29624/1‑148 (MQ=255)                                                                                                                                                                                             cGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGACCAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTc                         >  1:277361/1‑231 (MQ=255)                                                                                                                                                                                                aCATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGc  <  1:138187/251‑1 (MQ=255)                                                                                                                                                                                                     gCGTTCGCTGCAAAAACAGCGTTCGGCCAGTGCGCTGACGGTGATGGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACgg                                                                                                                                                    <  1:290527/100‑1 (MQ=255)                                                                                                                                                                                                              gcgaaaaCAGCGCTCGGCCAGTTCTCTGACGGTGACAGTTTGCGGCGTGCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTTAGACGGTGGAGCAATGGCATATAAAGCCGCAACAATATTATTTAACTATTCACTATTACTTCCGtat                                                                                       <  1:214818/149‑1 (MQ=255)                                                                                                                                                                                                              gCAAANCCAGCNTNNCGCCAGTTCGCTGACGGTGNCAGTTTGCGGCTTTCNGGCGNAAGGTNGCna                                                                                                                                                                             >  2:214818/1‑66 (MQ=255)                                                                                                                                                                                                                      aGCGTTCAGCCAGTTCGCTGACGGTGACAGGTTGCGGCTTTCCGACGGAAGGTTGACACAGCCGTTGATACTGGTTCGGCCGGTTGAGCAATCACATATATACCAGGAACACTATTATTTAACTATTCACTATTACTTCAGAATATATTAGGTGATaatc                                                                       >  1:485878/1‑157 (MQ=255)                                                                                                                                                                                                                       ncGTTCGGCCAGTTCGCTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTc                            >  1:133966/2‑202 (MQ=255)                                                                                                                                                                                                                                       cTGACGGTGACAGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGgc                     >  1:8865/1‑193 (MQ=255)                                                                                                                                                                                                                                             gtgacagtTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATaaa                                                                                                                                 >  1:321737/1‑79 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                         GCTCACCTGGGGCCAGTTGCGCCAGCTCCAGCACATCTTGCTGCTTTCCGGTTTCCAGTGCCTGTTCCATCATCGCATTGCGTAGCGATTCCGGCGTGACCAGAAAGCGTAATTGTCCGCGCTTTCCGCGTCCTGACTGCGCCTGCCACTCCAGCCATCCCGCCTCCTGTGCCTGACGCAACAGCGTACGAACATGGCGTTCGCTGCAAAAACAGCGTTCGGCCAGTTCGCTGACGGTGACGGTTTGCGGCTTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGC  >  NC_000913/467003‑467444

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