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

ACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACA  >  NC_000913/1793754‑1794032                                                                                                                                                |                                                                                                                                        gcTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTcc                                                                                                                                   <  1:540920/148‑1 (MQ=255)  cTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCa                                                                                                                                  <  2:444636/149‑1 (MQ=255)  cTGCGGTCGGCGCTGCCGCAATCTATTGGTGATACATCGGATGGAGTCCATCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCa                                                                                                                                  <  2:327617/149‑1 (MQ=255)        tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACt                                                                                                                            >  1:213761/1‑149 (MQ=255)                     aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGaa                                                                                                               >  2:211763/1‑149 (MQ=255)                     aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGaa                                                                                                               >  1:621027/1‑149 (MQ=255)                     aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGaa                                                                                                               >  1:310866/1‑149 (MQ=255)                     aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTGTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGaa                                                                                                               >  2:404447/1‑149 (MQ=255)                           tttttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAg                                                                                                         <  1:634367/149‑1 (MQ=255)                                  tACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                  >  1:258268/1‑149 (MQ=255)                                      gCGGAAGGCGTCTTTCTCCATTAACTTCAATCCTACTGAACATCGGGAACGTCACCCCCCAAGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaacc                                                                                              <  1:565128/149‑3 (MQ=255)                                            ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                        >  1:531428/1‑149 (MQ=255)                                                      gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCGCGGGAATCCCGGCACCATACAACCTCGATCGACCGAGGc                                                                              >  1:143483/1‑149 (MQ=255)                                                          ttAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCCGCACCATAAAACCTCCATCGACCCCGGCTTTc                                                                          >  2:308904/1‑149 (MQ=255)                                                          ttAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTc                                                                          >  2:97703/1‑149 (MQ=255)                                                              cTTCAATCTTACTGAACATCGGTAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                                      <  2:103077/149‑1 (MQ=255)                                                               ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                     <  2:620135/149‑1 (MQ=255)                                                               ttCAATCTCACTGCACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                     <  2:247814/149‑1 (MQ=255)                                                                tCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAata                                                                    >  1:18130/1‑149 (MQ=255)                                                                     cTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                               >  2:760840/1‑149 (MQ=255)                                                                      ttACTGAACATCGGGAACGTCACCCTCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc                                                              >  2:541170/1‑149 (MQ=255)                                                                      ttACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc                                                              >  1:665708/1‑149 (MQ=255)                                                                          tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCCTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                          >  2:796664/1‑149 (MQ=255)                                                                          tGAACATCGGGAACGTCACCCCCCATGTGATGGTACAGTAAGTTTTAATCTCTTCATCGCGGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAAGCCCAGCAGCATAAAACCGCGATCGACCCAGGCTTTCGTAATATTCTCCAACt                                                          >  1:562822/1‑149 (MQ=255)                                                                             aCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTACACGGGAATCCCAGCACCATAAAACCTCGAACGACCCAGGCTTTCTGAATATCCTCCAACTGCt                                                       >  2:260837/1‑149 (MQ=255)                                                                                tCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTcgt                                                    >  1:36182/1‑147 (MQ=255)                                                                                tCGGGAACGTCACCCCCCATGTAATGGTACAGTCAGAATTAGTCACTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCAGGCTTTCTGAATATGCTCCAACGCCTCat                                                    >  2:137813/1‑149 (MQ=255)                                                                                tCGGGAACGTCACAGACCTTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                                                    >  2:492298/1‑149 (MQ=255)                                                                                    gAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTg                                                <  2:310866/149‑1 (MQ=255)                                                                                        gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑GGTTCAAGAAAATGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTacggc                                            >  2:248601/1‑149 (MQ=255)                                                                                            ccccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                                        >  1:126420/1‑149 (MQ=255)                                                                                            ccccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTCCTCATATTTCGGCGTTa                                        >  2:65907/1‑149 (MQ=255)                                                                                              ccccATGTGGTGGTACAGTAAGTTTTAATCTCTCAATAATTACCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGCATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGACTATTCTCCAACTGCTCATATTGCGGCGTTAAg                                      >  1:389087/1‑149 (MQ=255)                                                                                                 cacggggttgtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCATCTGCTCATATTGCGGCGTTAAGCca                                   <  1:88658/144‑1 (MQ=255)                                                                                                       gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGTTTC‑TGTTCAAGATCCTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGCTCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                             <  1:107850/149‑1 (MQ=255)                                                                                                          gtACAGTAAGTTTTAATCTCTTCATCGCTGCCTGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                          <  2:450144/149‑1 (MQ=255)                                                                                                              aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGCGCCGa                      >  2:719374/1‑149 (MQ=255)                                                                                                               gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATATAACCTCGATCGACCCAGGCTTTCTGCATGTTATACAACTGCTCATATTGGGGGGTTAAGCAATAGTTTGATCCGAc                     >  2:550992/1‑149 (MQ=255)                                                                                                                       tAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACaa             >  1:468099/1‑149 (MQ=255)                                                                                                                           ctctTCACCGCTGCCCGGTTC‑TGTTACAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGGCTGGCTGCACTTTGAGGCGACATTGACAATCaa         <  2:174548/149‑1 (MQ=255)                                                                                                                              ttCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg      >  1:197794/1‑149 (MQ=255)                                                                                                                              ttCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg      >  1:594487/1‑149 (MQ=255)                                                                                                                              ttCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg      >  2:276020/1‑149 (MQ=255)                                                                                                                                cATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAGCCTCGATCGACCCAGGCTTTCAGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAacaaca    >  1:114665/1‑149 (MQ=255)                                                                                                                                cATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCAGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACGTTGACAATCAACAgga    >  1:136935/1‑147 (MQ=255)                                                                                                                                  gcgcgGCCCGGTTC‑TGTTCCAGAAACGGGTTGCACGGTAATCCCAGCACCAGAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCACCTGCTCATATTGCGGCGTTCGGCCACACTTTGAGACGACATTGACAATCAACAGcaca  <  2:405773/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                        ACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACA  >  NC_000913/1793754‑1794032

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