Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,127,049	A→G	T82A (ACG→GCG)	yceM →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,127,049	0	A	G	93.6%	143.6 / 0.9	47	T82A (ACG→GCG)	yceM	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (13/31);  total (14/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

GCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAA  >  NC_000913/1126952‑1127144                                                                                                  |                                                                                                gcgcATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCg                                                                                               <  1:936065/100‑1 (MQ=255)  cgcATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGt                                                                                              <  2:349050/100‑1 (MQ=255)   gcATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGtt                                                                                             <  2:381250/100‑1 (MQ=255)     aTTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTAc                                                                                           <  1:267960/100‑1 (MQ=255)              gCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCgg                                                                                  >  1:1116689/1‑100 (MQ=255)               ccGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCggg                                                                                 <  2:882342/100‑1 (MQ=255)                 gATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCggggg                                                                               <  2:120806/100‑1 (MQ=255)                    gatctATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTAc                                                                            <  1:254802/96‑1 (MQ=255)                      gTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCAACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACAt                                                                          >  2:549742/1‑100 (MQ=255)                        tATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGt                                                                        >  1:365156/1‑100 (MQ=255)                         aTCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACagatc                                                                       >  2:254799/1‑96 (MQ=255)                         aTCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTc                                                                       >  2:880904/1‑100 (MQ=255)                           ccAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCtg                                                                     >  2:575912/1‑100 (MQ=255)                           ccAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCtg                                                                     <  1:35724/100‑1 (MQ=255)                              gCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTgtgt                                                                  >  1:488746/1‑100 (MQ=255)                               ccTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTc                                                                 <  1:162702/100‑1 (MQ=255)                               ccTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTc                                                                 <  2:659773/100‑1 (MQ=255)                                  tGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGAt                                                              <  1:438674/100‑1 (MQ=255)                                   gccgccAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGAt                                                              <  2:511836/99‑1 (MQ=255)                                      gccAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAAc                                                          <  2:130926/100‑1 (MQ=255)                                      gccAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAAc                                                          <  1:898828/100‑1 (MQ=255)                                      gccAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAAc                                                          <  1:793210/100‑1 (MQ=255)                                           ttGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTg                                                     <  1:124525/100‑1 (MQ=255)                                             gCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGc                                                   <  2:350958/100‑1 (MQ=255)                                                        tttGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTg                                         <  1:1004202/99‑1 (MQ=255)                                                        tttGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTg                                         <  2:427343/99‑1 (MQ=255)                                                          tgtgCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgcgc                                      <  2:667252/100‑1 (MQ=255)                                                          tgtgCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgcgc                                      <  2:1156181/100‑1 (MQ=255)                                                          tgtgCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgcgc                                      <  1:1210507/100‑1 (MQ=255)                                                             gCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGAt                                   <  1:578735/100‑1 (MQ=255)                                                              cATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATg                                  <  2:145219/100‑1 (MQ=255)                                                              cATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATg                                  <  1:134432/100‑1 (MQ=255)                                                              cATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATg                                  >  2:361686/1‑100 (MQ=255)                                                               aTTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGc                                 <  1:1204523/100‑1 (MQ=255)                                                                 tCCAGCACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTg                               <  2:496555/100‑1 (MQ=255)                                                                  ccAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGa                              >  1:507038/1‑100 (MQ=255)                                                                      cACCGCCAGCCACTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACgg                          <  1:915747/100‑1 (MQ=255)                                                                          gccagccaCTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCtgg                      <  2:1172832/100‑1 (MQ=255)                                                                            cagccaCTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCtggtg                    <  2:714707/100‑1 (MQ=255)                                                                             agccaCTTTGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCtggtgg                   >  1:265496/1‑100 (MQ=255)                                                                                     tGACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTggcg           <  2:752347/100‑1 (MQ=255)                                                                                      gACGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTggcgg          >  1:862397/1‑100 (MQ=255)                                                                                        cGTGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGgcggcg        >  2:1071836/1‑100 (MQ=255)                                                                                          tGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGcgcg      <  2:861206/100‑1 (MQ=255)                                                                                          tGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGcgcg      >  1:852706/1‑100 (MQ=255)                                                                                          tGGTCAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGcgcg      <  2:961027/100‑1 (MQ=255)                                                                                              cAGTGCGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGaaa  >  2:1029196/1‑100 (MQ=255)                                                                                                  |                                                                                                GCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAA  >  NC_000913/1126952‑1127144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑