Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,976,997 G→A 6.0% L83L (TTG→TTA sprT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,976,9970GA 6.0% 223.0 / 3.4 67L83L (TTG→TTAsprTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

ATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTG  >  minE/1976927‑1977066
                                                                      |                                                                     
aTCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTg                                                                       <  1:4715035/71‑1 (MQ=255)
           cTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGaa                                                            <  1:3348623/71‑1 (MQ=255)
              ttGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACa                                                         >  1:4306076/1‑71 (MQ=255)
              ttGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACa                                                         >  1:21727/1‑71 (MQ=255)
               tGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACAt                                                        >  1:2574149/1‑71 (MQ=255)
               tGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACAt                                                        >  1:2199815/1‑71 (MQ=255)
                ggAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACAt                                                        <  1:152428/70‑1 (MQ=255)
                   aaaCAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCg                                                    >  1:3881665/1‑71 (MQ=255)
                       aGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGaa                                                            <  1:4093054/59‑1 (MQ=255)
                       aGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGaa                                                            <  1:4024303/59‑1 (MQ=255)
                       aGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGaa                                                            <  1:3591625/59‑1 (MQ=255)
                           aaGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCcgcg                                              >  1:111112/1‑69 (MQ=255)
                           aaGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCcgcg                                              >  1:816697/1‑69 (MQ=255)
                           aaGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGt                                             <  1:1148427/70‑1 (MQ=255)
                            aGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGt                                             <  1:2040273/69‑1 (MQ=255)
                             gCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTg                                                                       >  1:1728051/1‑42 (MQ=255)
                             gCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTa                                                                 >  1:2919628/1‑48 (MQ=255)
                             gCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTa                                                                 >  1:580552/1‑48 (MQ=255)
                              cTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCgg                                                   <  1:1917930/61‑1 (MQ=255)
                              cTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCgg                                                   <  1:2048920/61‑1 (MQ=255)
                                tttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGcc                                                 >  1:3873841/1‑61 (MQ=255)
                                tttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGcc                                                 >  1:1844139/1‑61 (MQ=255)
                                tttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGcc                                                 >  1:2906970/1‑61 (MQ=255)
                                 ttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGa                                                             <  1:373143/48‑1 (MQ=255)
                                 ttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGa                                                             <  1:1588279/48‑1 (MQ=255)
                                 ttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAg                                           >  1:3873781/1‑66 (MQ=255)
                                  tATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCcgcg                                              <  1:3394415/62‑1 (MQ=255)
                                    ttGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTa                                            >  1:2319156/1‑62 (MQ=255)
                                    ttGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTa                                            >  1:2486696/1‑62 (MQ=255)
                                      gaagaaGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATgg                                 >  1:3876871/1‑71 (MQ=255)
                                         gaagTGGTACCGCACGAACTGGCACATTTACTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCa                               >  1:1884457/1‑70 (MQ=255)
                                            gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAaga                           <  1:4734324/71‑1 (MQ=255)
                                            gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAaga                           <  1:197579/71‑1 (MQ=255)
                                             tGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAagag                          <  1:3780871/71‑1 (MQ=255)
                                               gTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCg                                                    >  1:2801596/1‑43 (MQ=255)
                                               gTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCg                                                    >  1:2515506/1‑43 (MQ=255)
                                               gTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCg                                                    >  1:643448/1‑43 (MQ=255)
                                                  ccGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgga                      <  1:1584355/70‑1 (MQ=255)
                                                  ccGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgga                      <  1:3034173/70‑1 (MQ=255)
                                                  ccGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgga                      <  1:584746/70‑1 (MQ=255)
                                                   cGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCaaa                             <  1:388420/62‑1 (MQ=255)
                                                   cGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCaaa                             <  1:2814216/62‑1 (MQ=255)
                                                   cGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCaaa                             <  1:1867977/62‑1 (MQ=255)
                                                        gAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgat               <  1:1273285/71‑1 (MQ=255)
                                                        gAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgat               <  1:809555/71‑1 (MQ=255)
                                                        gAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgat               >  1:3574785/1‑71 (MQ=255)
                                                         aaCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGTCAAagag                          >  1:3831423/1‑59 (MQ=255)
                                                         aaCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAagag                          >  1:1968273/1‑59 (MQ=255)
                                                         aaCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagtg                  <  1:3525448/67‑1 (MQ=255)
                                                           cTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAaga                           >  1:3264227/1‑56 (MQ=255)
                                                           cTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAaga                           >  1:3263711/1‑56 (MQ=255)
                                                           cTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgatgat            >  1:2069673/1‑71 (MQ=255)
                                                           cTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGgatgat            >  1:1895744/1‑71 (MQ=255)
                                                              gCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCcaca                                    >  1:2427704/1‑44 (MQ=255)
                                                              gCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATgg          <  1:1998718/70‑1 (MQ=255)
                                                                acaTTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGaa        >  1:2793139/1‑70 (MQ=255)
                                                                 caTTTACTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGaaaa      <  1:3811269/71‑1 (MQ=255)
                                                                 caTTTACTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGaaaa      <  1:4378543/71‑1 (MQ=255)
                                                                  aTTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagtg                  <  1:1169303/58‑1 (MQ=255)
                                                                  aTTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtggaagtg                  <  1:3637562/58‑1 (MQ=255)
                                                                  aTTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAAc     <  1:2985667/71‑1 (MQ=255)
                                                                  aTTTACTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAAc     >  1:3754060/1‑71 (MQ=255)
                                                                   tttGCTTGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACg    <  1:4412715/71‑1 (MQ=255)
                                                                   tttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGagtgg                       <  1:3970856/52‑1 (MQ=255)
                                                                   tttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAAc     <  1:2864348/70‑1 (MQ=255)
                                                                     tGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAAc     >  1:3682728/1‑68 (MQ=255)
                                                                     tGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTg  >  1:1086371/1‑71 (MQ=255)
                                                                     tGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTCGATGATGGAAAACGTg  >  1:4070986/1‑71 (MQ=255)
                                                                     tGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTAGAAGTGGATGATGGAAAAc     >  1:2907505/1‑68 (MQ=255)
                                                                      |                                                                     
ATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTG  >  minE/1976927‑1977066

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: