Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,760,383 A→C 29.8% V289G (GTG→GGG)  selA ← selenocysteine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,760,3830AC29.8% 33.2 / 2.2 29V289G (GTG→GGG) selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base C (2/7);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCT  >  NC_000913/3760146‑3760611
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCAATCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGCGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTATTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACCTCTCGCCGGGTAAAATCa                                                                                                                                                                                                                                      >  1:493258/1‑238 (MQ=255)
 ccTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCAc                                                                                                                                                                                                                                     <  2:57639/238‑1 (MQ=255)
    ggATTTTAATTACTTCTGGCGTGTGGTTAAGCAGGGGAAGGGTCGGTAATTTTTCACTCAGATCTTCAGGGTGTAAATAAAGAGCCAGCGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCATGCTGCAGCGGGTGGTTTGGCAGGCGGGCGATCATTTCTTTTTTACCAACAATAATTCCTGCCGGCGGCCCGCCTAACAGCTTGTCGCCGGAGAAACTCACCAGa                                                                                                                                                                                                                                 <  2:342541/239‑1 (MQ=255)
    ggATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCCAACAATAATTCCTTCCTGCTGCCCGCCTAACCACTTGTCGCCGGAGAAACTCCCCaa                                                                                                                                                                                                                                  >  2:499748/1‑237 (MQ=255)
           aaTGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCCGCCGGCCGCCCAACAAATTGTGGCCGGAGAAACCCACCAGACTGAcg                                                                                                                                                                                                                           >  2:390356/1‑238 (MQ=255)
                         tgcgGGTAAGCAGGTCAAGGGCCGGTATTTTTTCACTCAGAGCTACAGGGTGTAAATAAAGACGCAACGGGGCTCCCAGCGCCGCGGGGGTCATTTTACCCGCGTGTAATGCACGCTCCAGGGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAAAAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGCCGCCCGCCGCAATCaa                                                                                                                                                                                                             <  2:269345/238‑1 (MQ=255)
                         tgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAATTCCTTCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGAATGACGCCCGCCGCAATCaa                                                                                                                                                                                                             >  2:88736/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTACCAACTTGTCCCCGGAGAAACTCACCCGACTGCCCCCCGCCCCAATCAACt                                                                                                                                                                                                           >  2:93892/1‑238 (MQ=255)
                                                 gTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATTTCTTTTTTATCAACAATTAGTCCTCCCTGCGGCCCCCCTAACAACTTGTCGCCGGGGAAACTCACCAGACTGACGCCCGCCGGCATCAACTCCTCCCGCAGTGACCCtttacg                                                                                                                                                                                     >  2:152097/1‑235 (MQ=255)
                                                         tccgcaggagCTGCGGGGTAAAAAGAAAGACGCAACGTGGCTTACGGCGCCGGGGGGGTTATTTTATCCGCGTGTAAGCCACGCTTCAGGGGGTGGCTGTGCGGGGGGGCGATCATCTCTTTTTTACAAAAAATAATTCCTGCCTGCGCCCCGCCAAACACCTTGTCGCGGGAAAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:365581/232‑1 (MQ=255)
                                                         tcccgcagtcctgcggggTGTAAAAAGAGACGCAGCGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTCCAGGGGGGGGCTTGGCAGGCGGGCGATCATCTCTTTTTTCCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTGGTCGCCGGAAAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:500506/224‑1 (MQ=255)
                                                                   cTTCGGGGTAAAAAAAAAGCCGCAACGTGGCTCCCAGCGCCGCGAGGGTCTTTTTATCCGGGCGAAATGCCCGCTTCGGGGGGTGGCTTGGCGGGGGGGCGATCATCTCTTTTTTACCAACAATAATCCCTGCCGGCGCCCCGCCTAACAACTTGTCGCGGGAGAACCCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAg                                                                                                                                                                   <  2:349829/238‑1 (MQ=255)
                                                                                  aacggcGCAGCGTGGCTCCCGGCGCCGGGGGGTCTATTTTATCGGCGTTTAGTGCCCGCTTCAGGGGGTGTTTTTGCGGGCGGGCGTTCATCTTTTTTTTCCCAAAAAAAATCCCTCCCTGCGGCCCGCCAAACATTTTGTCGCCGGAAAACCCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                   <  2:304450/234‑1 (MQ=255)
                                                                                        gcacgggggCCTCCAGCGCCGAGAGGTGCATTTGACCCCCGTGTAACGCACGCTCCGGGGGGTGGCTTGGCAGGGGGGCACCCATCTTTTTTTTCCCAACAATAATCCCCCCCGGCGGCCCGCCTAAAACCTTGCCCCGGGAAAAACCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCt                                                                                                                                             <  2:493258/232‑1 (MQ=255)
                                                                                            cGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCCACAATAATTCCTTCCTGCGGCCCGCCCAACAAATTTTCGCCCGAGAAAATCCCCCGACTGACACCCGCCGCCATCAACTCCCGCTGCATTTGCTCGTTCGGACAAACGCAACGGCCAAAATCGACCAGccggccccccgccaaat                                                                                                                                         >  2:436655/1‑223 (MQ=255)
                                                                                                                                cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCCCCGCCTAAAACAATCACTACGGGAACCTCCCGCTCTTTGCCGAg                                                                                                       >  2:381849/1‑238 (MQ=255)
                                                                                                                                                        ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCCCTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCCAGCGCCCACCGTTCCGCCTCATCTAt                                                                               >  2:230588/1‑238 (MQ=255)
                                                                                                                                                                    gggcgATCATCTCTTTTTCACAAAAAATAATCCCGGCCTGCGGCCCGCCTAAAAACTTGTCGCCGGAAAACCTCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                   <  2:449297/238‑1 (MQ=255)
                                                                                                                                                                            atcTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAAt                                                                                                                                         >  1:166308/1‑160 (MQ=255)
                                                                                                                                                                             tttctTTTTTCCCAAAAAAACTCCCTGCCGGCGCCCCGCCAAACAACTTGCCGCCGGAGAACCCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATg                                                         <  2:165259/237‑1 (MQ=255)
                                                                                                                                                                               tctTTTTTACAAACAAAAATCCCGGCCGGCGCCCCGCCTAAAAATTTGTCGCGGAAGAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                       <  1:285478/239‑1 (MQ=255)
                                                                                                                                                                                ctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGact                                                          >  1:165259/1‑233 (MQ=255)
                                                                                                                                                                                    tttACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGGCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                                  >  1:456200/1‑239 (MQ=255)
                                                                                                                                                                                      caCCAACAAAAATCCCCGCCTGCGCCCCCCTTACCAACTGGCCGCGGGAGAACCTCACCAGACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  1:381849/238‑1 (MQ=255)
                                                                                                                                                                                                acttcccgccgcccgCCCCCCCAACACGTTGTCGCCGGAAAACCACACCAGTCGGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  1:244424/225‑1 (MQ=255)
                                                                                                                                                                                                aatCCCTGCCTGGGCCCACCTTAACAATTGGTCGCCGGAGAAACTCCCGACACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  1:34006/239‑1 (MQ=255)
                                                                                                                                                                                                aatCCCTGCCGGCGGCCCGCCAAACACCTTGCCGCCGAAAAACCGCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  1:502259/239‑1 (MQ=255)
                                                                                                                                                                                                aatACCTGCCGGCGCCCCCCCTAACACCTGGCCCCAGGGAAAACGCGCGAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  2:8078/239‑1 (MQ=255)
                                                                                                                                                                                                       gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccacca                                                                                               >  1:45874/1‑175 (MQ=255)
                                                                                                                                                                                                                                      gaaaaaCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCt  <  1:499748/235‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCT  >  NC_000913/3760146‑3760611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: