Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	982,632	A→C	8.5%	A26A (GCA→GCC)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	982,632	0	A	C	8.5%	226.4 / 19.2	94	A26A (GCA→GCC)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/46);  new base C (2/6);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

CTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCA  >  minE/982562‑982701                                                                       |                                                                      cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                                                       <  1:2137121/71‑1 (MQ=255)  tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                      >  1:617461/1‑71 (MQ=255)   aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                      >  1:3538746/1‑70 (MQ=255)   aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                     >  1:1185113/1‑71 (MQ=255)    tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                     <  1:360413/70‑1 (MQ=255)    tCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCAt                                                                     >  1:368021/1‑70 (MQ=255)      gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                                  >  1:2597968/1‑71 (MQ=255)      gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                                  >  1:2418339/1‑71 (MQ=255)         cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               >  1:690803/1‑71 (MQ=255)          ttGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                                  >  1:1072219/1‑67 (MQ=255)          ttGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               >  1:3167414/1‑70 (MQ=255)            ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                 <  1:1254250/66‑1 (MQ=255)            ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               <  1:3304546/68‑1 (MQ=255)            ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               <  1:2625073/68‑1 (MQ=255)            ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               <  1:3724591/68‑1 (MQ=255)            ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGt                                                                <  1:64645/67‑1 (MQ=255)             gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                                           <  1:1647967/71‑1 (MQ=255)                acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                         >  1:3407328/1‑70 (MQ=255)                   acGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                 <  1:1649955/59‑1 (MQ=255)                    cGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCa                                                    <  1:4121568/71‑1 (MQ=255)                     gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                     <  1:2877200/53‑1 (MQ=255)                     gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                     <  1:1924328/53‑1 (MQ=255)                     gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                   <  1:2559127/71‑1 (MQ=255)                         ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                               >  1:79402/1‑71 (MQ=255)                           ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                      >  1:3436327/1‑46 (MQ=255)                            tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                               <  1:2976596/68‑1 (MQ=255)                              gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                          <  1:3445543/71‑1 (MQ=255)                              gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                           >  1:2074359/1‑70 (MQ=255)                              gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                           >  1:1749456/1‑70 (MQ=255)                              gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                           >  1:1961839/1‑70 (MQ=255)                                 tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                                              <  1:2939933/64‑1 (MQ=255)                                 tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                          <  1:4012211/68‑1 (MQ=255)                                 tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                          <  1:4132219/68‑1 (MQ=255)                                   cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                                  >  1:2091483/1‑58 (MQ=255)                                    ggCGTGTTTCCAAACCAGGACTGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGa                                    <  1:902362/71‑1 (MQ=255)                                       gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                 <  1:3465116/71‑1 (MQ=255)                                       gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                 <  1:3696420/71‑1 (MQ=255)                                       gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                 <  1:3312984/71‑1 (MQ=255)                                       gtgtTTCCAAACAAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                                 <  1:2902657/71‑1 (MQ=255)                                       gtgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGaaac                                 <  1:3423917/71‑2 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               <  1:3185363/40‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                         <  1:1403516/46‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                         <  1:411044/46‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:2966535/71‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:3095275/71‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:857174/71‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:2794418/71‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:1176852/71‑1 (MQ=255)                                        tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:1925809/71‑1 (MQ=255)                                        tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAAcat                                              <  1:2597139/57‑3 (MQ=255)                                        tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAAcat                                              <  1:3212676/57‑3 (MQ=255)                                         gtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                                   <  1:660993/67‑1 (MQ=255)                                          tttCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                 <  1:113154/36‑1 (MQ=255)                                          tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                                   <  1:2608187/66‑1 (MQ=255)                                          tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATg                               <  1:3325922/70‑1 (MQ=255)                                          tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACCTTTAAAACCGaaaa                                 <  1:3151060/68‑1 (MQ=255)                                           ttCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                 <  1:3806118/35‑1 (MQ=255)                                           ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                   >  1:1860149/1‑49 (MQ=255)                                           ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                   >  1:3195584/1‑49 (MQ=255)                                           ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGt                             >  1:1701101/1‑71 (MQ=255)                                            tCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                               >  1:928099/1‑36 (MQ=255)                                            tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                                 >  1:2743020/1‑50 (MQ=255)                                            tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                   >  1:2187549/1‑48 (MQ=255)                                            tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                              >  1:2943252/1‑69 (MQ=255)                                            tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                              >  1:2938365/1‑69 (MQ=255)                                              ccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcg                          <  1:2021901/69‑2 (MQ=255)                                              ccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcg                          <  1:3405802/69‑2 (MQ=255)                                              cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:1902389/65‑1 (MQ=255)                                              cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:2615246/65‑1 (MQ=255)                                              cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                          <  1:3279710/71‑1 (MQ=255)                                              cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                          >  1:3640925/1‑71 (MQ=255)                                              aaaaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                          >  1:2602351/2‑71 (MQ=255)                                                aaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCg                        >  1:2528531/1‑71 (MQ=255)                                                 aCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                                             >  1:2418317/1‑49 (MQ=255)                                                 aCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGt                       >  1:4142994/1‑71 (MQ=255)                                                  ccAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcgcc                      <  1:1982902/71‑3 (MQ=255)                                                  ccAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                      >  1:2620667/1‑71 (MQ=255)                                                  ccAGGACGGCACGTTCACCGCAATGACGTATACAAAAAGCAAAACg                                               >  1:1177457/1‑46 (MQ=255)                                                     ggACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgct                   <  1:1588817/71‑1 (MQ=255)                                                      gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctg                  <  1:1125718/71‑1 (MQ=255)                                                      gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctg                  >  1:3832337/1‑71 (MQ=255)                                                         ggCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                                <  1:3469862/54‑1 (MQ=255)                                                          gaaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg              >  1:1052884/3‑71 (MQ=255)                                                           cACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                                 >  1:740026/1‑35 (MQ=255)                                                           cACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                                 >  1:1684374/1‑35 (MQ=255)                                                           cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                              >  1:1660293/1‑54 (MQ=255)                                                           cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctggctg              >  1:3328170/1‑70 (MQ=255)                                                            aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga            <  1:2083453/71‑1 (MQ=255)                                                             cGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa           <  1:1497759/71‑1 (MQ=255)                                                              gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                              >  1:1326102/1‑51 (MQ=255)                                                                 cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAAc                                       >  1:2146386/1‑39 (MQ=255)                                                                  aCCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa      >  1:708267/1‑71 (MQ=255)                                                                  aCCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa      >  1:3224748/1‑71 (MQ=255)                                                                   ccGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGCTGCCCGTCGCTGGCTGGAAAGAAAc     >  1:907464/1‑71 (MQ=255)                                                                     gCAATGACGTATACCAAAAGCAAAACGTTTAAAAc                                       >  1:2543582/1‑35 (MQ=255)                                                                     gCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt    >  1:360161/1‑70 (MQ=255)                                                                     gCAATGACGTATACCAAAAGCAAAACGGTTAAAAc                                       >  1:428864/1‑35 (MQ=255)                                                                      cAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCa  <  1:3163031/71‑1 (MQ=255)                                                                       aaTGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa      >  1:533357/1‑66 (MQ=255)                                                                       aaTGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga        <  1:2783499/64‑1 (MQ=255)                                                                       |                                                                      CTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCA  >  minE/982562‑982701

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