Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,781,964	C→G	G190G (GGC→GGG)	ydiS →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,781,964	0	C	G	98.0%	170.4 / NA	50	G190G (GGC→GGG)	ydiS	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (24/25);  minor base A (0/1);  total (24/26)
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.97e-01

AGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGCCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAGGATTCCATATCCT  >  NC_000913/1781828‑1782097                                                                                                                                         |                                                                                                                                      aGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCt                                                                                                                                     <  2:825135/139‑1 (MQ=255) aGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCt                                                                                                                                     <  2:222377/139‑1 (MQ=255)     tGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTcaca                                                                                                                                 >  2:919081/1‑139 (MQ=255)           tgatATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGaa                                                                                                                           <  1:555155/139‑1 (MQ=255)               atTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCt                                                                                                                                     >  2:816878/1‑125 (MQ=255)               atTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCt                                                                                                                                     <  1:816878/125‑1 (MQ=255)               atTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGa                                                                                                                       >  1:791108/1‑139 (MQ=255)                 tctcGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCt                                                                                                                                     <  1:183917/123‑1 (MQ=255)                 tctcGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATc                                                                                                                     >  1:717675/1‑139 (MQ=255)                   tcGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCaa                                                                                                                   >  2:729381/1‑139 (MQ=255)                       aGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               <  1:919081/139‑1 (MQ=255)                       aGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               >  2:429561/1‑139 (MQ=255)                       aGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               <  1:429561/139‑1 (MQ=255)                             tgtgGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCttt                                                                                                         >  1:630353/1‑139 (MQ=255)                             tgtgGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCttt                                                                                                         <  2:791108/139‑1 (MQ=255)                                 gtgATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGG‑TATTGGGCTCACACCAGAACAGATCAACGATCGCttt                                                                                                         >  2:430641/1‑134 (MQ=255)                                 gtgATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGG‑TATTGGGCTCACACCAGAACAGATCAACGATCGCttt                                                                                                         <  1:430641/134‑1 (MQ=255)                                              aTGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               >  1:56506/1‑116 (MQ=255)                                              aTGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               <  2:56506/116‑1 (MQ=255)                                                    ttAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               >  2:248154/1‑110 (MQ=255)                                                    ttAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               <  1:248154/110‑1 (MQ=255)                                                                 tGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCacac                                                                                                                                >  1:460528/1‑80 (MQ=255)                                                                 tGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCacac                                                                                                                                <  2:460528/80‑1 (MQ=255)                                                                       ctcgctGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGt                                                               >  1:779655/1‑139 (MQ=255)                                                                       ctcgctGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGt                                                               >  1:711559/1‑139 (MQ=255)                                                                        tcgctGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGtt                                                              >  2:611584/1‑139 (MQ=255)                                                                        tcgctGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGtt                                                              <  1:611584/139‑1 (MQ=255)                                                                            tgtgAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACTGGCGAGGAAGGTGCCGCCTGGCTGTGTGCCGGTTcccc                                                          <  2:1009833/136‑1 (MQ=255)                                                                              ggAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTTCCCtt                                                        <  1:1002235/139‑1 (MQ=255)                                                                                   ggTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGAc                                                   <  2:711559/139‑1 (MQ=255)                                                                                              ccGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATggg                                        <  2:717675/139‑1 (MQ=255)                                                                                               cGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGc                                       >  2:41235/1‑139 (MQ=255)                                                                                                 aTCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACg                                                                                                <  1:297142/80‑1 (MQ=255)                                                                                                 aTCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACg                                                                                                >  2:297142/1‑80 (MQ=255)                                                                                                     gCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGa                                 <  1:121948/139‑1 (MQ=255)                                                                                                      catcatTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCACCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGAt                                >  1:456657/1‑139 (MQ=255)                                                                                                         catTACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTtc                             <  1:41235/139‑1 (MQ=255)                                                                                                            tACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               >  1:264338/1‑54 (MQ=255)                                                                                                            tACGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGAt                                                                                                               <  2:264338/54‑1 (MQ=255)                                                                                                              cGCTGTTGGTGTTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCtat                        >  1:196572/1‑139 (MQ=255)                                                                                                                   ttGGTGTTAAAGAGGTTATTGGGCTCACACCAGaa                                                                                                                           >  2:1062418/1‑35 (MQ=38)                                                                                                                   ttGGTGTTAAAGAGGTTATTGGGCTCACACCAGaa                                                                                                                           <  1:1062418/35‑1 (MQ=38)                                                                                                                        gtTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAgg              >  2:463372/1‑139 (MQ=255)                                                                                                                        gtTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAgg              >  2:586929/1‑139 (MQ=255)                                                                                                                        gtTAAAGAGGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAgg              <  2:6085/139‑1 (MQ=255)                                                                                                                               agGTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAGGATTCCat       <  2:779655/139‑1 (MQ=255)                                                                                                                                 gTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTccc                                                           >  1:797735/1‑85 (MQ=255)                                                                                                                                 gTTATTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAGGATTCCatat     >  1:785814/1‑139 (MQ=255)                                                                                                                                 gTTATTGGACTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTccc                                                           <  2:797735/85‑1 (MQ=255)                                                                                                                                    aTTGGGCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAGGATTCCATATCCt  <  1:586929/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      AGGCTGGGGATGATATTCTCGAAGCGAATGTGGTGATTCTGGCTGATGGCGTTAACTCGATGCTTGGCCGCTCGCTGGGAATGGTTCCCGCTTCCGATCCGCATCATTACGCTGTTGGTGTTAAAGAGGTTATTGGCCTCACACCAGAACAGATCAACGATCGCTTTAATATTACGGGCGAGGAAGGTGCCGCCTGGCTGTTTGCCGGTTCCCCTTCTGACGGCCTGATGGGCGGGGGATTTCTCTATACCAACAAGGATTCCATATCCT  >  NC_000913/1781828‑1782097

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