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%	232.1 / NA	50	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (27/23);  total (27/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACC  >  NC_000913/1793750‑1794037                                                                                                                                                    |                                                                                                                                             cgcgACTGCGGTCGGCGCTGCGGAAAGCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGCCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtg                                                                                                                                             <  2:815816/149‑2 (MQ=255)           gTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAAc                                                                                                                                  <  2:737953/149‑1 (MQ=255)             cGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTg                                                                                                                                <  1:928255/149‑1 (MQ=255)               gcgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGt                                                                                                                              >  2:445324/1‑149 (MQ=255)                  cTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTCCTGAACATCGGGAACGTCACCCCCCATGTCGGGGTACAGTAAGTATTAATCTCTTCATCGCTGCCCGCTTC‑TGTTCCAGAAACTGGTTGc                                                                                                                           >  2:127081/1‑149 (MQ=255)                           tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                  >  2:755971/1‑149 (MQ=255)                           tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGGCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTGC‑TGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                  >  1:247001/1‑149 (MQ=255)                             aatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGCTCCAGAAACTGGTTGCACGGGAATacc                                                                                                                >  2:521549/1‑149 (MQ=255)                                      tACAGCGGATGGCGTCCTTCCCCATTAACTTCACTCTTACTGAACATCGGGAGCGTCCCCCCCCGTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                       <  2:394684/149‑1 (MQ=255)                                      tACAGCGGATGGCCTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                       <  2:1067816/149‑1 (MQ=255)                                                ggcggcCTTCGCCACTATCTTCAATCTCACCTAACATCGGGACCGTTCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTCCATCGCTGCCCGGTTC‑TGTTCCAGAAACTTGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                             <  2:264469/149‑1 (MQ=255)                                                ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                             <  1:143442/149‑1 (MQ=255)                                                cgCGTCCGTCTCGATTAACTTCAATCTTACTGAACATCGGGAACGTCCCCCCCCATGTGTTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGCGAATCCCAGCACCATAAAACCtcgat                                                                                             <  2:1168372/148‑1 (MQ=255)                                                   gTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcga                                                                                          <  2:864545/149‑1 (MQ=255)                                                    tCCTGCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCCCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAc                                                                                         <  1:1039370/149‑1 (MQ=255)                                                       ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGCAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGACCGCCCCa                                                                                      >  2:1179392/1‑149 (MQ=255)                                                       ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                      >  2:112630/1‑149 (MQ=255)                                                       ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCACGTGGTGGTACAGTAAGTTTTAATCTCATCATCGCTCCCCGGCTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATACACCCa                                                                                      >  2:851433/1‑149 (MQ=255)                                                                   ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                          <  1:1120279/149‑1 (MQ=255)                                                                   ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCACAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                          <  2:1248072/149‑1 (MQ=255)                                                                         cTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCGCTTCATCGCTGCCCGGTTC‑TGGTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTGCTGAATATtctc                                                                    <  2:219664/149‑1 (MQ=255)                                                                                     cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCCCata                                                        >  2:179075/1‑149 (MQ=255)                                                                                     cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCATCAAAGCTGCTGGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCCGGCTTTCTGAATATTCTCCAACTGCCCata                                                        >  1:256140/1‑149 (MQ=255)                                                                                     cGGGAACGTCACCCCCCATGTGGTGGTACAATAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCata                                                        >  1:1048087/1‑149 (MQ=255)                                                                                          aCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                   >  1:1203165/1‑149 (MQ=255)                                                                                            gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                 <  1:445324/149‑1 (MQ=255)                                                                                            gTCACCCCCCATGTGGTGGTACAGAAAGCTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                 >  1:968668/1‑149 (MQ=255)                                                                                            gTCACCCCCCATGGGGTGGTACAGTAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                 <  2:1288007/149‑1 (MQ=255)                                                                                                 cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGTCGTTaa                                            >  1:39937/1‑149 (MQ=255)                                                                                                 cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                            <  2:897871/149‑1 (MQ=255)                                                                                                  ccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCCGCACCATAAATCCTCGATCGTCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                           <  1:755971/149‑1 (MQ=255)                                                                                                   cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCACCTGCTCATATTGCGGCGTTAAGc                                          <  2:1128973/149‑1 (MQ=255)                                                                                                    ccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATCAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGCCGTTATGcc                                         >  2:706408/1‑149 (MQ=255)                                                                                                    ccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCCCATATTGCGGCGTTAAGcc                                         >  2:1053175/1‑149 (MQ=255)                                                                                                     cATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                                        <  1:477161/149‑1 (MQ=255)                                                                                                         tggtggTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATCTTGCGGCGCTAAGCCCCACt                                    >  2:420082/1‑149 (MQ=255)                                                                                                             ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                <  1:970923/149‑1 (MQ=255)                                                                                                                   gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGAc                          <  2:751662/149‑1 (MQ=255)                                                                                                                       gTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGCCGACATTg                      >  2:192351/1‑149 (MQ=255)                                                                                                                         tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAc                    >  1:166996/1‑149 (MQ=255)                                                                                                                         tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAc                    >  1:988588/1‑149 (MQ=255)                                                                                                                           tAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACaa                  >  2:1206874/1‑149 (MQ=255)                                                                                                                            aaTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAAt                 <  2:817454/149‑1 (MQ=255)                                                                                                                             aTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATc                >  2:414324/1‑149 (MQ=255)                                                                                                                                tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc             >  1:440166/1‑149 (MQ=255)                                                                                                                                    cATCGCTCCCGCGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCCGCACCATAAAACCTCGATCGCCCCAGGCTTTCCGAATCTTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca         <  2:90343/149‑1 (MQ=255)                                                                                                                                    cATCCCTCCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca         <  2:90455/149‑1 (MQ=255)                                                                                                                                        gCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGACGTTAAGCCACCCTTTGAGGCGACACTGACAATCACCAGCCCAtt     >  2:20185/1‑149 (MQ=255)                                                                                                                                          tGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAc   >  1:474823/1‑149 (MQ=255)                                                                                                                                          tGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCCGCATCAACAATCAACAGCACATTAc   >  2:61734/1‑149 (MQ=255)                                                                                                                                           gCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACGCTTTGAGGCGACATTGACAATCAACCGCACATTAcc  >  1:520743/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                             CGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACC  >  NC_000913/1793750‑1794037

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: ‑