Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,415	C→T	A31A (GCC→GCT)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,415	0	C	T	96.0%	61.7 / ‑5.2	25	A31A (GCC→GCT)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (15/9);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

AATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGT  >  NC_000913/4617172‑4617630                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         aaTTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGtc                                                                                                                                                                                                                  <  2:26497/251‑1 (MQ=255)                    cATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGt                                                                                                                                                                                               >  1:136003/1‑251 (MQ=255)                                      ttGTACGTAGATATCCTTCATTGTGCCGTATCTCGCAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGCATTTTATTGCCGACAAGCCAGGAGAATGAAACGCCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATCGATGGACCTGACCACCCTGAATGACGCCGACACCGACGAGAAAGTAATTGCT‑CTATGTCATCACGCCCAAACCCCGGTCGGCACTACCGCCGCTAt                                                                                                                                                                             >  2:109881/1‑251 (MQ=255)                                       tGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATc                                                                                                                                                                            >  2:23897/1‑251 (MQ=255)                                         tAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTg                                                                                                                                                                          >  2:103105/1‑251 (MQ=255)                                                             gaaggtaaaTCGAAGTGGGTTGCGGTGAAGAATTGAGAAAAATAAAAAACTCGAAAGGTGATTTTGATGGGCGAAAAGCCAGGAAAATGAAATTACTGATCTAAAAGCAAACAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTAAAGAACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATcc                                                                                                                                                      <  1:185664/243‑1 (MQ=255)                                                               aTGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg                                                                                                                                                    >  1:176271/1‑251 (MQ=255)                                                                          ggtgtgtTGCGGAGTAGATGTTAAAAAACTAACAAACTCGCAAGTTAAATTTTATTGGCGACAAGCCAGGAGAATAAAATGACTAATCTGAAAGCAAGCAGCCGGCGTGCATTGAAATTGAGGGACCGGCCCACCCTAAATGACGCCGACACCGACGAGAAAGT‑ATTTCCTCTGTTTCATCAGGCCAAACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCaa                                                                                                                                         <  1:114402/250‑1 (MQ=255)                                                                                          gATGTTAGAATACTAACAAACTCGCAAGGTGATTTTTATTGGCGACAAGCCAGGAGCATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                           >  2:185852/1‑249 (MQ=255)                                                                                                          aaaaCTCGAAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATTACTGATCTGAAAGCAATAAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAATTGTT‑TTGTGTCATCGGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATcc                                                                                                         <  1:21014/250‑1 (MQ=255)                                                                                                                    aGGTGAATTTTCTTGGCGACACGCCAGGAGAATGAAATGACGGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGCCACCGACGAGAATGTAATTGAT‑CTGTGTCATCAGGCACAAACCCCGGTCGACAATACCGCCGCTATCTGTATCAATACTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGCCAATCCGTAttgatac                                                                                               >  2:197886/1‑245 (MQ=255)                                                                                                                    aGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTCAAAGAAAGCAGCCTGCGTGCACTGCAATTGATGGACCTGACCAACCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATAAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTAc                                                                                               >  1:172401/1‑251 (MQ=255)                                                                                                                                        aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                           >  2:134937/1‑251 (MQ=255)                                                                                                                                            cAGGAGAATGAAATGACTGATCTAAAAGCAAGCAGCCGGCGTGCACTAAAATTGATGGACCTGACCCCCCTAAATGACGACGACACCGACGAAAAAGTAATTCCT‑CTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTaa                                                                       <  1:23897/251‑1 (MQ=255)                                                                                                                                                   atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGACCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCCCCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                >  1:130476/1‑251 (MQ=255)                                                                                                                                                              aatcggaaagcaaaaggcatGCGGGCACTGAAATTGATGGACCTGCCCCCCCTAAAAGACGCCACCCCCGACAAGAAAGTAATGGCT‑CTGTGTCAACAGGCAAAACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACATCCCACACGGTAACGATgacatcgaaatc                                                       <  2:136003/230‑1 (MQ=255)                                                                                                                                                                    aaaGCAACAAGCCTGCGTGCCGGGAAATTGATGGACCTGCCCCCCCTGAATGACAACGACACCGCCGAGAAAGTAATTGCT‑CTGTTTCACCAGCCCAAACCCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                               <  2:15109/251‑1 (MQ=255)                                                                                                                                                                           gcagcCTGCGTGCACTGAAATTGCTGGACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGCCATCGCCATCGCGCTTGCAGAAAcc                                        >  1:12427/1‑251 (MQ=255)                                                                                                                                                                                 tGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGCACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                                  >  1:226463/1‑251 (MQ=255)                                                                                                                                                                                       cACTGAAAAGGATGGACCGGCCAACCCGGAATGCCGACGACACCGACTAGAAAGTAATTGCT‑CGGTGCCATCAGGCCAACTCCCGGGGCGCCAAAACCGCCGCTATCTGTATCTATCCGCGCGTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                            <  2:12427/251‑1 (MQ=255)                                                                                                                                                                                                     gACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                          <  1:223821/223‑1 (MQ=255)                                                                                                                                                                                                     gACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                          >  2:223821/1‑223 (MQ=255)                                                                                                                                                                                                     gACCTGACAACCCTGAATGACGACGACACCGAAGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGCCGGCAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGACAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccgg              >  2:163162/1‑248 (MQ=255)                                                                                                                                                                                                          gACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATACGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatga           >  1:41477/1‑249 (MQ=255)                                                                                                                                                                                                                 cgaaaTGACGACAACACCGACGAAAAAGTAATTGCT‑CTGGGTCATCAGGCAAAACCCCCGGTCGGCAAAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGt  <  1:130736/248‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         AATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGT  >  NC_000913/4617172‑4617630

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