Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,281,111 G→A W160* (TGG→TGA)  radD → putative DNA repair helicase RadD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,281,1110GA96.5% 91.1 / ‑3.8 29W160* (TGG→TGA) radDputative DNA repair helicase RadD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (14/14);  total (15/14)
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.15e-01

TCACCGTATTGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGA  >  NC_000913/2280988‑2281239
                                                                                                                           |                                                                                                                                
tCACCGTATTGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCAtt                                                                                                                  >  2:82773/1‑140 (MQ=255)
    cGTATTGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCa                                                                                                              >  2:214146/1‑140 (MQ=255)
      tATTGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACg                                                                                                            >  2:6069/1‑140 (MQ=255)
        ttGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGc                                                                                                          <  1:214146/140‑1 (MQ=255)
         tGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCCGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCa                                                                                                         <  1:366728/140‑1 (MQ=255)
           gTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATg                                                                                                       >  1:98064/1‑140 (MQ=255)
                     agagagCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACgcggcg                                                                                             <  2:146367/140‑1 (MQ=255)
                           ccAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTAcg                                                                                                  <  2:39939/129‑1 (MQ=255)
                           ccAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTAcg                                                                                                  >  1:39939/1‑129 (MQ=255)
                           ccAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATgaga                                                                                       <  2:98064/140‑1 (MQ=255)
                           ccAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATgaga                                                                                       >  1:73525/1‑140 (MQ=255)
                           ccAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATgaga                                                                                       <  1:82773/140‑1 (MQ=255)
                                        aTCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCg                                                                          >  2:184498/1‑140 (MQ=255)
                                         tCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGt                                                                         >  2:98237/1‑140 (MQ=255)
                                           ctcactcaCCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGa                                                                       >  2:398818/1‑140 (MQ=255)
                                                         aaaaGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGc                                                         >  2:75163/1‑140 (MQ=255)
                                                                 aTCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTg                                                 <  2:73525/140‑1 (MQ=255)
                                                                 aTCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTg                                                 <  1:98237/140‑1 (MQ=255)
                                                                     ccACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGcc                                                     >  2:105973/1‑132 (MQ=255)
                                                                     ccACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGcc                                                     <  1:105973/132‑1 (MQ=255)
                                                                      cACTTACGGCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTta                                            <  2:254816/140‑1 (MQ=255)
                                                                           aCGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGa                                       >  1:142842/1‑140 (MQ=255)
                                                                             gCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTg                                                 >  2:346653/1‑128 (MQ=255)
                                                                             gCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTg                                                 <  1:346653/128‑1 (MQ=255)
                                                                                                         tCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcgcc                 <  2:163187/132‑1 (MQ=255)
                                                                                                         tCGATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGAcgcc                 >  1:163187/1‑132 (MQ=255)
                                                                                                           gATTGGGCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGa       <  1:6069/140‑1 (MQ=255)
                                                                                                                ggCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGa  <  2:142842/140‑1 (MQ=255)
                                                                                                                ggCAAAGGCTGAATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGa  >  1:40191/1‑140 (MQ=255)
                                                                                                                           |                                                                                                                                
TCACCGTATTGGTGACGATGAAGAGAGCCAGTATCAGCAAATCCTCACTCACCTGACAAAAGTGAATCCCCACTTACGCCTGCTGGGGCTGACTGCCACGCCTTTTCGATTGGGCAAAGGCTGGATCTACCAGTTTCATTATCACGGCATGGTACGCGGCGATGAGAAAGCCCTTTTCCGTGACTGCATTTATGAGCTGCCGCTGCGTTATATGATTAAACACGGCTATCTGACGCCGCCAGAACGACTGGA  >  NC_000913/2280988‑2281239

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