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.96.9% 137.4 / ‑2.9 32coding (213/552 nt)btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (16/15);  total (17/15)
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.98e-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.

ACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793754‑1794041
                                                                                                                                               |                                                                                                                                                
aCTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTcc                                                                                                                                            <  2:688401/149‑1 (MQ=255)
aCTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTcc                                                                                                                                            <  2:111561/149‑1 (MQ=255)
    cGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGaa                                                                                                                                        >  2:590900/1‑149 (MQ=255)
    cGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTGCTGTTCCAGaa                                                                                                                                        >  2:177371/1‑149 (MQ=255)
               gCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACCTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCAc                                                                                                                             >  1:736088/1‑149 (MQ=255)
               gCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCAc                                                                                                                             <  1:687577/149‑1 (MQ=255)
                cGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGATTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACg                                                                                                                            >  1:398852/1‑149 (MQ=255)
                 ggCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTTTCCAGAAACTCGTTGCACgg                                                                                                                           >  1:327424/1‑149 (MQ=255)
                      tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                      >  2:516510/1‑149 (MQ=255)
                                  aCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                          >  2:632730/1‑149 (MQ=255)
                                        gATGGCGTCCTTCGCCATTACCTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCATAGGCTGGTTGCACGGTAATCCCAGCACCATAAAACCtc                                                                                                    <  2:383329/149‑1 (MQ=255)
                                            gCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatc                                                                                                <  1:803718/149‑1 (MQ=255)
                                                  ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGACCGACCCa                                                                                          >  2:284569/1‑149 (MQ=255)
                                                                           aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                                                                 >  2:275214/1‑149 (MQ=255)
                                                                           aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACCGc                                                                 >  2:285101/1‑149 (MQ=255)
                                                                            aCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAACCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                                                <  1:627430/149‑1 (MQ=255)
                                                                                cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCata                                                            >  1:646865/1‑149 (MQ=255)
                                                                                 gggAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCatat                                                           <  1:375455/149‑1 (MQ=255)
                                                                                    aaCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgc                                                        <  1:565540/149‑1 (MQ=255)
                                                                                            cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                                <  2:776392/149‑1 (MQ=255)
                                                                                             ccccATTGGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATACAACCTCGAGCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                               <  2:379091/149‑1 (MQ=255)
                                                                                             ccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGGCCCAGGCTTTCTGAATATTCTCCAACTCCTCATATTGCGGCGTTAAg                                               >  2:298982/1‑149 (MQ=255)
                                                                                              cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                              >  1:866795/1‑149 (MQ=255)
                                                                                                  tgtgGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCcaca                                          >  2:615746/1‑149 (MQ=255)
                                                                                                          tACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTCCCCGATACTGGCTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCATGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAgg                                  <  1:800129/149‑1 (MQ=255)
                                                                                                             aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGa                               <  1:516510/149‑1 (MQ=255)
                                                                                                              gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCACGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGGTAGGCCACACTGTGAgacgac                              >  2:689657/1‑149 (MQ=255)
                                                                                                               aaaGTTCTAACCTCTTCATCGCTGCCCGGTTCTGGTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACa                             <  2:454297/148‑1 (MQ=255)
                                                                                                                  gTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTAAGGCGCCAtcg                          >  2:320810/1‑147 (MQ=255)
                                                                                                                     ttAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCCTATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACa                       <  1:877318/149‑1 (MQ=255)
                                                                                                                            ctTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCCTAAAACCTCGATCGCCCCAGGCTTTCTGAATATTCGCCAACGGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa                <  1:632730/149‑1 (MQ=255)
                                                                                                                                          ggTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCGGCACCATAAAACCTCGATCGACCCAGGCTATCTGCAGATTCTCCAACTGCTCATATTGCGGCATTGAGACTCACATGAAGGCGACATTGACAATCAACAGCACCTTACCggag  >  2:456975/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793754‑1794041

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: