Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,102	T→C	G123G (GGT→GGC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,102	0	T	C	81.8%	59.4 / 5.8	34	G123G (GGT→GGC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (12/15);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

GAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAT‑GCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAA  >  NC_000913/3874‑4341                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGACCAAAAGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTgg                                                                                                                                                                                                                            >  1:236159/1‑250 (MQ=255)                     gctgTGCGAATATCCTCTCGGCGTGTATTGTTTATGAAACCCCGCTGGAAATCCTGGAATAGCGCGTGCGCTCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTTTTCCACGTGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCATTGCCCATTCTGACCGCGTCCTCCGGTGATAcc                                                                                                                                                                                                        <  2:141496/248‑1 (MQ=255)                      cagtgtgaagattcgtgcggggTTTATTGGGGATGATACCCCGGAGGAATGCCTGGGAGAGCGCGTGCGCGCTGCGTTTGCTTGTCCGGCACCGGTGGCCAATGGTGAAAGCGATGTCGGTTGTCTGTAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCGCAAATGCTGACCCATATGGCGTGCGATAATCCATTGACCAGTCGGACCTGGACCTCCTGTTATACCg                                                                                                                                                                                                       <  2:245948/230‑1 (MQ=255)                          gCGAAGACCCTCTCTGCGTTGATGGGTGATGAAAGCCCGCAGGAAATCCTGGAAGAGCTCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTGGTCTGGAATGGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCTACCTCCGGTGATACCGGAg                                                                                                                                                                                                    <  2:33102/250‑1 (MQ=255)                           cGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCg                                                                                                                                                                                                  <  1:325085/251‑1 (MQ=255)                                aTCCTCGCGGCGTTTAGTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCAGATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCGCCGGTGATACCGGAGCGGCAGt                                                                                                                                                                                             <  2:149506/251‑1 (MQ=255)                                       cggggTTTTTTGGTGACGAAAGCCCGCAGGAACTTCTGGAAGAGCTCGTGCGCGCGGCGTTTGCTGTCCCGGCTCCGGGCGCCAATGTTGAAAGCTATGTCGGTTGTCTGGAATTGTTCCTCGGGCCAACGCTGGCATTTAAAGATTTCTGCGGGCTCTTTATGGCACAAATGCTGACCCATGTTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATCCCGGAGCGGCAGTGGCTCAt                                                                                                                                                                                      <  2:137976/247‑1 (MQ=255)                                       cGGCGTTTATTGGTTATGAAATCCCACAGGAAATCCTGGAAGAGCGCGGGCTCGCGGCGTTTGCCTTCCCGGCTTCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCAt                                                                                                                                                                                      <  2:209003/251‑1 (MQ=255)                                           gTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAACTGCTGACCCATATTGCGGGCGATAAGCCATTGACCATTCTTAACGCGACCTCCGGTGATACCGGAGCGGCAGTTGCTCATGCtt                                                                                                                                                                                  >  1:275149/1‑251 (MQ=255)                                             ttATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCCTATTGCGGGTGATCAGCCCGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCCCATGCTTTc                                                                                                                                                                                >  1:132830/1‑251 (MQ=255)                                                ttGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGGGGCTCATGCTTTCTAc                                                                                                                                                                             >  1:308764/1‑251 (MQ=255)                                                      atgaAATCCCACAGGAAATCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTa                                                                                                                                                                       <  2:132830/251‑1 (MQ=255)                                                                     aaaTCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGGGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAggg                                                                                                                                                        >  1:204166/1‑249 (MQ=255)                                                                      aaTCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTgg                                                                                                                                                       <  1:108166/251‑1 (MQ=255)                                                                                                   ttGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAg                                                                                                                           <  1:50316/250‑1 (MQ=255)                                                                                                   gtcgccctgcctagctggTCGCCAATGTGGAAAGTGATGACGGTTGGCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATATCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGt                                                                                                                          <  2:186107/235‑1 (MQ=255)                                                                                                      ccTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTTAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTc                                                                                                                         <  1:199774/249‑1 (MQ=255)                                                                                                      ccTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATATAGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGCCCCATCTTGCGGGGGATAAGCCAGTGACCATTCTGACCGCCAACTCCGGGGGTACCGGAGCGGCAGGGGATAAAGCCTTCTTCGGTTTAACGGATGGTAAAGTGGTTAACCTCCAACCACGGGGCAAAATCAggca                                                                                                                        >  2:199774/1‑247 (MQ=255)                                                                                                                   ggTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGGGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGGATGTGAAAGTGGTTATCCTCTATCCACGGGGCAAAATCAGtagataggaaaagacc                                                                                                           >  1:186107/1‑235 (MQ=255)                                                                                                                     tCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTTTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGAGAAAGTGGTTATCCTCTATCCACGAGGGAAAATCAGTCCAATGCAAGgaaaaatg                                                                                                        >  1:30669/1‑248 (MQ=255)                                                                                                                              ttGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCAGTTAAAGATTTCGGCGGTCGCGTTATTGCCCAAATGCTGACCCCTATTGCGGGCGATAAGCAAGTGAACATTCTGACCGCTACCTCCGGTGATACCGGGACGGCAGGGGCGCCTGCTTTCTACGGTTTACCGAATGTGAAAGGGGTTAACCCATATCCACGGGGCAAAATACGACCACTGCAAACAAAACTGTTCTGTAc                                                                                                >  2:290531/1‑250 (MQ=255)                                                                                                                                gAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGAGGCTCATGCTTTCTACGGATTACCGAATGTGCAAGTGGGTATCCTCTATCCACGAGGCAAAATCAGTACACTGCAAGAAAAACTGTTCTGTAcctt                                                                                             >  1:342634/1‑248 (MQ=255)                                                                                                                                    gCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGggc                                                                                         <  1:121456/251‑1 (MQ=255)                                                                                                                                     cGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCTAGAAAAACTGTTCTGTACATTGggcg                                                                                        <  1:290531/251‑1 (MQ=255)                                                                                                                                           cGGTTGTCTGGATGCGGTCCACGGGCCAAAGCTGGCATATAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATCTTGCGGGCTATAAGCATGTGACCATTCTGTCCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCCCATGCTGTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCTCTGTAAGAAAATCTGTTCTGTACATTTGGCGGCAat                                                                                   <  2:342634/250‑1 (MQ=255)                                                                                                                                               tGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAACCAGTGACCATTCTGACCGAGACCTCCGGGGATACCGGAGCGGCAGTGGCTCATGCATTCTAAGGTTTACCGAATGTGGAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTTTGTACATTGGGCGGCAATATCg                                                                               >  1:380182/1‑250 (MQ=255)                                                                                                                                                                         gCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCg                                                    >  1:185463/1‑251 (MQ=255)                                                                                                                                                                                  tAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAT‑GCCTGTc                                          >  1:285941/1‑251 (MQ=255)                                                                                                                                                                                               ggTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGtt                                                                       >  1:277429/1‑210 (MQ=255)                                                                                                                                                                                               ggTCGCGTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCAAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGtt                                                                       <  2:277429/210‑1 (MQ=255)                                                                                                                                                                                                gTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAT‑GCCTGTCAGGCGCTGGTGacg                            >  1:295107/1‑249 (MQ=255)                                                                                                                                                                                                    cTTTATGACACAAATGCTGACCCATATTGCGGGTGCTAAGCCAGTGCCCATTCTGACCGCGAACTCCGGTGATACCGGAGCGGAAGTGGCTCATGCTATCTCAGGCTTACCGAATGTGAACATGGTTATCCTCTATCCACGCGGCAAAATCAGTCCTATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGACACTGTTGACATCGACGGCGTTTTC‑ATAGCCTGTCAGGCGCTGGAGaatccgg                        >  2:157115/1‑246 (MQ=255)                                                                                                                                                                                                        atggaccaAATGCGGACCCATATTGTGGGCTATAAGCCAGGGACCATTCTGACCGCGGCCTCCGGTGATACCGGGGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAG‑GCCTGTCAGGCGCTGGTGAAGCGGGCGt                     <  2:204166/244‑1 (MQ=255)                                                                                                                                                                                                                          cATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAT‑GCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGaa  <  1:337956/251‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAT‑GCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAA  >  NC_000913/3874‑4341

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: ‑