Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 148,490 T→C 11.3% L105L (TTG→CTG)  galE → UDP‑glucose 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730148,4900TC11.3% 109.8 / 4.9 47L105L (TTG→CTG) galEUDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (2/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ATATGATGCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCATTATCACCATATGCAATAGA  >  CP000730/148354‑148629
                                                                                                                                         |                                                                                                                                           
atatGATGCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGtt                                                                                                                                         <  1:441808/141‑1 (MQ=255)
    gatgCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTCAGTGTGGTAGaa                                                                                                                                     >  1:83037/1‑141 (MQ=255)
     atgCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGaaa                                                                                                                                    <  1:65497/141‑1 (MQ=255)
     atgCAGTTGAACAAAGGATGAATACATAGCAATTTGATTAGGGTAGTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGaaa                                                                                                                                    >  2:458465/1‑141 (MQ=255)
       gtaGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATc                                                                                                                                  <  1:419312/139‑1 (MQ=255)
            tGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACAAAGGCAAGATATTCAGATTCTTAGAAATCATTaa                                                                                                                             <  2:424089/141‑1 (MQ=255)
               aCAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTCTCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTaaaaa                                                                                                                          >  2:115504/1‑141 (MQ=255)
               aCAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTaaaaa                                                                                                                          >  2:397210/1‑141 (MQ=255)
               aCAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAACATTGTTCGAAATCATTaaaaa                                                                                                                          >  2:317217/1‑141 (MQ=255)
                 aaaTCATGAATACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAat                                                                                                                        <  2:55439/141‑1 (MQ=255)
                  aaTCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAata                                                                                                                       >  1:209292/1‑141 (MQ=255)
                   aTCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAatat                                                                                                                      >  1:268975/1‑141 (MQ=255)
                      aTGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGACGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAAGAATTAAAAAAGataaa                                                                                                                   >  1:121178/1‑140 (MQ=255)
                          aGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATACTc                                                                                                               >  2:249344/1‑141 (MQ=255)
                           gACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCa                                                                                                              >  1:40712/1‑141 (MQ=255)
                            aCATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTATTGGTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCat                                                                                                             >  1:90010/1‑141 (MQ=255)
                                    aTTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACg                                                                                                     <  2:265577/141‑1 (MQ=255)
                                         aTTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTa                                                                                                >  2:409827/1‑141 (MQ=255)
                                         aTTATGTAATACATTTAGCAGCATTATTTAGTGTTGCTGAGTAGGTTGAGAAACCTATCTTATATCAAGAACT‑AAAAGTCGTAGCAACATTAAGATTGTTAGAAATCATTCAAAAATATCATAATCATATAAAACGTTGTa                                                                                                >  2:202458/1‑141 (MQ=255)
                                                            agcgTTAGTTTGTGTTGCTGAATCGGTTGAGAAACCTATCATATCTCAAGAAAT‑ATACGTCGTAGCAACATTAATTTTGTAAAAAAAAAAAAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCagcagc                                                                             <  2:264145/141‑1 (MQ=255)
                                                                ttagttagTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGtt                                                                         >  1:151118/1‑141 (MQ=255)
                                                                     tagtGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAAACAATAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATgg                                                                    <  2:135470/141‑1 (MQ=255)
                                                                               gAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCt                                                          >  1:124635/1‑141 (MQ=255)
                                                                               gAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCt                                                          <  1:397210/141‑1 (MQ=255)
                                                                               gAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCt                                                          <  1:115504/141‑1 (MQ=255)
                                                                                 gTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGa                                                        <  1:140372/141‑1 (MQ=255)
                                                                                      ttGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGc                                                   >  2:122251/1‑141 (MQ=255)
                                                                                      ttGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCAGTTTATGGAGAACTTCCTGATTTGc                                                   >  2:171630/1‑141 (MQ=255)
                                                                                       tGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGcc                                                  <  2:40712/141‑1 (MQ=255)
                                                                                       tGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGcc                                                  <  2:398190/141‑1 (MQ=255)
                                                                                       tGAGAAACCTAGCTTATCTCAAGAAATAAAACATC‑TGGTAAAAATACGAAGAAATAAAATCAATAAAAAAGATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGcc                                                  >  2:444346/1‑141 (MQ=255)
                                                                                          gaAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTaa                                               <  2:363306/141‑1 (MQ=255)
                                                                                                tatcttatctCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTCGATCTCCTTAAAAAATATAATAATCATATCAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGa                                         <  2:209292/141‑1 (MQ=255)
                                                                                                    ttatctCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGatcaa                                     <  1:409827/141‑1 (MQ=255)
                                                                                                    ttatctCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTACCTTTGCTTCGTCAGAAGCTGTTTATGGTGATATTCCTGATTTGCCTAAAAGTGataaa                                     >  2:304717/1‑141 (MQ=255)
                                                                                                       tctcAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGGGATCTTCCTGATTTGCCTAAAAGTGatcaatca                                  >  1:115323/1‑141 (MQ=255)
                                                                                                       actcAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGatcaatca                                  <  2:81940/140‑1 (MQ=255)
                                                                                                           aaGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTaa                              >  2:205725/1‑141 (MQ=255)
                                                                                                             aaaaT‑AAACGTCGTACCAACATTAAGACTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATc                            <  2:237424/140‑1 (MQ=255)
                                                                                                              aaaT‑AAACGTCGTAGCAACATTAAGATTGTTAAAAAAAATAAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCt                           <  2:251523/141‑1 (MQ=255)
                                                                                                                aT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTa                         >  1:381814/1‑141 (MQ=255)
                                                                                                                   aaaCGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTAcc                       <  2:421616/141‑1 (MQ=255)
                                                                                                                     aCGTCGTAGCAACTGTAAGATTGGTAGAAATCAGTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCAt                     >  2:458476/1‑141 (MQ=255)
                                                                                                                          gtAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCATTATCa                <  1:363895/141‑1 (MQ=255)
                                                                                                                           tAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGCGCTATCATTAATCTTACCATTATCAc               <  1:458465/141‑1 (MQ=255)
                                                                                                                                 cATTAAGATCGTTATAAATCATTAAAAAATATCATACCTATATATTTCGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGACCTTCCTGATTTGCCTAAAAGTGATCCATCATTCATCTTACCACTATCACCATATg         <  2:95910/141‑1 (MQ=255)
                                                                                                                                        aTTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCATTATCACCATATGCAATAGa  <  1:143618/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
ATATGATGCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAAT‑AAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCATTATCACCATATGCAATAGA  >  CP000730/148354‑148629

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: ‑