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,6810CA96.8% 107.1 / ‑6.2 31L84I (CTA→ATA) hxlApossible 3‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (14/16);  minor base G (1/0);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGT  >  CP000730/629551‑629812
                                                                                                                                  |                                                                                                                                   
tttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                                                                           <  1:123664/141‑1 (MQ=255)
  tACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCag                                                                                                                         <  2:46536/141‑1 (MQ=255)
      agcagTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCagaaga                                                                                                                     <  2:11801/141‑1 (MQ=255)
           ttAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCAt                                                                                                                >  2:161745/1‑141 (MQ=255)
                           caacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAt                                                                                                >  1:273364/1‑141 (MQ=255)
                             acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTg                                                                                              <  1:226944/141‑1 (MQ=255)
                                     aaTGAAAAAGTATTAGCAGACATGAAAAATATAGATGCAGCTGCTTATGAAGTTAGCAAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCt                                                                                      <  2:7203/141‑1 (MQ=255)
                                       tGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCa                                                                                    >  2:83542/1‑141 (MQ=255)
                                         tAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATa                                                                                  >  1:257604/1‑141 (MQ=255)
                                                  tAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                                                                         <  1:200029/141‑1 (MQ=255)
                                                   aGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAAtaataa                                                                        >  1:112568/1‑141 (MQ=255)
                                                    gCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATaaa                                                                       >  2:174578/1‑141 (MQ=255)
                                                          aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTa                                                                 >  2:237807/1‑141 (MQ=255)
                                                           tGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTAc                                                                <  1:161745/141‑1 (MQ=255)
                                                            gAAACTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAAATAAAATTGGAGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACt                                                               >  2:232651/1‑141 (MQ=255)
                                                                     aGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGatat                                                      >  1:294383/1‑141 (MQ=255)
                                                                      gATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCCCATAAAAATAATAAACAATTACTAGTTGATATg                                                     >  2:334610/1‑141 (MQ=255)
                                                                                  tATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCaa                                         >  2:276077/1‑141 (MQ=255)
                                                                                   aggAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAGTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAATGCTGTTCAAg                                        >  2:234009/3‑141 (MQ=255)
                                                                                   aTGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAg                                        >  2:214258/1‑141 (MQ=255)
                                                                                   aTGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAg                                        >  1:359059/1‑141 (MQ=255)
                                                                                      aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAtt                                     <  1:202862/141‑1 (MQ=255)
                                                                                      aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAtt                                     <  1:110921/141‑1 (MQ=255)
                                                                                                   aaTTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGc                        <  2:349381/141‑1 (MQ=255)
                                                                                                   aaTTAAATTTGGCGCGGATGTAATTAAAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAAAAATAAACAATTACTAGTTGATATGATTGCTGTTAAAAAAAAAAAAAAACGTGc                        <  2:125129/141‑1 (MQ=255)
                                                                                                     ttAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaa                      <  1:107835/141‑1 (MQ=255)
                                                                                                            tGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACt               >  2:274379/1‑141 (MQ=255)
                                                                                                             ggCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTa              <  2:257604/141‑1 (MQ=255)
                                                                                                                  ggATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGa         <  1:237807/141‑1 (MQ=255)
                                                                                                                         aTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGt  <  1:296828/141‑1 (MQ=255)
                                                                                                                         aTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGt  <  1:296820/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
TTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGT  >  CP000730/629551‑629812

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