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_0009131,793,8970T.97.6% 179.4 / ‑2.0 41coding (213/552 nt)btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (16/24);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793749‑1794040
                                                                                                                                                    |                                                                                                                                               
gcgcGTCTGCGGTCGGCGTTGCGGGAATCACGTTATCATAAAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCt                                                                                                                                                 <  1:564877/149‑1 (MQ=255)
gcgcGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCt                                                                                                                                                 >  2:644693/1‑149 (MQ=255)
                       gCAATCAATTTTTGATACTGCGGAGGGCGTCCTGCGCCATTAACTTTAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACggg                                                                                                                         <  2:47646/149‑1 (MQ=255)
                        caatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGa                                                                                                                        >  2:131554/1‑149 (MQ=255)
                           tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                     >  1:726326/1‑149 (MQ=255)
                            caatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATcc                                                                                                                    >  2:230989/1‑149 (MQ=255)
                              atTTTTGATACAGCGGATGGTCTCCTTCGCCAGGCACTGCAATCTTACTGCACATCGGGATCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCGCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCa                                                                                                                  <  2:353385/149‑1 (MQ=255)
                               tttttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAg                                                                                                                 >  1:667718/1‑149 (MQ=255)
                               tttttGACACAGCGGATGGCGTCCTTCGCCTTTAACCTCAATCTTACCCAACATCGGGACCGTCACCCCCCCTGTGGTGGTACAGTAAGTTTTAACCTCTTCATCGCTTCCCGGTTCTT‑TTCCAGAAACTGGTTCCACGGGAATCCCat                                                                                                                 <  2:197962/149‑2 (MQ=255)
                                     aTACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGGTGCACGGGAATCCCCGCACCAt                                                                                                           >  2:286377/1‑149 (MQ=255)
                                       aCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                         >  2:330125/1‑149 (MQ=255)
                                         ccccgATGGCGTCCTTCGCCATTAACTTCAATCTTACTCAACATCGGGACCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                       <  2:663677/145‑1 (MQ=255)
                                         aGCGTATTGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                       <  2:546089/149‑1 (MQ=255)
                                         aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGTGAATCCCAGCACCATaaaa                                                                                                       <  1:531681/149‑1 (MQ=255)
                                             ggtGGCGGCCTTCGCCAGTAAGTTTCACCTTGCTGAACAGCGGGCCAGTCACCTCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                                                                                   <  1:714346/147‑1 (MQ=255)
                                               tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTAATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                                                                                 <  1:387960/149‑1 (MQ=255)
                                                ggCTTCCTTCGCCAGTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGCTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                                <  1:188202/149‑1 (MQ=255)
                                                    tCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAc                                                                                            <  2:719324/149‑1 (MQ=255)
                                                                 aCTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGa                                                                               >  2:453731/1‑149 (MQ=255)
                                                                 aCTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGa                                                                               <  1:472318/149‑1 (MQ=255)
                                                                     cAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGAGGTACAGTAAGTTCTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCGGCACCATAAAACCTCGATCGACCCAGGCTTTCGGAatat                                                                           >  2:447098/1‑149 (MQ=255)
                                                                         cctACTAAACAACCGGAACCTCACCCCCCCAGTGGTGGTACAGTAAGTTTTAATCCCTTCATCGCCGCCCGGTTC‑TGTCCCACAAACTGGTTGCACGGGAATCCCAGCACCACAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                                       <  2:292948/147‑1 (MQ=255)
                                                                          ttACTGAACATCGCGACCTTCACCCCCCCTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTACAGAAACTGGTTGCACGTTAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc                                                                      <  2:628664/149‑1 (MQ=255)
                                                                               gAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTg                                                                 <  1:662868/149‑1 (MQ=255)
                                                                                     cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCCCCAACTGCTCata                                                           >  2:697541/1‑149 (MQ=255)
                                                                                               aCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCNGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGtt                                                 <  1:730002/149‑1 (MQ=255)
                                                                                                 cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATACTCCAACTGCTCATATTGCGGCGTTaa                                               >  1:511173/1‑149 (MQ=255)
                                                                                                   cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                             <  2:205708/149‑1 (MQ=255)
                                                                                                          ggtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACtt                                      >  2:371927/1‑149 (MQ=255)
                                                                                                           gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                                     <  2:60417/149‑1 (MQ=255)
                                                                                                            tggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTg                                    >  1:492167/1‑149 (MQ=255)
                                                                                                             ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                   <  1:541502/149‑1 (MQ=255)
                                                                                                             ggtACAGCAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTCTGa                                   <  2:487300/149‑1 (MQ=255)
                                                                                                              gtACAGTAAGTCTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTACAGAAACTGATTACACGGGAATCCAAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCAACACTTTAAg                                  >  2:299029/1‑149 (MQ=255)
                                                                                                                  aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                              <  1:93027/149‑1 (MQ=255)
                                                                                                                          ttAAGCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACa                      <  2:649113/149‑1 (MQ=255)
                                                                                                                          atAAGCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACa                      <  2:78680/148‑1 (MQ=255)
                                                                                                                               ctctTCATCGCTGCCCGGTTC‑TGTGCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                 <  1:644693/149‑1 (MQ=255)
                                                                                                                                 ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               <  2:725480/149‑1 (MQ=255)
                                                                                                                                  ttCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCGCCATAAAACCTCGGTCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAAAg              >  2:458769/1‑149 (MQ=255)
                                                                                                                                     aTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGcac           >  2:690682/1‑149 (MQ=255)
                                                                                                                                          tGCCCGGTTC‑TGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAc      <  1:667170/149‑1 (MQ=255)
                                                                                                                                              cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc  >  2:433166/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793749‑1794040

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