Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,096,632 T→A 8.2% S103C (AGT→TGT) ‡ garK ← glycerate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,096,6320TA 8.2% 166.3 / 6.3 61S103C (AGT→TGT) ‡garKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base A (2/3);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCA  >  minE/2096564‑2096699
                                                                    |                                                                   
cAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                                                   <  1:3979976/71‑1 (MQ=255)
      cAGCCCACTGGCGGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCa                                                             <  1:2768444/71‑1 (MQ=255)
          ccACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACttt                                                         <  1:506739/71‑1 (MQ=255)
             cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtctc                                                      >  1:3926480/1‑71 (MQ=255)
              tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                                                     >  1:4799295/1‑71 (MQ=255)
              tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                                                     >  1:3593843/1‑71 (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)
                    ccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCa                                                 <  1:526871/69‑1 (MQ=255)
                     cgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCATCg                                              <  1:5037709/71‑1 (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)
                                  aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCa                                 <  1:4018471/71‑1 (MQ=255)
                                    cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                                           >  1:4705379/1‑43 (MQ=255)
                                    cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                                           >  1:4415560/1‑43 (MQ=255)
                                     gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                                       >  1:922753/1‑46 (MQ=255)
                                     gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                                       >  1:1380603/1‑46 (MQ=255)
                                      cgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCg                             <  1:2474521/71‑1 (MQ=255)
                                          tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACttt                                                         <  1:2655758/39‑1 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAg                                                >  1:976787/1‑47 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGcc                                           >  1:4816747/1‑52 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa                          <  1:660612/69‑1 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa                          <  1:4347880/69‑1 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGAc                         <  1:1676490/70‑1 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGAc                         <  1:3019747/70‑1 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg                        >  1:3011316/1‑71 (MQ=255)
                                           ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTATCGCCCAGCGGCCCTGTAACCCAGGCGTGa                          <  1:5476708/69‑1 (MQ=255)
                                           ttGCCAGCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg                        >  1:737594/1‑71 (MQ=255)
                                            tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGt                       >  1:1211212/1‑71 (MQ=255)
                                            tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGt                       >  1:5274956/1‑71 (MQ=255)
                                             gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg                        >  1:592681/1‑69 (MQ=255)
                                             gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGtt                      >  1:4407448/1‑71 (MQ=255)
                                                 tCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCg                             <  1:102826/60‑1 (MQ=255)
                                                   gCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTaa                                     >  1:2216563/1‑50 (MQ=255)
                                                      ggAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCt             <  1:4734987/71‑1 (MQ=255)
                                                      ggAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCt             <  1:1456238/71‑1 (MQ=255)
                                                       gagaTCCCCCAACTGGCATTCACTTTCTCGCCCATCGGCCCTGt                                       <  1:5020226/44‑1 (MQ=255)
                                                       gagaTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTg            >  1:3797718/1‑71 (MQ=255)
                                                         gaTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa                          >  1:2952483/1‑55 (MQ=255)
                                                          aTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGt         <  1:1443069/71‑1 (MQ=255)
                                                              cccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGt                            <  1:851536/48‑1 (MQ=255)
                                                              cccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGt                            <  1:4788619/48‑1 (MQ=255)
                                                              cccAACTGGCATTCACTTTCTCGCCCAGAGGCCCTGTAACCCAGGCGt                            <  1:3133575/48‑1 (MQ=255)
                                                                cAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTaa                                     >  1:112187/1‑37 (MQ=255)
                                                                cAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGc   <  1:3056159/71‑1 (MQ=255)
                                                                 aaCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCATCCCCCTGGGTGGCTGCa  <  1:5485661/71‑1 (MQ=255)
                                                                  aCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCa                    >  1:4794662/1‑52 (MQ=255)
                                                                  aCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCa                    >  1:2803914/1‑52 (MQ=255)
                                                                  aCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGccc                <  1:4021149/56‑1 (MQ=255)
                                                                   cTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAccc                                  >  1:2997395/1‑37 (MQ=255)
                                                                    |                                                                   
CAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCA  >  minE/2096564‑2096699

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