Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,032,832	C→T	R6C (CGC→TGC)	yqfG →	uncharacterized protein YqfG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,032,832	0	C	T	96.9%	88.4 / NA	32	R6C (CGC→TGC)	yqfG	uncharacterized protein YqfG
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (15/16);  total (16/16)
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.92e-01

TCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGCGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATCCGCTCAAAACGGGCGGTTGTCGATAAACGCTCACTTGGTTAATCATTTCACTCTTCAATTATCTATAATGATGAGTGATCAGAATTACATGTGAGAAATTATGCAAACGG  >  NC_000913/3032601‑3033074                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    tcttttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATTGCTATTAGCATCATGCTGAGCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCtgt                                                                                                                                                                                                                                   <  2:163171/250‑1 (MQ=255)          aGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATCAAGCTTTAAAAGTAAGGACAGTACGTATGTGGATTCTGATGCGCAAGCTGAAGAAAACTGAGCCTGGAGACTAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCtttt                                                                                                                                                                                                                         >  2:439436/1‑251 (MQ=255)          aGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAAGTTTTTAAGGTGCGCTATATTCAGTCTGTTGTTGCtttt                                                                                                                                                                                                                         <  1:41885/251‑1 (MQ=255)           gTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCAAGCTGACCTGGTGAAGGAATCCGTATCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCttttt                                                                                                                                                                                                                        >  1:446415/1‑251 (MQ=255)                                  tgctACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTc                                                                                                                                                                                                 >  2:153533/1‑251 (MQ=255)                                       cATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGAAAGTGCGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATGTTTTAATGTGCGCTATATTCAGACTGTTGTTGCTTTTTACTCCCTCTATTTCTGGCATTTCTTACt                                                                                                                                                                                            <  2:116270/251‑1 (MQ=255)                                             ggCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTg                                                                                                                                                                                      >  1:421292/1‑251 (MQ=255)                                                 ttATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGATATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATACAGATTTAAAAGTCAGGACAGTACGGATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAAt                                                                                                                                                                                  >  1:140596/1‑251 (MQ=255)                                                                 tAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGGTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTATTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGc                                                                                                                                                                  >  2:299891/1‑251 (MQ=255)                                                                       cGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCaaa                                                                                                                                                            <  2:306910/251‑1 (MQ=255)                                                                       cGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCaa                                                                                                                                                             >  1:110538/1‑250 (MQ=255)                                                                           gtgCTGAACGACTTCTTTAATGGCTATCAGCATCATGCTGACCGGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAatatt                                                                                                                                                        >  2:123871/1‑250 (MQ=255)                                                                              cTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGACGGCCTCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGACAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCt                                                                                                                                                                                                                                     >  1:200391/1‑171 (MQ=255)                                                                              cTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCt                                                                                                                                                                                                                                     <  2:200391/171‑1 (MQ=255)                                                                                gAACGTCTGCTTTAATGGCTAGCAGCATCATGCTGGCTTGGTGAAGGAATCCGTCTCTGAGAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAGTAATATGAATTGTTGAATGTGCTCTATATTCAGTCGGTTGTTGTTTTTTAGGCTCTCTATGCCTTTCTGTTCTGACGGTGGTGTAAGGGCCATTGAGATTTGCTTCaagtttattata                                                                                                                                                   <  2:280658/251‑11 (MQ=255)                                                                                    gTCTTATTTAATGGATATCAGAATCATGCTGACCTGGTAAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTCCGTATGTGGCTTCTGATGCGCAAGATGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTACTGCGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTATCTATTCATGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCtg                                                                                                                                               >  2:311278/1‑249 (MQ=255)                                                                                           ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTTCGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGGGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACtttttt                                                                                                                                        <  1:139742/250‑1 (MQ=255)                                                                                           ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACtttttt                                                                                                                                        <  1:88052/250‑1 (MQ=255)                                                                                           ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACtttttt                                                                                                                                        <  1:153533/250‑1 (MQ=255)                                                                                                      cAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCAGTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGttttt                                                                                                                             <  1:435075/250‑1 (MQ=255)                                                                                                               gCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAAGTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATGCCTGTCATTTCTGACTGTGTTGCAATGGCCAGTGAAAGTCGCTTCAAATATGTGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGc                                                                                                                   <  1:46514/251‑1 (MQ=255)                                                                                                                  gACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTAAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAACTATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATc                                                                                                                >  2:155245/1‑251 (MQ=255)                                                                                                                   aCCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATcc                                                                                                               <  1:439436/251‑1 (MQ=255)                                                                                                                   aCCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTAGTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATcc                                                                                                               <  1:403868/251‑1 (MQ=255)                                                                                                                               aaTCAGTCTCTGATACAGCTTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGACTAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCACTGGCCATTGAAAGTCGCTTAAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATCCGCTCAAAACggg                                                                                                   >  2:121906/1‑251 (MQ=255)                                                                                                                                 tCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATCCGCTCAAAACGGGCg                                                                                                 >  1:178870/1‑251 (MQ=255)                                                                                                                                    gtcgctGATAAAGACTTAAAAGTCAGGACAGTACGTAGGTGGCGTCTGATGCGTCAGCTGAAGAAAAATGATCAGGTAGAGTAATATTAATTTTTTAATGTGCGCTATATTCAGGTTGTTGTTGCTTGTTACTCTCGTTATTCCTGTTATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCGGCAAATA‑GTGAATGATTCGTTTTGAAAGGGTTGTTGCCCGTCGGCATCCGCTCAAAACGGGCGGta                                                                                              <  2:407514/247‑2 (MQ=255)                                                                                                                                                    tAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAATTATatga‑t                                                                                                                                                   >  2:255287/1‑182 (MQ=255)                                                                                                                                                    tAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATatga‑t                                                                                                                                                   <  1:255287/182‑1 (MQ=255)                                                                                                                                                                cAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGTGCGCTATATTCAGTCTATTGTTGCTTTTTTCTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCATTGGCCTTTGACAGTCGCTTCAAATATATGA‑TGCTACTTTTTTACATGGTTTTTACCTGTCGGCATCCCCTCAAAACGGGCGGTTTTCGATAAACGCTCACTTGGTTAATCAt                                                                  >  1:317671/1‑251 (MQ=255)                                                                                                                                                                                                                            aaTTTTTTAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATCCGCTCAAAACGGGCGGTTGTCGATAAACGCTCACTTGGTTAATCATTTCACTCTTCAATTATCTATAATGATGAGTGATCAGAATTACATGTGAGAAATTATGCaa      <  1:88051/251‑1 (MQ=255)                                                                                                                                                                                                                                 tttAATGTGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATCCGCTCAAAACGGGCGGTTGTCGATAAACGCTCACTTGGTTAATCATTTCACTCTTCAATTATCTATAATGATGAGTGATCAGAATTACATGTGAGAAATTATGCAAACgg  >  1:403263/1‑250 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    TCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTAATGCGCGCTATATTCAGTCTGTTGTTGCTTTTTACTCTCTCTATTCCTGTCATTTCTGACTGTGTTGCAATGGCCATTGAAAGTCGCTTCAAATATATGA‑TGCTACTTTTTTAAATGGTTTTTACCTGTCGGCATCCGCTCAAAACGGGCGGTTGTCGATAAACGCTCACTTGGTTAATCATTTCACTCTTCAATTATCTATAATGATGAGTGATCAGAATTACATGTGAGAAATTATGCAAACGG  >  NC_000913/3032601‑3033074

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