Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,494,331 A→G 100% V113V (GTA→GTG ibpA → heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,494,3310AG81.4% 99.5 / 12.0 64V113V (GTA→GTGibpAheat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (3/9);  new base G (0/52);  total (3/61)
Fisher's exact test for biased strand distribution p-value = 5.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGCTGTATATCGATCTCGAACGCGTGAT  >  minE/2494263‑2494366
                                                                    |                                   
gCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACGTGGTa                                     >  1:1874063/1‑69 (MQ=255)
gCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTa                                     >  1:884788/1‑69 (MQ=255)
gCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTa                                     >  1:2123563/1‑69 (MQ=255)
gCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTa                                     >  1:1153449/1‑69 (MQ=255)
gCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTa                                     >  1:2707769/1‑69 (MQ=255)
gCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTa                                     >  1:12619/1‑69 (MQ=255)
      cGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTaaa                                   <  1:1489775/65‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:3209615/66‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:3139133/66‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:3398715/66‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:111023/66‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:1355661/66‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:731889/66‑1 (MQ=255)
             ttGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:1311878/66‑1 (MQ=255)
                   gCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGCTGTATATCg                  <  1:2246706/69‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTTCTGTATATCGATCTCGAAc         <  1:2145159/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2470830/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:3158321/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:3238434/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:3130527/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2938271/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:3258053/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:271027/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:848692/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2688506/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:25871/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2509396/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2496558/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:1030271/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:3376702/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:378658/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2070315/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2021994/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2018113/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:1981571/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:1502923/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:2192341/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:1901243/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAAc         <  1:1612116/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCAGCTGTATATCGATCTCGAAc         <  1:2439223/65‑1 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGCTGTATATCGATCTCGAACGCGt     >  1:3368880/1‑69 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGCTGTATATCGATCTCGAACGCGt     >  1:1448808/1‑69 (MQ=255)
                                aGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGCTGTATATCGATCTCGAACGCGt     >  1:1679795/1‑69 (MQ=255)
                                 gCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGc                           <  1:1961136/46‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTTTTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1760590/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:706447/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:665187/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:602162/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1521709/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1257270/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:567861/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:802622/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:817144/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:500702/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1173193/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1466038/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1089280/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:3367822/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:3313112/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:3293862/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1494605/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1545526/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:3042379/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:3009332/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:2988773/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:2788468/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1910193/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:2623882/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:2323544/69‑1 (MQ=255)
                                   tGAGAACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:2257952/69‑1 (MQ=255)
                                    gagaACATTCATGTTCGTGGTGCTAACCTGGTGAATGGTCTGCTGTATATCGATCTCGAACGCGTGAt  <  1:1466336/68‑1 (MQ=255)
                                                                    |                                   
GCTGAACGCAACTTTGAACGCAAATTCCAGTTAGCTGAGAACATTCATGTTCGTGGTGCTAACCTGGTAAATGGTTTGCTGTATATCGATCTCGAACGCGTGAT  >  minE/2494263‑2494366

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