Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	793,581	G→T	19.8%	L70I (CTC→ATC)	modF ←	ABC family protein ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	793,581	0	G	T	19.8%	45.8 / 3.0	25	L70I (CTC→ATC)	modF	ABC family protein ModF
Reads supporting (aligned to +/- strand):  ref base G (12/8);  new base T (2/3);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGT  >  NC_000913/793440‑793716                                                                                                                                              |                                                                                                                                        tgcaACGCGGTGCATCCTTTACTTAATCCTGAATGATCTCAGCCGTAGTGCGTCAGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTTCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagaga                                                                                                                                        >  1:527918/1‑143 (MQ=255)       gcgGTGCATCCTTTACTTCATCCTGTATGATCCCTCCCTTGGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTg                                                                                                                                  <  2:250171/143‑1 (MQ=255)              aTCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGGGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCCATTGCTCGATGGAGAGACGAGAGAtgtgcg                                                                                                                           >  1:163329/1‑141 (MQ=255)                 cTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAGTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGaaa                                                                                                                        <  2:327423/143‑1 (MQ=255)                       ttCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGCGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAACTGCTCGCAGGCGAGACCACCGAAGAGGGCAAACTgcc                                                                                                                  >  2:243052/1‑141 (MQ=255)                        gcATCCTGAATGAGCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGGGCATATCGGTGTTATTCCGGTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTGGTATGAGATACGAGTGATGTGGGAAAACtcgtc                                                                                                                 <  2:314288/142‑5 (MQ=255)                             cTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGcc                                                                                                            >  1:334132/1‑143 (MQ=255)                                      cgcAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCtt                                                                                                   <  1:298260/141‑1 (MQ=255)                                       tcAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTGTGCCGTGCACCttt                                                                                                  >  2:562731/1‑143 (MQ=255)                                             gTAGTGCGTCCGGTGTCATCTTCGCCATGGCCGAGCAGATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGCAAAACTGGCTTTGCCGTTCACCTTTCAAAAg                                                                                            <  2:378017/143‑1 (MQ=255)                                                 tGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGATACGAGCGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCggc                                                                                        >  2:218526/1‑142 (MQ=255)                                                     tCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTGGAAGGAGAGACGAGTGATGTGGGAAAACGGGCTTTGCCGTTCACCTTTCAAAAGCGGCAGtt                                                                                    >  1:288472/1‑143 (MQ=255)                                                        ggTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGAAATTGCTTGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTccc                                                                                 <  1:509655/143‑1 (MQ=255)                                                            tCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGCCGAGTGATGTGGGAACACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGc                                                                             >  2:306774/1‑143 (MQ=255)                                                            gcATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGATACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGc                                                                             <  1:246368/142‑1 (MQ=255)                                                               tCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGGCGAGGTTTTGCAATTGCTCGAAGGAGATACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAg                                                                          >  1:241745/1‑143 (MQ=255)                                                                cTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAgc                                                                         >  1:338030/1‑143 (MQ=255)                                                                 ttCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAgcg                                                                        >  1:61699/1‑143 (MQ=255)                                                                 ttCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTCCAAAAGCGGAAGTTCCCCCGCCAgcg                                                                        >  1:71478/1‑143 (MQ=255)                                                                    gCCCGGGCCGTGCGTATCGGTGTTATTCCGCTGCCATGCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGGGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAgcgcgc                                                                     <  2:238571/143‑1 (MQ=255)                                                                    gCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAgcgcgc                                                                     <  2:239019/143‑1 (MQ=255)                                                                                cATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGGGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAGGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCg                                                         >  2:37945/1‑143 (MQ=255)                                                                                    tCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAACACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCgctct                                                     >  1:397468/1‑139 (MQ=255)                                                                                               ccGCTCCCAGTCGTCGCTGACGAGGGTGTGCCATTGCTCGAAGGAGAGACGCGTGCTGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCAt                                          <  1:525376/143‑1 (MQ=255)                                                                                                                       tttttGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCAAAACAAATTGGTGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGCCACACGCCCACCtct                  >  2:525536/1‑141 (MQ=255)                                                                                                                                       gaGGAGATACGAGTGACGGGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGt  <  2:241745/142‑1 (MQ=255)                                                                                                                                              |                                                                                                                                        TGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAACTGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGT  >  NC_000913/793440‑793716

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