Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,704,167	G→A	P48S (CCA→TCA)	ydgJ ←	putative oxidoreductase YdgJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,704,167	0	G	A	96.5%	81.9 / NA	29	P48S (CCA→TCA)	ydgJ	putative oxidoreductase YdgJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/7);  minor base C (0/1);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCAACCCAACACGGATGTTGTCG  >  NC_000913/1704033‑1704303                                                                                                                                       |                                                                                                                                         gACCACATGTTTCCCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTCCGGCTCAGAGACAACCGTAACCGTTGACCAGt                                                                                                                                     <  1:691867/140‑1 (MQ=255)      caTGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCt                                                                                                                                >  2:97633/1‑140 (MQ=255)              cccGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACttt                                                                                                                        <  1:97633/140‑1 (MQ=255)              cccGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACAGTACCCGTTGACCCGTCGGCTTTTACttt                                                                                                                        >  1:283250/1‑140 (MQ=255)                 gCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTAGGCTTTTACTTTTGt                                                                                                                     >  1:244484/1‑140 (MQ=255)                  ccTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAGCCGTTGACCAGTCGGCTTTTACTTTTGtt                                                                                                                    >  1:315603/1‑140 (MQ=255)                                        cGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCg                                                                                              >  2:123596/1‑140 (MQ=255)                                         gggAAATGGGTATCGTGGGGTGTAGGAATGCTAGTCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGCCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGc                                                                                             <  2:347211/140‑1 (MQ=255)                                                   tATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACAGCAGGAGATGACCGCCAGTTCCTGc                                                                                   >  1:305055/1‑140 (MQ=255)                                                          gggTGTAGGAATGACAACCAGGTCTATGTTGGGATCGTTTAACAGATGTTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCg                                                                            <  2:315603/140‑1 (MQ=255)                                                          gggTGTAGGAATGACAAACAGGTCTATGTTGGGATTGTTAAACAGATGTTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCg                                                                            <  2:301771/140‑1 (MQ=255)                                                           ggTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGCTTGCCGCCAGTTCCTGCCCGGGCGt                                                                           >  1:301771/1‑140 (MQ=255)                                                               tAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGccc                                                                       >  2:484084/1‑140 (MQ=255)                                                                       aCAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAg                                                               >  2:684332/1‑140 (MQ=255)                                                                            cAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGATGGAGATTCCCGCCAGATCCTGCCCGGGCGGGCCCGCAATCAGGGgcg                                                          >  2:241253/1‑140 (MQ=255)                                                                                        gggATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGttt                                              >  1:5574/1‑140 (MQ=255)                                                                                         ggATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACTGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGtttt                                             >  2:407600/1‑140 (MQ=255)                                                                                         ggATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGtttt                                             >  2:118094/1‑140 (MQ=255)                                                                                         ggATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGtttt                                             >  2:144627/1‑140 (MQ=255)                                                                                         ggATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGtttt                                             >  2:712710/1‑140 (MQ=255)                                                                                          agtCGTTATACAGATGCTTCGGCTCAGAGAATCTCGTAAGCGTTGACCAGTGGGCTTTTACTTTTTTTTCATCACTGCTGGAGTTTCCCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTg                                            <  2:358087/138‑1 (MQ=255)                                                                                            tCGTTAAACAGATGCTTGGGCTCAGAGACAATCGTGACCGTTGACCAGCCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCt                                          <  2:305055/140‑1 (MQ=255)                                                                                                       aTGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCGATCAGGGACGCATGAAAGGTTGTGCTCGCATAAccct                               >  1:358087/1‑138 (MQ=255)                                                                                                       aTGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCAt                               >  1:633020/1‑140 (MQ=255)                                                                                                       aTGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCAt                               >  2:392134/1‑140 (MQ=255)                                                                                                        tGCTTCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATa                              >  1:435010/1‑140 (MQ=255)                                                                                                            tCGGCTCAGAGACAACCTTACCCGTTGACCATTCGGCTTTTACTTTTGCTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCa                          <  2:244484/140‑1 (MQ=255)                                                                                                            tCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCa                                <  1:61241/134‑1 (MQ=255)                                                                                                            tCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCa                                >  2:61241/1‑134 (MQ=255)                                                                                                            tCGGCTCAGAGACAACCGTAACCGTTGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCa                          >  1:89599/1‑140 (MQ=255)                                                                                                                                    ttGACCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCAACCCAACACGGATGTTGTCg  <  1:684332/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         GACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCAACCCAACACGGATGTTGTCG  >  NC_000913/1704033‑1704303

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