Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,127,334	G→A	G177S (GGC→AGC)	yceM →	putative oxidoreductase YceM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,127,334	0	G	A	93.5%	78.0 / NA	32	G177S (GGC→AGC)	yceM	putative oxidoreductase YceM
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (20/9);  total (22/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTG  >  NC_000913/1127101‑1127575                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  aTCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCACTAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTg                                                                                                                                                                                                                                  >  1:349956/1‑251 (MQ=255)     gcgggATGCTGAACAGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTg                                                                                                                                                                                                                                  <  2:51705/243‑1 (MQ=255)               taccggttgggggAACTGGGGGCGCGGAAAAACCTGCCGTTGATTGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCCCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTAc                                                                                                                                                                                                                       <  2:179212/237‑1 (MQ=255)                          ggAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTAGGGCCACACGATCTTTATTTCACGTGGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGAGGGTGTCGAGCTACAAAGCCTCTATGGATGGCGGTACGCTACTGACTAacg                                                                                                                                                                                                         >  1:54455/1‑251 (MQ=255)                                     cgtggAACAAAATGACGCTGATGGTAGGTTTTAACCGTCGTTACGCACCAATCGACGGTGAGTTAAAAACGCAACTCGCCACAGCAGCCTCGCTAAGAATGGCTATACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATCTCACGTGGCTGGATGATTATCTGCATGTGGTGGCAACCGCGACGTGGTTGTCGGGGTACCAAGCCTCTCTGGTTGGCGGTATGCCACTGAATATAGCCGCTGGCGa                                                                                                                                                                                               >  2:167080/4‑250 (MQ=255)                                        ggAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATg                                                                                                                                                                                           >  1:278071/1‑251 (MQ=255)                                               aCTGACGTTGATGGTCGGTTTTAACCGTTGTTTCGCTCCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCTCACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTg                                                                                                                                                                                    <  1:312597/251‑1 (MQ=255)                                                cTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCCCCACTCGAAGGTAAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTACGAATGGATAAACATCGGAGCAATAGCGTCGGACCACACGATCTTTCTTTCACTTTGCTGGATGATTATCTGCATGTGGTGGATACGGCGCTGTGGATTTCGAGAGACAAACCCTATCTGGAAGGCGGGACGCTACTGACCGACGCCGCTGAGGAAATGCTGTTTgt                                                                                                                                                                                   >  2:281110/1‑250 (MQ=255)                                                       tGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCAc                                                                                                                                                                            >  1:57879/1‑251 (MQ=255)                                                                    tAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTCGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTTCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTCCGCTACTGACTAACGGCGCTGGCGAAAt                                                                                                                                                                                            >  2:249588/1‑222 (MQ=255)                                                                    tAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAAt                                                                                                                                                                                            <  1:249588/222‑1 (MQ=255)                                                                         gtcgtTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGCTTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCt                                                                                                                                                          >  1:57151/1‑251 (MQ=255)                                                                         gtcgtTTCGCACAACTCTAAGGTGAGTTACAAACGCAACTCGCCACCTCAGCCTCGCTAAGAATGGATAAACATCGAAGCAATAGCGTCGGGCAACACGATCTTTATTTCACGTTGCTGGATGCTTGTCTGAATGTGGTGGATCCCGCTCTGTGATTGGAGAGCGACAAACCCTCTCTGGATGGCGGTACGCTACTGACTAGCGACGCTAGCGACATTCTGTTTGACGTGCACCATTTTTCGGCTGGtact                                                                                                                                                          >  2:76161/1‑248 (MQ=255)                                                                                              gTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAAATCACCACCTGTAt                                                                                                                                      >  1:79674/1‑250 (MQ=255)                                                                                                 aGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACCCGCTCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCAt                                                                                                                                  >  1:186264/1‑251 (MQ=255)                                                                                                   ttAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCg                                                                                                                                >  1:157920/1‑251 (MQ=255)                                                                                                                            ccTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAAcc                                                                                                       >  1:357838/1‑251 (MQ=255)                                                                                                                                 cTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGc                                                                                                   >  1:101090/1‑250 (MQ=255)                                                                                                                                         ggttAAACATCGTAGCAATAGCGTCGGGCCACAAGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGAGGTGGATACCGCGCTGTGGTTGTAGAGCGACAAAGCCACTATGGATGGCGGTACGCTACTGACTAACGACGCTGGAGAGATGCTGTTTGCCGAGCACAATTTTTCGGATGGACCATTGCAGATCACCACCTGCATGCATCGCCGTGTCGGAGGTCGGCGTGAAAcaggcaaggacgtgg                                                                                          >  2:89212/4‑237 (MQ=255)                                                                                                                                          gctAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGCTCCCGCGCTGTGGTTGTCGAGCGCCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCCCCACCTGTATGCATCTCCGGTACGGATGTCTGCTTGACACCGTGCAGGACGttac                                                                                         >  2:173190/3‑248 (MQ=255)                                                                                                                                                              cGTCGGGCCACACGATCTTTATTTCACGTTGATGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGACACCGTGCAGGCCGTGACTGCAGGGGCGCTCATCGAc                                                                      >  1:99951/1‑250 (MQ=255)                                                                                                                                                                 cGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAAACGTGCAGGCCGTGACTGACGGTGCGCTCATCGACAtt                                                                   >  1:17804/1‑250 (MQ=255)                                                                                                                                                                                     tCACGTTCCGGGATGATTATCTGAAGGTGGGGGATACCTCGCTGGGGTTGTCGATCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTACCGACGCTGGCGAAATGCGGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCTTGAAACCGTGCAGGCCGTGACTGGCGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCg                                               <  2:350168/250‑1 (MQ=255)                                                                                                                                                                                             cTGGATGATTATCTGCATGTGGTGGATACAGCGCTGTGGTTGTCGCGCGACAAAGCCTCTCTGGATGGCGGTACGCGACTGACTAACGACGCTGGAGAAATGCTGTTTGCCGAGCCCCATTTTTCGGCTGGTCCTTTGCAGATCACCCCCTGTATGCCTCGCCGTGCCGGAAGTCAGCGTTAACCCGTTCAGGCCTTGACTGACGGTGCGGTCCTCTACATTCAGGATATGCGCGCATGGCGTcaggcacg                                      >  2:328804/1‑247 (MQ=255)                                                                                                                                                                                               ggATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAgcgc                                     <  1:269489/250‑1 (MQ=255)                                                                                                                                                                                                      tATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAggg                             >  1:332271/1‑251 (MQ=255)                                                                                                                                                                                                           gcctgtgGTGGATACCGCGCTGGGGTTGTCGAGCGACAAAGCCTTCCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAg                        <  2:307351/248‑1 (MQ=255)                                                                                                                                                                                                           gCATGTGTTGGATTCCGTGTTGTGTTTGGCGGGCGACAAAGCCTTGCTTGATGGCAGTTCGCTTCTGACTTACGACGCTGGCGAAATGCTGTTTGCCGGGTACCCGTTTTCGGCTGGTCCTTTGCATATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGATACCGTTCAGGCCGTGACTGACGGTGCGCTCATCTACAGTACGGATATGCTCGGATGGCGTGAGGAGCGCGGGCAGGGCGTAg                        <  2:357838/251‑1 (MQ=255)                                                                                                                                                                                                                   tggATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATaa                 >  1:360226/1‑250 (MQ=255)                                                                                                                                                                                                                    ggagaCCGCGCTGTGGGTGACGAGCGACAAAGCCTCTCTGGATTGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATaaa                <  2:101090/246‑1 (MQ=255)                                                                                                                                                                                                                               tgtgGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCCCCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGt    >  1:227664/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTg  <  1:285153/250‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  ATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTG  >  NC_000913/1127101‑1127575

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