Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,096,591 T→C 8.3% E116E (GAA→GAG garK ← glycerate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,096,5910TC 8.3% 163.9 / 4.9 60E116E (GAA→GAGgarKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base C (2/3);  total (28/32)
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.24e-01

GCGTGAAGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTG  >  minE/2096522‑2096661
                                                                     |                                                                      
gcgTGAAGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCAttt                                                                       <  1:1031675/71‑1 (MQ=255)
gcgTGAAGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCAttt                                                                       >  1:643035/1‑71 (MQ=255)
  gTGCAGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCACCAGCCCACTGGCCGCCGCCATTTCa                                                                     <  1:1717968/71‑1 (MQ=255)
     aaGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAAt                                                                   >  1:4348168/1‑70 (MQ=255)
     aaGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAAt                                                                   >  1:5198934/1‑70 (MQ=255)
       gTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCa                                                                     <  1:4720065/66‑1 (MQ=255)
       gTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCa                                                                     <  1:2078999/66‑1 (MQ=255)
            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)
                       cgcgTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCa                                                                     >  1:3487636/1‑50 (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)
                                ccGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATcc                                      >  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)
                                                                     |                                                                      
GCGTGAAGTGGTCACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTG  >  minE/2096522‑2096661

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