Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 629,681 C→A L84I (CTA→ATA)  hxlA → possible 3‑hexulose‑6‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730629,6810CA94.4% 38.8 / ‑2.9 18L84I (CTA→ATA) hxlApossible 3‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (11/6);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

TTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTG  >  CP000730/629551‑629816
                                                                                                                                  |                                                                                                                                       
tttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                                                                               <  2:168697/141‑1 (MQ=255)
                    tGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                                                                           >  2:410971/1‑141 (MQ=255)
                    tGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                                                                           >  2:116185/1‑141 (MQ=255)
                    tGGAAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                                                                           >  2:304576/1‑141 (MQ=255)
                             aaataagaaaTATAAAAGTAATAGCAAACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTg                                                                                                  <  1:8365/133‑1 (MQ=255)
                                         tAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATa                                                                                      >  2:132049/1‑141 (MQ=255)
                                                                            gCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCt                                                   <  2:334909/141‑1 (MQ=255)
                                                                             cTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTg                                                  <  2:397389/141‑1 (MQ=255)
                                                                                      aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAtt                                         >  2:102696/1‑141 (MQ=255)
                                                                                         ttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAg                                      >  2:107538/1‑141 (MQ=255)
                                                                                            gCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGaaa                                   >  2:82198/1‑141 (MQ=255)
                                                                                            gCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGATGATATGCTTGCGGTTCAAGATTTAGaaa                                   >  2:41401/1‑141 (MQ=255)
                                                                                              cAAGCAATTAAACTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGATGATGCATCAATTCAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTTAAGATTTAGaaaaa                                 <  2:340069/141‑1 (MQ=255)
                                                                                                    aTTAAATTTGGCGCGGATGTAATTACAATAAGAGGTGTTGCAGAAGATGAATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCa                           <  2:132458/141‑1 (MQ=255)
                                                                                                                 cgGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATg              >  2:434782/1‑141 (MQ=255)
                                                                                                                   gATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGaa            >  2:163864/1‑141 (MQ=255)
                                                                                                                     tGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGAGGaaag          >  2:460308/1‑140 (MQ=255)
                                                                                                                         aTTACAATAATAGGTGTGGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGt      <  2:413636/141‑1 (MQ=255)
                                                                                                                             aaataatAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTg  <  2:107035/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
TTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTG  >  CP000730/629551‑629816

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