Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,898	Δ1 bp	coding (212/552 nt)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,897	0	T	.	95.0%	80.1 / ‑2.9	20	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (14/5);  minor base A (1/0);  total (15/5)
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

GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAC  >  NC_000913/1793753‑1794026                                                                                                                                                 |                                                                                                                                   gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                <  1:376568/149‑5 (MQ=255)             cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑Gtt                                                                                                                  >  1:117228/1‑149 (MQ=255)                                      cgCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATaaaa                                                                                         >  2:327158/2‑149 (MQ=255)                                      aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATaaaa                                                                                         <  2:5100/149‑1 (MQ=255)                                      aGCGGATGGCGCCCTTCGCCATCAACTTCAAGCTTACAGGACCTCGCGAACGTCACCCCCCATGGGGTCGTACGATAAGTTTTAATCTCTTCATCGCTGCCCCGTT‑ATGGTCCAGAGACTGCGTTG‑ACGTGAATGCCAGCACCATaata                                                                                         >  2:241177/1‑147 (MQ=255)                                                       gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGACACTG‑GTTGCACGGGAATCCCAGCACCATAAAACATCGTTCGACCCACGc                                                                        >  2:268229/1‑149 (MQ=255)                                                       gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGc                                                                        >  1:339575/1‑149 (MQ=255)                                                          aTTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCttg                                                                     >  1:179999/1‑148 (MQ=255)                                                               cTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCTACCCAGGCTTTCTGaa                                                                <  2:442211/149‑1 (MQ=255)                                                               cTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                                >  1:406890/1‑149 (MQ=255)                                                                 tCAAACTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAACCTCTTCATCGCCGCCCGGTTC‑TGCTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAata                                                              <  2:55306/149‑1 (MQ=255)                                                                   aaTCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATAtt                                                            <  2:301945/149‑1 (MQ=255)                                                                      cTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                         <  1:177266/149‑1 (MQ=255)                                                                             aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                                                  >  1:133889/1‑149 (MQ=255)                                                                                                    tgtgGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAACACCTCGATCGCCCCGGGCTTTCTACATGTTATACCACGGCTCACATTGCGGCGGTAAGCcacc                           >  1:170482/1‑148 (MQ=255)                                                                                                        gtggtACAGTAAGTTTTAATCTCTTCATCGATGCCCGGTTC‑TGTTCCAGAAACTG‑GTTACACGGGAATCTCAGCACCATAAAACCTCGATCGACCCAGGATTTATGAATATTCTTCAACTGGTCATATTGGGGCGTTTAGcccaaattt                       >  2:269180/1‑142 (MQ=255)                                                                                                                     ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGCGGCGACATTGa          >  1:124424/1‑149 (MQ=255)                                                                                                                     ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa          >  1:427714/1‑149 (MQ=255)                                                                                                                     ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGATTTCTGACTATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTGTGAGGCGAAATTgc          >  2:125081/1‑148 (MQ=255)                                                                                                                     ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGAGTATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCGCACTTTGAGGCGCCATTgc          >  1:349490/1‑148 (MQ=255)                                                                                                                             tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc  >  1:351300/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                   GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTG‑GTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAC  >  NC_000913/1793753‑1794026

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