Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,398,935 C→T 5.6% E566K (GAG→AAG)  uvrD ← DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,398,9350CT 5.6% 278.2 / 3.8 89E566K (GAG→AAG) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base C (38/46);  new base T (2/3);  total (40/49)
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.95e-01

CTTCATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCACCT  >  minE/2398866‑2399005
                                                                      |                                                                      
cTTCATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTTCCATACCAACGATAAACACCTGCGGGAACTc                                                                        <  1:3245684/71‑1 (MQ=255)
 ttCATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTTCCATACCAACGATAAACACCTGCGGGAACTc                                                                        <  1:4648738/70‑1 (MQ=255)
  tCATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAg                                                                     >  1:2539370/1‑71 (MQ=255)
   cATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAgg                                                                    >  1:686058/1‑71 (MQ=255)
   cATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAgg                                                                    >  1:4335218/1‑71 (MQ=255)
    aTCCAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGc                                                                   >  1:4628640/1‑71 (MQ=255)
      ccAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCa                                                                      >  1:2613351/1‑66 (MQ=255)
      ccAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCa                                                                      >  1:4402618/1‑66 (MQ=255)
      ccAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCa                                                                      >  1:4388795/1‑66 (MQ=255)
         tcGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTc                                                              >  1:4443449/2‑71 (MQ=255)
         gCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTc                                                              <  1:229470/71‑1 (MQ=255)
             cATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCg                                                          <  1:1064174/71‑1 (MQ=255)
             cATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCg                                                          <  1:3836506/71‑1 (MQ=255)
              aTCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGa                                                         <  1:4616214/71‑1 (MQ=255)
              aTCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGa                                                         <  1:3859583/71‑1 (MQ=255)
              aTCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCAGGAACTCCAGGCCTTTCGCCg                                                          <  1:2700054/70‑1 (MQ=255)
               tCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTTCAGGCCTTTCGcc                                                           <  1:2136372/68‑1 (MQ=255)
               tCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCtt                                                                <  1:1547900/63‑1 (MQ=255)
                   gCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAgtgt                                                     <  1:3944070/70‑1 (MQ=255)
                   gCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAg                                                        <  1:2289197/67‑1 (MQ=255)
                   gCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTa                                                    <  1:3431095/71‑1 (MQ=255)
                       gggAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCg                                                 <  1:3985704/70‑1 (MQ=255)
                        ggAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGt                                                <  1:2699141/70‑1 (MQ=255)
                           aCATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:3559474/71‑1 (MQ=255)
                           aCATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:4496791/71‑1 (MQ=255)
                           aCATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:2440755/71‑1 (MQ=255)
                           aCATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcat                                             >  1:1417298/1‑70 (MQ=255)
                           aCATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcat                                             >  1:1344609/1‑70 (MQ=255)
                           aCATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcat                                             >  1:4498732/1‑70 (MQ=255)
                                cccTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTc                                                              >  1:997116/1‑48 (MQ=255)
                                cccTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTg                                       <  1:4691528/71‑1 (MQ=255)
                                  ctctTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCa                                     >  1:3995093/1‑71 (MQ=255)
                                   tctTCCATACCAACGATAAACACCTGCGGGAACTTCAg                                                                     <  1:3237463/38‑1 (MQ=255)
                                   tctTCCATACCAACGATAAACACCTGCGGGAACTTCAg                                                                     <  1:3782865/38‑1 (MQ=255)
                                    ctTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCtt                                                                >  1:805704/1‑42 (MQ=255)
                                    ctTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTc                                                              >  1:2932281/1‑44 (MQ=255)
                                      ttCCATACCAACGATAAACACCTGCGGGAACTCCa                                                                      >  1:3451677/1‑35 (MQ=255)
                                      ttCCATACCAACGATAAACACCTGCGGGAACTCCAGGCGTTTCGCCGAGTGTAGCGTCATCAACTGCACCg                                  >  1:1133556/1‑71 (MQ=255)
                                      ttCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAAc                                         >  1:2893230/1‑64 (MQ=255)
                                      ttCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCg                                  >  1:1189787/1‑71 (MQ=255)
                                        ccATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCAcc                                   <  1:639514/68‑1 (MQ=255)
                                        ccATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCa                                >  1:1737541/1‑71 (MQ=255)
                                          aTACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATc                              <  1:2278140/71‑1 (MQ=255)
                                          aTACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATc                              <  1:2063099/71‑1 (MQ=255)
                                           tACCAACGATAAACACCTGCGGGAACTCCAGGCCtt                                                                <  1:2495751/36‑1 (MQ=255)
                                           tACCAACGATAAACACCTGCGGGAACTCCAGGCCtt                                                                <  1:551660/36‑1 (MQ=255)
                                           tACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtc                                               <  1:211152/53‑1 (MQ=255)
                                           tACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCa                                <  1:1515800/68‑1 (MQ=255)
                                           tACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCa                                <  1:4661710/68‑1 (MQ=255)
                                           tACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATcc                             >  1:681333/1‑71 (MQ=255)
                                            aCCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCt                            >  1:3952628/1‑71 (MQ=255)
                                             ccAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGa                                                         >  1:854351/1‑41 (MQ=255)
                                             ccAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCg                                                 <  1:1942177/49‑1 (MQ=255)
                                               aaCGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatca                                           <  1:4050397/53‑1 (MQ=255)
                                                aCGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtca                                              >  1:325341/1‑49 (MQ=255)
                                                aCGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCa                                <  1:17413/63‑1 (MQ=255)
                                                aCGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGc                          >  1:1217024/1‑69 (MQ=255)
                                                 cGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:530259/50‑1 (MQ=255)
                                                  gATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:625458/49‑1 (MQ=255)
                                                  gATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAgg                      >  1:2611243/1‑71 (MQ=255)
                                                   aTAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAgg                      >  1:1883819/1‑70 (MQ=255)
                                                    tAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:1390460/47‑1 (MQ=255)
                                                    tAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:1359118/47‑1 (MQ=255)
                                                    tAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTa                    <  1:104661/71‑1 (MQ=255)
                                                     aaaCACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCaa                                          >  1:1328215/1‑48 (MQ=255)
                                                     aaaCACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTAt                   >  1:2830937/1‑71 (MQ=255)
                                                     aaaCACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTAt                   >  1:1561750/1‑71 (MQ=255)
                                                      aaCACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAAc                                         >  1:4519551/1‑48 (MQ=255)
                                                      aaCACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGc                                      <  1:4290243/51‑1 (MQ=255)
                                                      aaCACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCa                                <  1:205685/57‑1 (MQ=255)
                                                       acacCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGtcatc                                            <  1:283924/44‑1 (MQ=255)
                                                       acacCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATcc                 <  1:3669156/71‑1 (MQ=255)
                                                        cacCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCg                >  1:877681/1‑71 (MQ=255)
                                                         acCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGc               <  1:435655/71‑1 (MQ=255)
                                                           cTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCt             >  1:3940927/1‑71 (MQ=255)
                                                             gCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATc                              >  1:1501149/1‑52 (MQ=255)
                                                             gCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGc           >  1:2081446/1‑71 (MQ=255)
                                                                ggAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCg                                  <  1:1633914/45‑1 (MQ=255)
                                                                ggAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCt        <  1:1871309/71‑1 (MQ=255)
                                                                  aaCTTCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCtt       >  1:2999698/1‑70 (MQ=255)
                                                                  aaCTTCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCtt       >  1:3406293/1‑70 (MQ=255)
                                                                  aaCTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTc      <  1:1529144/71‑1 (MQ=255)
                                                                  aaCTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTc      <  1:3048908/71‑1 (MQ=255)
                                                                  aaCTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTAACCGCCTGCCCTTc      <  1:68940/71‑1 (MQ=255)
                                                                   aCTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTc      <  1:1899822/70‑1 (MQ=255)
                                                                   aCTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCa     >  1:1923991/1‑71 (MQ=255)
                                                                    cTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTa                    <  1:1206046/55‑1 (MQ=255)
                                                                    cTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATcc                 <  1:2569599/58‑1 (MQ=255)
                                                                     tCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCtt       <  1:3343603/67‑1 (MQ=255)
                                                                     tCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCAcc   >  1:3667140/1‑71 (MQ=255)
                                                                     tCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCAcc   >  1:1621652/1‑71 (MQ=255)
                                                                      ccAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATc                  >  1:323546/1‑55 (MQ=255)
                                                                      ccAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGcc              >  1:203613/1‑59 (MQ=255)
                                                                      ccAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCACCt  >  1:2161224/1‑71 (MQ=255)
                                                                      ccAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCACCt  >  1:1978110/1‑71 (MQ=255)
                                                                      ccAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCACCt  >  1:4563394/1‑71 (MQ=255)
                                                                      ccAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCACCt  >  1:4318078/1‑71 (MQ=255)
                                                                      |                                                                      
CTTCATCCAGCGACATCTGGCTTGGGAACATGCCCTCTTCCATACCAACGATAAACACCTGCGGGAACTCCAGGCCTTTCGCCGAGTGTAGCGTCATCAACTGCACCGCATCCTGCCAGGTATCCGCCTGCCCTTCACCT  >  minE/2398866‑2399005

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: