Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,848,406	T→C	T91A (ACG→GCG)	ydjA ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,848,406	0	T	C	95.2%	140.1 / ‑2.0	42	T91A (ACG→GCG)	ydjA	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/21);  total (20/22)
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

ATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAA  >  NC_000913/1848271‑1848535                                                                                                                                        |                                                                                                                                  aTGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGcgatg                                                                                                                               <  2:330606/140‑1 (MQ=255)  tGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGcgatga                                                                                                                              >  2:686376/1‑140 (MQ=255)  tGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGcgatga                                                                                                                              >  2:271208/1‑140 (MQ=255)      ccACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAg                                                                                                                          >  1:1018262/1‑140 (MQ=255)      ccACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAg                                                                                                                          <  1:686376/140‑1 (MQ=255)        aCTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCg                                                                                                                        <  2:141790/140‑1 (MQ=255)             gcCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTgcg                                                                                                                   <  1:574693/140‑1 (MQ=255)                 gATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGa                                                                                                               >  2:250790/1‑140 (MQ=255)                  aTGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGaa                                                                                                              <  2:1195321/140‑1 (MQ=255)                     ccgccAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACgg                                                                                                           <  2:1018262/140‑1 (MQ=255)                      cgccAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGgc                                                                                                          >  2:638611/1‑140 (MQ=255)                      cgccAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGgc                                                                                                          >  2:952409/1‑140 (MQ=255)                                cTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACggg                                                                                                >  1:1218068/1‑140 (MQ=255)                                  cggaaaCTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCt                                                                                              <  1:266948/136‑1 (MQ=255)                                        ctgctgCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCg                                                                                        >  1:405215/1‑140 (MQ=255)                                                aTTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCttt                                                                                <  1:638611/140‑1 (MQ=255)                                                aTTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCAGAACGGCGCATTACGGGCTTTGTCGATAGCttt                                                                                <  1:765547/140‑1 (MQ=255)                                                    gCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtca                                                                            <  2:384401/140‑1 (MQ=255)                                                          ccATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTa                                                                      >  2:1086485/1‑140 (MQ=255)                                                              gACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCgg                                                                  >  1:1030626/1‑140 (MQ=255)                                                                   cgcATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGGTTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGgcagca                                                             >  2:1098926/1‑140 (MQ=255)                                                                       tCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCg                                                         <  1:447722/140‑1 (MQ=255)                                                                           gCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcccc                                                     <  1:952409/140‑1 (MQ=255)                                                                           gCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcccc                                                     <  2:1030626/140‑1 (MQ=255)                                                                           gCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcccc                                                     <  1:268260/140‑1 (MQ=255)                                                                           gCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcccc                                                     <  1:271208/140‑1 (MQ=255)                                                                            cAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccccc                                                    <  2:1218068/140‑1 (MQ=255)                                                                            cAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTCCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccccc                                                    <  1:1078193/140‑1 (MQ=255)                                                                            cAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCACATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccccc                                                    <  1:445696/140‑1 (MQ=255)                                                                                  tttCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTcc                                              >  2:582451/1‑140 (MQ=255)                                                                                             cAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCt                                   >  2:867891/1‑140 (MQ=255)                                                                                             cAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCt                                   >  2:321514/1‑140 (MQ=255)                                                                                              aGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccccc                                                    >  2:154772/1‑122 (MQ=255)                                                                                              aGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccccc                                                    <  1:154772/122‑1 (MQ=255)                                                                                              aGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTa                                          <  2:1222122/132‑1 (MQ=255)                                                                                              aGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTa                                          >  1:1222122/1‑132 (MQ=255)                                                                                                cgcgGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGaa                                >  1:857160/1‑140 (MQ=255)                                                                                                     gACTTTATGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCt                           >  1:590883/1‑140 (MQ=255)                                                                                                            tGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGccc                    <  2:590883/140‑1 (MQ=255)                                                                                                            tGATTCTCTTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGccc                    <  2:857160/140‑1 (MQ=255)                                                                                                                 ctctTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGccc                    <  1:170571/135‑1 (MQ=255)                                                                                                                 ctctTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGccc                    >  2:170571/1‑135 (MQ=255)                                                                                                                 ctctTCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTGTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttccc               >  2:763745/1‑140 (MQ=255)                                                                                                                    ttCGCATTTCGCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttcccctt            >  1:329613/1‑140 (MQ=255)                                                                                                                              gCCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCCCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACaaa  >  2:855324/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                  ATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAA  >  NC_000913/1848271‑1848535

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