Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,096,597 A→C 8.2% F114L (TTT→TTG) ‡ garK ← glycerate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,096,5970AC 8.2% 149.7 / 7.4 61F114L (TTT→TTG) ‡garKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base C (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAC  >  minE/2096534‑2096665
                                                               |                                                                    
cACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATaaa                                                                    >  1:4661834/1‑66 (MQ=255)
 aCGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGt                                                              <  1:1169047/71‑1 (MQ=255)
    aGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGttt                                                            >  1:3821196/1‑70 (MQ=255)
      tGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                         <  1:522963/71‑1 (MQ=255)
      tGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                         <  1:4657247/71‑1 (MQ=255)
       ggATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                         <  1:979523/70‑1 (MQ=255)
       ggATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                         >  1:2143747/1‑70 (MQ=255)
        gATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGcc                                                        <  1:5144728/70‑1 (MQ=255)
          tCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATc                                                     >  1:146437/1‑71 (MQ=255)
                ttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGcc                                                  >  1:2997059/1‑68 (MQ=255)
                 tttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCg                                                 >  1:3010229/1‑68 (MQ=255)
                 tttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGga                                               <  1:5280595/70‑1 (MQ=255)
                  ttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTg                                                          >  1:1095900/1‑58 (MQ=255)
                  ttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCg                                                 >  1:2439995/1‑67 (MQ=255)
                   tCCGCAGGTACCTGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGgaga                                             <  1:1582598/70‑1 (MQ=255)
                    ccGCAGGTACCAGCTCCAGCCCACTGGCCGCCG‑CATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATcc                                          >  1:787884/1‑71 (MQ=255)
                     cGCATGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATc                                           <  1:3388872/70‑1 (MQ=255)
                      gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGt                                                              <  1:675312/50‑1 (MQ=255)
                      gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGt                                                              <  1:1759723/50‑1 (MQ=255)
                      gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccc                                         >  1:4005276/1‑71 (MQ=255)
                           tACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                    <  1:602647/71‑1 (MQ=255)
                             ccAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGcc                                                        <  1:814367/49‑1 (MQ=255)
                              cAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                 <  1:3979976/71‑1 (MQ=255)
                               aGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGcc                                                        <  1:4865336/47‑1 (MQ=255)
                                 cTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATcc                                          >  1:1091799/1‑59 (MQ=255)
                                 cTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATc                                           >  1:1253572/1‑58 (MQ=255)
                                  tCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                     >  1:2957591/1‑63 (MQ=255)
                                    cAGCCCACTGGCGGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCa                           <  1:2768444/71‑1 (MQ=255)
                                        ccACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACttt                       <  1:506739/71‑1 (MQ=255)
                                          aCTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGgaga                                             <  1:755891/47‑1 (MQ=255)
                                          aCTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCa                                      <  1:2198443/54‑1 (MQ=255)
                                           cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                       <  1:4667717/52‑1 (MQ=255)
                                           cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                       <  1:1745855/52‑1 (MQ=255)
                                           cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtctc                    >  1:3926480/1‑71 (MQ=255)
                                           cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCCGAGATccccc                                       <  1:179563/52‑1 (MQ=255)
                                            tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                   >  1:3593843/1‑71 (MQ=255)
                                            tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                   >  1:4799295/1‑71 (MQ=255)
                                             ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                     >  1:3957323/1‑52 (MQ=255)
                                             ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                     >  1:4515923/1‑52 (MQ=255)
                                             ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                    >  1:2219952/1‑53 (MQ=255)
                                             ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGc                  >  1:2411976/1‑71 (MQ=255)
                                              gctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                 >  1:786951/4‑71 (MQ=255)
                                              gctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                 >  1:697326/4‑71 (MQ=255)
                                               ctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                 <  1:5535532/68‑1 (MQ=255)
                                               ccgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                 <  1:395542/54‑1 (MQ=255)
                                               ccgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGccc                <  1:3151299/71‑1 (MQ=255)
                                                  ccgccATTTCAATAAACGCGGTTTTGCCATCGCCg                                                 <  1:4194555/35‑1 (MQ=255)
                                                  ccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCa               <  1:526871/69‑1 (MQ=255)
                                                   cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                       <  1:5388770/44‑1 (MQ=255)
                                                   cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                       <  1:3167975/44‑1 (MQ=255)
                                                   cgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCATCg            <  1:5037709/71‑1 (MQ=255)
                                                    gccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATc                                           >  1:605773/1‑39 (MQ=255)
                                                    gccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                    >  1:5606833/1‑46 (MQ=255)
                                                     ccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGc          >  1:433689/1‑71 (MQ=255)
                                                        tttCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCt       >  1:5273455/1‑71 (MQ=255)
                                                         ttCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTg      >  1:5572796/1‑71 (MQ=255)
                                                          tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCAtt                             <  1:5260805/47‑1 (MQ=255)
                                                          tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                     >  1:5192297/1‑55 (MQ=255)
                                                          tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtc                      >  1:3790962/1‑54 (MQ=255)
                                                             aTAAACGCGGTTTTGACATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAc  >  1:3663530/1‑71 (MQ=255)
                                                             aTAAACGCGGTTTTGACATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAc  >  1:3670230/1‑71 (MQ=255)
                                                              tAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAg              >  1:3800223/1‑58 (MQ=255)
                                                               |                                                                    
CACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAC  >  minE/2096534‑2096665

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