Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 629,681 C→A 100% 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,6810CA88.1% 65.5 / 1.1 25L84I (CTA→ATA) hxlApossible 3‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (8/14);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGA  >  CP000730/629549‑629817
                                                                                                                                    |                                                                                                                                        
ggTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGtt                                                                                                                                  <  1:226427/141‑1 (MQ=255)
ggTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGtt                                                                                                                                  >  1:155064/1‑141 (MQ=255)
  tttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGc                                                                                                                                <  1:197111/141‑1 (MQ=255)
   ttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCa                                                                                                                               >  1:49773/1‑141 (MQ=255)
          cagTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATg                                                                                                                        <  1:25616/141‑1 (MQ=255)
            gTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCa                                                                                                                      <  2:107764/141‑1 (MQ=255)
              tAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATc                                                                                                                    <  1:37899/141‑1 (MQ=255)
                  cATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAtt                                                                                                                <  1:326225/141‑1 (MQ=255)
                      tGGCAGACAACATTAGTNATGTAAAAGTATTAGCAGANATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACTATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                                                                            >  1:233353/1‑141 (MQ=255)
                               acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGTCTAATAAGATAACCAAGATAGAAAAATTGGCGCGTATGTAAGTTCAATAATTGGTGTTGCAGGAGATGCATCAAGGAAAGCAGCGATTg                                                                                                   >  1:197819/1‑141 (MQ=255)
                                          gTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCTAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                                                                        >  2:153916/1‑141 (MQ=255)
                                                            aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTa                                                                      >  2:183117/1‑141 (MQ=255)
                                                                                tGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATTATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGt                                                  <  1:24176/141‑1 (MQ=255)
                                                                                    tATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCaa                                              <  1:135084/141‑1 (MQ=255)
                                                                                         aGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAttt                                         >  2:292422/1‑141 (MQ=255)
                                                                                          gTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTAATATGATTGCTGTTCAAGATTTa                                        <  1:48171/141‑1 (MQ=255)
                                                                                                cAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGaaaaa                                  <  1:151572/141‑1 (MQ=255)
                                                                                                   gCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGt                               <  1:264952/141‑1 (MQ=255)
                                                                                                     aaTTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGc                             <  2:224912/141‑1 (MQ=255)
                                                                                                        tAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaaa                          <  1:27910/141‑1 (MQ=255)
                                                                                                             ttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAAc                     <  1:253491/141‑1 (MQ=255)
                                                                                                              tGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACt                    >  2:287274/1‑141 (MQ=255)
                                                                                                                      aTGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGaaa            >  1:138567/1‑141 (MQ=255)
                                                                                                                       tGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAAt           <  1:321909/141‑1 (MQ=255)
                                                                                                                                aaTACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGa  >  1:139991/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                        
GGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGA  >  CP000730/629549‑629817

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