Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,791,067 A→G C93C (TGT→TGC tdh ← threonine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,791,0670AG95.0% 57.2 / ‑5.7 20C93C (TGT→TGCtdhthreonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (11/8);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTACCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGCCC  >  NC_000913/3790826‑3791313
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
gaccTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGat                                                                                                                                                                                                                                                >  1:399815/3‑250 (MQ=255)
      cACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCt                                                                                                                                                                                                                                         >  2:10120/1‑251 (MQ=255)
                     caTCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGCAGGGATCAAAAATGGACGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGAcgcg                                                                                                                                                                                                                           >  2:162685/1‑250 (MQ=255)
                        cctCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCg                                                                                                                                                                                                                       >  1:44232/3‑251 (MQ=255)
                              caaccaGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATc                                                                                                                                                                                                                 >  1:188864/3‑251 (MQ=255)
                                             aCAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACt                                                                                                                                                                                                  >  2:126116/1‑251 (MQ=255)
                                                       aTGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGCTCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGAccg                                                                                                                                                                                         >  2:161221/1‑249 (MQ=255)
                                                                                                                                                                     cgggacggtgaacgccggcGCCGATCGTGTTGCGCCAGAGAGGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCAGGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGc                                                                           <  2:47798/232‑1 (MQ=255)
                                                                                                                                                                          gggttaacgccgacgcggaTCGGGTTGCGGCACAGAGGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCAGCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGtttt                                                                     <  1:162589/233‑1 (MQ=255)
                                                                                                                                                                            gtcaccgccaccgccgacggtgtGGCGGCACAGATGGTTACGCCCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTAc                                                                   <  2:332659/232‑1 (MQ=255)
                                                                                                                                                                               aacgccgaCGCCGATCGTGTTGGGGCAAAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGa                                                                <  1:318795/244‑1 (MQ=255)
                                                                                                                                                                                acgccgaccccgtggtgtttGCGGAACATATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGAt                                                               <  1:138650/233‑1 (MQ=255)
                                                                                                                                                                                        gccgaTCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAAtt                                                        >  1:351726/5‑249 (MQ=255)
                                                                                                                                                                                          cgaTCGTGTTGCGGCACAGATGGGTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATtagc                                                     >  1:7300/3‑251 (MQ=255)
                                                                                                                                                                                                  tcgcgccacaggtgcatcctgccgccacGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTTCTACACCGACATATTCAGGACCGTCAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATTAGCAGCTCGt                                              <  2:436395/226‑1 (MQ=255)
                                                                                                                                                                                                                ggTACGACCGCCACGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGACTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATTAGCAGATCGTTATGACCCAGTTCCg                               >  1:190518/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggtccgccacGACAGTTGCGACAGGGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATTAGCAGATCGTTATGACCCAGTTCCggctc                           <  2:290646/245‑4 (MQ=255)
                                                                                                                                                                                                                       ccgcctcgaaaGTTGCGACAGTGACCACAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCCCTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATTAGCAGATCGTTATGACCCAGTTCCggctcagg                        <  2:349989/241‑7 (MQ=255)
                                                                                                                                                                                                                             cGACAGTTGCGACAGTGACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATTAGCAGATCGTTATGACCCAGTTCCGGCTCAGGAACGTc                  >  1:406818/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gACCGCAGGTGATATGGCCTTCTCCAGAGACGCGATCGCCAATCTTAAAGCCTTTCACTTCCTGACCGATGCCTACTACTTCACCGACATATTCATGACCGACAACCATCGGCACCGGGATAGTTTTTTGCGACCACTCATCCCAGTTGTAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATTAGCAGATCGTTATGACCCAGTTCCGGCTCAGGAACGTCGGTCATCCAGATGccc  <  2:163124/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
GGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACGACAGTTGCGGCAATGACCACAGGTGATATGGCCTTCGCCAGAAACGCGATCGCCGATCTTGAAGCCTTTCACTTCCTGACCAATACCTACCACTTCACCGACATATTCATGGCCCACGACCATCGGCACCGGGATGGTTTTTTGCGACCACTCATCCCAGTTATAGATGTGAACGTCAGTCCCGCAGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGCCC  >  NC_000913/3790826‑3791313

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