Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,661,680 G→A W375* (TGG→TGA ynfF → putative selenate reductase YnfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,661,6800GA96.2% 79.3 / NA 26W375* (TGG→TGAynfFputative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/11);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCTGGGCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGGGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAG  >  NC_000913/1661543‑1661810
                                                                                                                                         |                                                                                                                                  
atTCTGGGCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTgagg                                                                                                                                  >  2:534157/1‑140 (MQ=255)
     gggCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCcac                                                                                                                             >  1:139977/1‑140 (MQ=255)
     gggCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCcac                                                                                                                             >  2:469136/1‑140 (MQ=255)
      ggCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGTGATCGGTTGAGCTACACCTGCTTATATTTGATAGGGTTGAGGGCcaca                                                                                                                            >  2:287700/1‑140 (MQ=255)
      ggCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCcaca                                                                                                                            >  2:569561/1‑140 (MQ=255)
      ggCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCcaca                                                                                                                            >  1:643889/1‑140 (MQ=255)
           aGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAACCCTGCTTATATGTGCCAGGGTGGAGGGCCACAACGAc                                                                                                                       >  1:271168/1‑140 (MQ=255)
                 tGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCa                                                                                                                 <  2:139977/140‑1 (MQ=255)
                    cGGCATAGCTAAAACGCCGGAATGGGCAGCAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACgg                                                                                                             <  1:408119/140‑1 (MQ=255)
                       cATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCg                                                                                                           <  1:8599/140‑1 (MQ=255)
                                    ccGGAATGGGCAGCATAAATCACCAGCATACCGGCAGTAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCcg                                                                                              <  1:287700/140‑1 (MQ=255)
                                    ccGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCcg                                                                                              <  2:643889/140‑1 (MQ=255)
                                      ggAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTTTCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGTCCACAACGACACTCAAACGGCCAGCAAACATCCcgcg                                                                                            >  1:292148/1‑140 (MQ=255)
                                                    aaaTCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGGGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCtt                                                                              >  2:57413/1‑140 (MQ=255)
                                                     aaTCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCttt                                                                             <  1:605364/140‑1 (MQ=255)
                                                     aaTCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCttt                                                                             <  1:112035/140‑1 (MQ=255)
                                                          ccAGCATCCCGGCAGAAAAAATCATCCAGTTTGCAGGAGCGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGtt                                                                        <  2:271168/140‑1 (MQ=255)
                                                              cATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCa                                                                    >  2:477663/1‑140 (MQ=255)
                                                                tCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGtt                                                                        >  2:339818/1‑134 (MQ=255)
                                                                tCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGtt                                                                        <  1:339818/134‑1 (MQ=255)
                                                                                            aCGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACt                                      >  1:545207/1‑140 (MQ=255)
                                                                                                             cAAAACCTGCTTATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAggg                     >  2:531266/1‑139 (MQ=255)
                                                                                                                       ttATATTTGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATc           >  2:371583/1‑140 (MQ=255)
                                                                                                                        tataTTTGTCAGGGTTGGGGGCTACAACGACATTCCAACAGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCt          <  2:243739/140‑1 (MQ=255)
                                                                                                                             ttGTCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTggggg     <  1:617170/140‑1 (MQ=255)
                                                                                                                             ttGTCAGGGTTGAGGGCCACAAAGACATTCCTACGTCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTggggg     <  2:292892/140‑1 (MQ=255)
                                                                                                                                tCAGGGTTGAGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAg  >  1:444373/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                  
ATTCTGGGCGAAGGGCCTGACGGCATAGCTAAAACGCCGGAATGGGCAGCAAAAATCACCAGCATCCCGGCAGAAAAAATTATCCAGTTGGCACGAGAGATCGGTTCAGCAAAACCTGCTTATATTTGTCAGGGTTGGGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAG  >  NC_000913/1661543‑1661810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: