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

ATATGATGCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCATTATCACCATATGCAATAGA  >  CP000730/148354‑148629
                                                                                                                                         |                                                                                                                                           
atatGATGCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGtt                                                                                                                                         <  1:441808/141‑1 (MQ=255)
    gatgCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTCAGTGTGGTAGaa                                                                                                                                     >  1:83037/1‑141 (MQ=255)
     atgCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGaaa                                                                                                                                    <  1:65497/141‑1 (MQ=255)
     atgCAGTTGAACAAAGGATGAATACATAGCAATTTGATTAGGGTAGTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGaaa                                                                                                                                    >  2:458465/1‑141 (MQ=255)
       gtaGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATc                                                                                                                                  <  1:419312/139‑1 (MQ=255)
            tGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACAAAGGCAAGATATTCAGATTCTTAGAAATCATTaa                                                                                                                             <  2:424089/141‑1 (MQ=255)
               aCAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTCTCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTaaaaa                                                                                                                          >  2:115504/1‑141 (MQ=255)
               aCAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTaaaaa                                                                                                                          >  2:397210/1‑141 (MQ=255)
               aCAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAACATTGTTCGAAATCATTaaaaa                                                                                                                          >  2:317217/1‑141 (MQ=255)
                 aaaTCATGAATACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAat                                                                                                                        <  2:55439/141‑1 (MQ=255)
                  aaTCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAata                                                                                                                       >  1:209292/1‑141 (MQ=255)
                   aTCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAatat                                                                                                                      >  1:268975/1‑141 (MQ=255)
                      aTGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGACGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAAGAATTAAAAAAGataaa                                                                                                                   >  1:121178/1‑140 (MQ=255)
                          aGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATACTc                                                                                                               >  2:249344/1‑141 (MQ=255)
                           gACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCa                                                                                                              >  1:40712/1‑141 (MQ=255)
                            aCATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTATTGGTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCat                                                                                                             >  1:90010/1‑141 (MQ=255)
                                    aTTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACg                                                                                                     <  2:265577/141‑1 (MQ=255)
                                         aTTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTa                                                                                                >  2:409827/1‑141 (MQ=255)
                                         aTTATGTAATACATTTAGCAGCATTATTTAGTGTTGCTGAGTAGGTTGAGAAACCTATCTTATATCAAGAACTAAAAGTCGTAGCAACATTAAGATTGTTAGAAATCATTCAAAAATATCATAATCATATAAAACGTTGTa                                                                                                >  2:202458/1‑141 (MQ=255)
                                                            agcgTTAGTTTGTGTTGCTGAATCGGTTGAGAAACCTATCATATCTCAAGAAATATACGTCGTAGCAACATTAATTTTGTAAAAAAAAAAAAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCagcagc                                                                             <  2:264145/141‑1 (MQ=255)
                                                                ttagttagTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGtt                                                                         >  1:151118/1‑141 (MQ=255)
                                                                     tagtGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAAACAATAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATgg                                                                    <  2:135470/141‑1 (MQ=255)
                                                                               gAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCt                                                          >  1:124635/1‑141 (MQ=255)
                                                                               gAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCt                                                          <  1:397210/141‑1 (MQ=255)
                                                                               gAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGACTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCt                                                          <  1:115504/141‑1 (MQ=255)
                                                                                 gTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGa                                                        <  1:140372/141‑1 (MQ=255)
                                                                                      ttGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGc                                                   >  2:122251/1‑141 (MQ=255)
                                                                                      ttGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCAGTTTATGGAGAACTTCCTGATTTGc                                                   >  2:171630/1‑141 (MQ=255)
                                                                                       tGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGcc                                                  <  2:40712/141‑1 (MQ=255)
                                                                                       tGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGcc                                                  <  2:398190/141‑1 (MQ=255)
                                                                                       tGAGAAACCTAGCTTATCTCAAGAAATAAAACATCTGGTAAAAATACGAAGAAATAAAATCAATAAAAAAGATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGcc                                                  >  2:444346/1‑141 (MQ=255)
                                                                                          gaAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTaa                                               <  2:363306/141‑1 (MQ=255)
                                                                                                tatcttatctCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTCGATCTCCTTAAAAAATATAATAATCATATCAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGa                                         <  2:209292/141‑1 (MQ=255)
                                                                                                    ttatctCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGatcaa                                     <  1:409827/141‑1 (MQ=255)
                                                                                                    ttatctCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTACCTTTGCTTCGTCAGAAGCTGTTTATGGTGATATTCCTGATTTGCCTAAAAGTGataaa                                     >  2:304717/1‑141 (MQ=255)
                                                                                                       tctcAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGGGATCTTCCTGATTTGCCTAAAAGTGatcaatca                                  >  1:115323/1‑141 (MQ=255)
                                                                                                       actcAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGatcaatca                                  <  2:81940/140‑1 (MQ=255)
                                                                                                           aaGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTaa                              >  2:205725/1‑141 (MQ=255)
                                                                                                             aaaaTAAACGTCGTACCAACATTAAGACTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATc                            <  2:237424/140‑1 (MQ=255)
                                                                                                              aaaTAAACGTCGTAGCAACATTAAGATTGTTAAAAAAAATAAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCt                           <  2:251523/141‑1 (MQ=255)
                                                                                                                aTAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTa                         >  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)
                                                                                                                                         |                                                                                                                                           
ATATGATGCAGTTGAACAAATCATGAAGACATATCAATTTGATTATGTTATTCATTTAGCAGCATTAGTTAGTGTTGCTGAGTCGGTTGAGAAACCTATCTTATCTCAAGAAATAAACGTCGTAGCAACATTAAGATTGTTAGAAATCATTAAAAAATATAATAATCATATAAAACGTTTTATCTTTGCTTCGTCAGCAGCTGTTTATGGTGATCTTCCTGATTTGCCTAAAAGTGATCAATCATTAATCTTACCATTATCACCATATGCAATAGA  >  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: