Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,238,307 C→T 100% Q141* (CAA→TAA)  USA300HOU_1162 → glycerone kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,238,3070CT100.0% 52.6 / NA 18Q141* (CAA→TAA) USA300HOU_1162glycerone kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/7);  total (11/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAAAAATATTGAAAGTGAATCTGAAATTAATTCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGCAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACTTATTAGCTGTACTTAAAGAAGTTGGTGTTGTT  >  CP000730/1238170‑1238438
                                                                                                                                         |                                                                                                                                   
gTAAAAATATTGAAAGTGAATCTGAAATTAATTCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAg                                                                                                                                  >  2:46723/1‑141 (MQ=255)
         ttGAAAGTGAATCTGAAATTAATTCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAg                                                                                                                         >  1:117193/1‑141 (MQ=255)
          tGAAAGTGAATCTGATATTAATTCAAGATTGTTATCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGCAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGa                                                                                                                        <  1:173823/141‑1 (MQ=255)
                           ttaattCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATAcag                                                                                                       >  2:146101/1‑139 (MQ=255)
                           ttaattCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAACCTGCAATAGAAAAAGCAAATAATACTg                                                                                                       >  2:144982/1‑141 (MQ=255)
                             aattCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACGGaa                                                                                                     <  1:26391/141‑1 (MQ=255)
                                     aTTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCCATAGAAAAAGCAAATAATACTGAAGATTGtat                                                                                             <  1:222533/141‑1 (MQ=255)
                                      ttGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGtata                                                                                            <  2:284887/141‑1 (MQ=255)
                                                       tCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACat                                                                           >  1:188056/1‑141 (MQ=255)
                                                                tgtTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTaaa                                                                  >  2:248461/1‑141 (MQ=255)
                                                                 gtTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTaaaa                                                                 <  1:132080/141‑1 (MQ=255)
                                                                         ggCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGa                                                         >  1:282135/1‑141 (MQ=255)
                                                                           cATATAAAGCTGTTATGAAACCAGTCGACAGTACAATACTTACAATGGCAAAAGATGCTACGTAAGCTGCAATAGAAACGACAACAACGCCGGAACCACGTCTAGAATGAATGGAGGAAATTATTGTaaagagcgatgcac                                                       >  2:104441/1‑130 (MQ=255)
                                                                                                   ttGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACttat                               >  2:169033/1‑141 (MQ=255)
                                                                                                               tACTTACAGTTGCAAAAGATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACTTATTAGCTGTACTTa                   <  1:46723/141‑1 (MQ=255)
                                                                                                                              aaGATGCTGCGTAAGCTGCAATAGAAAAAGCAAGTAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACTTATTAGCTGTACTTAAAGAAGTTGGtgttt    >  1:13576/1‑140 (MQ=255)
                                                                                                                                gATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACTTATTAGCTGTACTTAAAGAAGTTGGTgttgtt  <  1:548/141‑1 (MQ=255)
                                                                                                                                gATGCTGCGTAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACTTATTAGCTGTACTTAAAGAAGTTGGTgttgtt  >  2:548/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                   
GTAAAAATATTGAAAGTGAATCTGAAATTAATTCAAAATTGTTAGCTGAAAGTTTTCAAGCTGGTGTTGAAACGGCATATAAAGCTGTTATGAAACCAGTTGAAGGTACAATACTTACAGTTGCAAAAGATGCTGCGCAAGCTGCAATAGAAAAAGCAAATAATACTGAAGATTGTATAGAATTAATGGAGTACATTATTGTAAAAGCCAATGAATCACTTGAAAACACACCAAACTTATTAGCTGTACTTAAAGAAGTTGGTGTTGTT  >  CP000730/1238170‑1238438

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