Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,793,897	Δ1 bp	100%	coding (213/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	.	100.0%	137.3 / NA	31	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (15/16);  total (15/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAC  >  NC_000913/1793753‑1794036                                                                                                                                                 |                                                                                                                                            gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGGCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgctc                                                                                                                                         >  2:376948/1‑145 (MQ=255)                     caatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGa                                                                                                                    <  1:70944/149‑1 (MQ=255)                         caatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCCTACTGAACGTCGGGAACGGCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATcc                                                                                                                <  2:38331/149‑1 (MQ=255)                           atTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGGATCCCa                                                                                                              >  1:736510/1‑149 (MQ=255)                           atTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCa                                                                                                              <  2:60232/149‑1 (MQ=255)                           atTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCa                                                                                                              >  1:480127/1‑149 (MQ=255)                           atTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTCGTTGTCCGGGAATCCCa                                                                                                              >  2:195354/1‑149 (MQ=255)                                    aCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                     >  1:98829/1‑149 (MQ=255)                                         ggATGGCGTCCTTCGCCATTAACTTCCATCTTACTGAACATCGGGGGCAACGTACCGGATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCt                                                                                                <  1:24834/149‑1 (MQ=255)                                             ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                            <  2:733129/149‑1 (MQ=255)                                                gTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcga                                                                                         <  2:316491/149‑1 (MQ=255)                                                    ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                     <  2:536996/149‑1 (MQ=255)                                                       gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCCGCACCATAAAACCTCGATCGACCCAGGc                                                                                  >  2:548048/1‑149 (MQ=255)                                                          aTTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTCATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCttt                                                                               >  1:573191/1‑149 (MQ=255)                                                            tAACCTCAATCTTCCTGAACAGCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                             <  2:19736/149‑1 (MQ=255)                                                              aCTTCCATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGCCTTTCTGa                                                                           <  2:612708/149‑1 (MQ=255)                                                               tctCAATCCTACTGCACATCGTGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCCACATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                                          <  1:201758/147‑1 (MQ=255)                                                                ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                         >  2:791756/1‑149 (MQ=255)                                                                           tGAACATCGGGAACGTCACCCCCCATGTGGTGGTCCAGGAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGCTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCAATCGACCCAGGCTTTATGAAAATTCACCAACt                                                              >  2:433112/1‑149 (MQ=255)                                                                                   gggAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGGATCCCAGCACCATAAAACCTGGATCGACCCAGGCTTTCGGAATATTCTCCAACTGCTCctat                                                      >  2:17538/1‑149 (MQ=255)                                                                                       aCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                  <  1:395407/149‑1 (MQ=255)                                                                                                     gtggtgGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGccgccc                                    >  1:107806/1‑145 (MQ=255)                                                                                                     gtggtgGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCacac                                    >  1:544152/1‑149 (MQ=255)                                                                                                         tggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCCACTGCTCATATTCCGGCCTTACGCCACACTTTg                                >  2:161837/1‑149 (MQ=255)                                                                                                                 tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAGTATTCTCCAACTGCTCATATTGCGGCGTTAAGACACACTTTGAGGCGCCa                        >  2:56950/1‑149 (MQ=255)                                                                                                                  aaGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGCAAATGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACAATTTGAGGCGACAt                       <  1:583211/149‑1 (MQ=255)                                                                                                                   aGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTACACGGGAATTCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATCTTGCGGCGATAAGCCACACATTGAGGCGACAtt                      >  1:381115/1‑149 (MQ=255)                                                                                                                    gTTTTAATCTCTTCCTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTg                     <  2:765504/149‑1 (MQ=255)                                                                                                                     ttttAATTCCTTCATCGCGGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                    <  1:29510/149‑1 (MQ=255)                                                                                                                       ttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGCCGACATTGACa                  >  2:53586/1‑149 (MQ=255)                                                                                                                               ttCATCGCTGCCCGGTTC‑TGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg          <  2:734879/149‑1 (MQ=255)                                                                                                                                       tGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATACTACTACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAc  <  1:376948/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAC  >  NC_000913/1793753‑1794036

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