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.92.9% 110.7 / ‑0.5 28coding (213/552 nt)btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (15/11);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGG  >  NC_000913/1793753‑1794039
                                                                                                                                                |                                                                                                                                              
gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                            <  1:487936/149‑5 (MQ=255)
   tGCGGTCTGCGCTGCCGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAGCTTCATTCTTACTGAACATCCGGACCGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAg                                                                                                                                        <  2:80249/149‑1 (MQ=255)
     cGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGaa                                                                                                                                      <  1:156225/149‑1 (MQ=255)
     cGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGaa                                                                                                                                      >  2:16576/1‑149 (MQ=255)
                                   aCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                        >  1:700311/1‑149 (MQ=255)
                                     aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                      >  1:116695/1‑149 (MQ=255)
                                              cGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcg                                                                                             <  1:385524/149‑1 (MQ=255)
                                              cGTCCTTCGCCATTAACTTCAATCTTACTGAACATCCGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcg                                                                                             <  1:435554/149‑1 (MQ=255)
                                                     cGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCTgg                                                                                      >  2:551780/1‑149 (MQ=255)
                                                      gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGc                                                                                     >  1:745133/1‑149 (MQ=255)
                                                      gCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGGTGCACGGGAATCCCAGCACCATAACACCTCGATCGACCCAGGc                                                                                     >  2:121820/1‑149 (MQ=255)
                                                        cATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCtt                                                                                   >  1:704516/1‑149 (MQ=255)
                                                              cTTCAATCTTACTGAACATCGGGAACGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                                             <  1:641420/149‑1 (MQ=255)
                                                                   aTCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                                                        <  1:708604/149‑1 (MQ=255)
                                                                                 cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTCCTCata                                                          >  2:223360/1‑149 (MQ=255)
                                                                                 cGGGAACGTCACCCCCCATGTGGAGGTAAAGTAAGGTTTAATCTCATAATCGGTCCCCACCTCTTTCCCAGCAACCGGTTGCCCCGGTTTCCCAGTACCCTACAACCTCGCTCGACCCCGGCTTTCTCAATATTCTCCAACTGCTCatc                                                          >  2:299304/1‑148 (MQ=255)
                                                                                        gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATCTTCTCCAACTGACAATATTgcggc                                                   >  2:564041/1‑149 (MQ=255)
                                                                                        gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                   >  2:269081/1‑149 (MQ=255)
                                                                                          cagaccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTCCTGTTCCCGAAACTGCTTGCACGGGAATCACAGCACCATAAAACCTCGACCGAGCCGGGCTGTCTGAATATTCTCCAACTGCCCATATTGCAGCGt                                                 >  2:299142/5‑149 (MQ=255)
                                                                                                ccATGTGGTGGTACAGTAAGTTTTAATCTCTGCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGcc                                           <  2:306349/149‑1 (MQ=255)
                                                                                                           tACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAgg                                <  1:16576/149‑1 (MQ=255)
                                                                                                            aCAGTAAGTTTTACTCTCCTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGc                               <  2:258392/149‑1 (MQ=255)
                                                                                                                        aaTCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACCAt                   <  2:721026/149‑1 (MQ=255)
                                                                                                                           ctctTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCGTATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACacacaa                >  2:149263/1‑144 (MQ=255)
                                                                                                                                      tGCCCGGTTCTGTTCCAGAAACTGGTTGCGCGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAACCAACAGCACATTAc     >  1:491181/1‑149 (MQ=255)
                                                                                                                                       gCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAcc    <  1:299142/149‑1 (MQ=255)
                                                                                                                                        cccGGTTCTGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATCGACAATCCACAGCACATAACCg   >  2:581732/1‑149 (MQ=255)
                                                                                                                                         ccGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCgg  >  1:625044/1‑149 (MQ=255)
                                                                                                                                         ccGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCgg  >  2:643436/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGG  >  NC_000913/1793753‑1794039

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: