Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,430,584	C→T	P241S (CCA→TCA)	cycA →	D‑serine/alanine/glycine/:H(+)symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,430,584	0	C	T	88.5%	55.9 / ‑0.1	26	P241S (CCA→TCA)	cycA	D‑serine/alanine/glycine/:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (16/7);  total (18/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 = 8.52e-01

AATCGTCGCCATTGTGTCGCTGATTGTC‑GTCGGCCTGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTG  >  NC_000913/4430367‑4430829                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      aaTCGTCGCCATTGTGTCGCTGATTGTC‑GTCGGCCTGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCa                                                                                                                                                                                                                      <  1:138238/251‑1 (MQ=255)     gtcgCCATTGTGTCGCTGATTGTC‑GTCGGCCTGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGGTGTGGCAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGGTTAAGTGGCTTCTTTGCCGGATTCCAGATAGAGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGGAATCACTGCCACGCGCGATTAACTCCATTcc                                                                                                                                                                                                                  >  1:140976/1‑251 (MQ=255)       ccctctcgtgtCGCTGATTGTCTGT‑GTCCTGGACATGGTGGCGATGC‑CTGTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACTGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGAGAGCGGTTTGCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTTCGGAAACCAAAGATTCATAGAAATCACTGCCACGCGCGATTAACTCCAGTCCGa                                                                                                                                                                                                                <  2:358574/244‑1 (MQ=255)                     tGATTGTC‑GTCGGCCTGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCAGATCCGTATCATTATg                                                                                                                                                                                                  >  1:230916/1‑251 (MQ=255)                                     tGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTT‑AAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCgcgc                                                                                                                                                                                  >  1:68876/1‑251 (MQ=255)                                     tGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGATGTGGACGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGATTAAGTGGATTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGAATGAGCTGGTAGGTACAACAGATGCGGAAACCAAAGATTCAGAGCAATCAATGCCACGCGCGATTAAATCCAATCCGATCCGTATCATTATTTTCTAAGTCTTagcg                                                                                                                                                                                   >  2:382899/1‑251 (MQ=255)                                                                 cACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACCGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACgacg                                                                                                                                                        >  1:10978/1‑248 (MQ=255)                                                                 cACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGcc                                                                                                                                                         <  1:148702/250‑1 (MQ=255)                                                                      aCTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTTGGGATTGAGCTGGTAGGTACAACAGCTGCGGCAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTgg                                                                                                                                                     >  1:125218/1‑249 (MQ=255)                                                                                   cGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTTGGTAGTCCCg                                                                                                                                      >  2:121436/1‑251 (MQ=255)                                                                                               ttttGTGGAAAGACGGCGGCTGGTTCCCGAAAGGTTTAAGGGGCTTCTTTGCAGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTTAGCTGGTAGGTACAACGGCTGCGGAAACCACAGATTCAGAGGAATCACTGACACGCGCGGTTAAATCCATTCCGATCCGACTCATTATGTTCTACGTCTTCACGCTGCTTGTGATTATTTGCGAGACGCCGTGGAAGTCGGTAGTCCCGGGGaaaatccgg                                                                                                                          >  2:370708/2‑246 (MQ=255)                                                                                               aTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTc                                                                                                                                                                                                                   <  1:304030/162‑1 (MQ=255)                                                                                               aTTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGGAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTc                                                                                                                                                                                                                   >  2:304030/1‑162 (MQ=255)                                                                                                tttGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAAAAGCTGCGGAAACCAAAGATTCAGAGAAATAACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGt                                                                                                                         >  1:305771/1‑251 (MQ=255)                                                                                                   gtgGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGtttt                                                                                                                      >  1:92990/1‑250 (MQ=255)                                                                                                             ggcggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCTGCAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGAATCATTATGTTCTACGTCTTAGCGCTGATTGTGCTTATGTCCGTGACGCCGTGTATTTCGTTAGTACCGGAGAAAGACCCGTTTGTTGAACTGtt                                                                                                            >  2:206478/1‑251 (MQ=255)                                                                                                                                gTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACCGCTGCGGCAACAACAGATTCAGAGACATCACTGCCACGCGCGATTAACTCCATTCTGATCCGTATCATTATGTTCTACGTCTTCGCGCTGCTTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGctgcct                                                                                         >  1:196089/1‑251 (MQ=255)                                                                                                                                gTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAACGCCCGTTTGTTGAACTTGTCGTGTTGGTAGGGatgcat                                                                                         >  2:94363/1‑249 (MQ=255)                                                                                                                                                     gATTCCAGATAGCGGTTTTCGCTTTCGTGGG‑ATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGCGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGGGATCAACt                                                                   >  2:168884/1‑251 (MQ=255)                                                                                                                                                          cAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCACGGAGAAAAGCCCGTTTGTTGCACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTtgt                                                               >  1:24883/1‑251 (MQ=255)                                                                                                                                                            gATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTgtgg                                                             >  1:115648/1‑251 (MQ=255)                                                                                                                                                             ttAGCGGTTTTCTTTTTCGTGGGGATTGTGCTGGTAGGTTCATCAGCTGCGGAACCCAAAGATTCATAGAAATTACTGCCTCGCTCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCCCGCTGATTGTTATTATGTCCTTGACTCCGTGGTGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTTGGAGGGCTGCCCGCTTCCTCATTCTTGATCAACTTTgtgg                                                             <  2:230916/249‑1 (MQ=255)                                                                                                                                                                        cgttttgttGGGGATTGAGCTGGTAGGTACAACAGCTGGGGAAACCAAAGATTCAGAGAAATCTCTGCCCCGCGCGATTAACTTCATTCCGAGCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTTACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGACCTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACctct                                                  <  2:335305/242‑1 (MQ=255)                                                                                                                                                                         gCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTgc                                                >  1:290334/1‑251 (MQ=255)                                                                                                                                                                                                                   aCCAAAGATTCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTAt      <  1:121436/251‑1 (MQ=255)                                                                                                                                                                                                                       aaGATTCAGAGAAATCCCTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTg  <  1:357328/251‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      AATCGTCGCCATTGTGTCGCTGATTGTC‑GTCGGCCTGGTCATGGTGGCGATGCACT‑TTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATGTCCGTGACGCCGTGGAGTTCGGTAGTCCCGGAGAAAAGCCCGTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTG  >  NC_000913/4430367‑4430829

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