Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 415,960 A→G 6.0% I253I (ATT→ATC ybeZ ← hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE415,9600AG 6.0% 301.9 / 8.2 100I253I (ATT→ATCybeZhypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (45/49);  new base G (4/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 111 total aligned reads displayed.

TGAGTTAAAACCGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTT  >  minE/415891‑416025
                                                                     |                                                                 
tGAGTTAAAACCGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCgaga                                                                   <  1:202153/70‑1 (MQ=255)
tGAGTTAAAACCGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGaa                                                                  >  1:1530854/1‑71 (MQ=255)
     tAAAACCGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGAt                                                              >  1:4357761/1‑70 (MQ=255)
      aaaaCCGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATa                                                             >  1:2623437/1‑70 (MQ=255)
          ccGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAAcgc                                                        >  1:2976042/1‑71 (MQ=255)
            gATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGcgt                                                      <  1:2846761/71‑1 (MQ=255)
                cGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGcgtcgt                                                   <  1:4644514/70‑1 (MQ=255)
                  ggTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCa                                                <  1:179884/71‑1 (MQ=255)
                      agtaaCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGt                                            <  1:2819391/68‑1 (MQ=255)
                      aGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGt                                            <  1:3031544/71‑1 (MQ=255)
                      aGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGt                                            <  1:3891284/71‑1 (MQ=255)
                      aGGAACATCTTCATCTGCTCGACGGTCGTGTTCTGGGCTTCGTCCAGGATCACGAACGCGTGGTTc                                                 >  1:3525706/1‑66 (MQ=255)
                      aGGAACATCTTCATCTGCTCGACGGTCGTGTTCTGGGCTTCGTCCAGGATCACGAACGCGTGGTTCAgg                                              >  1:1431861/1‑68 (MQ=255)
                       ggAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGt                                            <  1:370660/70‑1 (MQ=255)
                        gAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTa                                           >  1:2739599/1‑70 (MQ=255)
                        gAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTa                                           >  1:3005192/1‑70 (MQ=255)
                        gAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTa                                           >  1:4744366/1‑70 (MQ=255)
                        gAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTAc                                          >  1:2057428/1‑71 (MQ=255)
                        gAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTAc                                          >  1:3465459/1‑71 (MQ=255)
                         aaCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTAc                                          >  1:4577611/1‑70 (MQ=255)
                         aaCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACg                                         <  1:3546143/71‑1 (MQ=255)
                          aCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGaa                                                                  >  1:2467472/1‑45 (MQ=255)
                          aCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGaa                                                                  >  1:2556841/1‑45 (MQ=255)
                          aCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGaa                                                                  >  1:2735683/1‑45 (MQ=255)
                          aCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGaa                                                                  >  1:2204159/1‑45 (MQ=255)
                          aCATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACg                                         <  1:1101927/70‑1 (MQ=255)
                           cATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATa                                                             <  1:4088829/49‑1 (MQ=255)
                           cATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATa                                                             <  1:4115012/49‑1 (MQ=255)
                            aTCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCTTCGAGAATGATa                                                             <  1:3324702/48‑1 (MQ=255)
                               ttCATCTGTTCGATGGTGGTGTTCTGGCTCTCTTCGAGAATGATAAACGCGTCGTTCAGCGTAc                                          >  1:1895402/1‑64 (MQ=255)
                               ttCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATaa                                                            >  1:3380947/1‑46 (MQ=255)
                               ttCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATaa                                                            >  1:2655258/1‑46 (MQ=255)
                               ttCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTa                                           >  1:2112583/1‑63 (MQ=255)
                               ttCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTAc                                          >  1:3156393/1‑64 (MQ=255)
                               ttCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTAc                                          >  1:2840705/1‑64 (MQ=255)
                                 cATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCAc                                    <  1:531802/68‑1 (MQ=255)
                                  aTCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCa                                 <  1:3825442/70‑1 (MQ=255)
                                   tCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:1646717/1‑71 (MQ=255)
                                     tGTTCGATGGTGGTGTTCTGGCTCTCATCGAGGATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:1186726/1‑69 (MQ=255)
                                      gTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:302678/1‑68 (MQ=255)
                                       ttCGATGGTTGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               <  1:165374/67‑1 (MQ=255)
                                       ttCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACg                                         <  1:3809886/57‑1 (MQ=255)
                                       ttCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               <  1:2727805/67‑1 (MQ=255)
                                       ttCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               <  1:4196811/67‑1 (MQ=255)
                                       ttCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               <  1:4695784/67‑1 (MQ=255)
                                       atCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAACCGCGTCGTTCAGCGTACGACCACGCata                               <  1:3172875/66‑1 (MQ=255)
                                         cGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGCCCACGCATATAGGc                          >  1:588577/1‑70 (MQ=255)
                                         cGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGc                          >  1:3962830/1‑70 (MQ=255)
                                         cGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGc                          >  1:1185104/1‑70 (MQ=255)
                                          gATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGcc                         <  1:1113756/70‑1 (MQ=255)
                                          gATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGcc                         <  1:2045746/70‑1 (MQ=255)
                                           aTGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGcc                         <  1:4052627/69‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:1696876/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:3428267/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:2253320/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:2613602/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:4197899/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:2255534/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      <  1:4750944/71‑1 (MQ=255)
                                            tggtggTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGc                      >  1:1865383/1‑71 (MQ=255)
                                            cggtcgtGTTCTGGGCTTCGTCCAGGATCACGAACGCGTGGTTCAGGGTGCGCCCGCGCATATAGGCCAGc                      <  1:654269/70‑1 (MQ=255)
                                            cggtcgtGTTCTGGGCTTCGTCCAGGATCACGAACGCGTGGTTCAGGGTGCGCCCGCGCATATAGGCCAGc                      <  1:474649/70‑1 (MQ=255)
                                               tggtgTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGt                   <  1:840300/71‑1 (MQ=255)
                                               tggtgTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGt                   <  1:2768528/71‑1 (MQ=255)
                                               gtgtgtTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCAc                                    <  1:1532803/52‑1 (MQ=255)
                                                ggtgTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGt                   >  1:3322729/1‑70 (MQ=255)
                                                ggtgTTCTGGCTCTAATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGt                   >  1:557101/1‑70 (MQ=255)
                                                 gtgtTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGc                 >  1:849650/1‑71 (MQ=255)
                                                 gtgtGCTGGCTCTCATCGAGAATGATAAACGCGGCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGc                 <  1:3044325/71‑1 (MQ=255)
                                                   gtTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAg                            >  1:4430137/1‑58 (MQ=255)
                                                      cTGGCTCTCATCGAGAATGATAAACGCGTCGTTCa                                                >  1:3556052/1‑35 (MQ=255)
                                                      cTGGCTCTCATCGAGAATGATAAACGCGTCGTTCa                                                >  1:3487815/1‑35 (MQ=255)
                                                      cTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:842776/1‑52 (MQ=255)
                                                      cTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:206683/1‑52 (MQ=255)
                                                      cTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:4221569/1‑52 (MQ=255)
                                                      cTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACtt            >  1:760054/1‑71 (MQ=255)
                                                          ctctCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCat                                >  1:2979153/1‑47 (MQ=255)
                                                          ctctCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCaa         >  1:2452422/1‑70 (MQ=255)
                                                           tctcATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAAt        >  1:2887041/1‑70 (MQ=255)
                                                             tcATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAAc     >  1:4115116/1‑71 (MQ=255)
                                                             tcATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAAc     >  1:320400/1‑71 (MQ=255)
                                                              cATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCAtata                             >  1:2515331/1‑46 (MQ=255)
                                                              cATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCAtata                             >  1:3731579/1‑46 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCata                               >  1:2784660/1‑43 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATaa      <  1:2098809/68‑1 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATaa      <  1:3880971/68‑1 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAAc     <  1:2729676/69‑1 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAAc     <  1:1161111/69‑1 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACg    <  1:1342408/70‑1 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGt   >  1:560368/1‑71 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGt   >  1:3830188/1‑71 (MQ=255)
                                                               aTCGAGAATGATAAACGCGTCGTTCAGCGTACGACAACGCATATAGGCCAGCGGTGCGACTTCAATAACGt   >  1:3218378/1‑71 (MQ=255)
                                                                tCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCa                                 <  1:1090200/40‑1 (MQ=255)
                                                                tCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACt             <  1:638351/60‑1 (MQ=255)
                                                                tCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACt             <  1:1515911/60‑1 (MQ=255)
                                                                tCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGtt  <  1:1814320/71‑1 (MQ=255)
                                                                  gagaATGATAAACGCGTCGTTCAGCGTACGACCACGCa                                 <  1:2504673/38‑1 (MQ=255)
                                                                  gagaATGATAAACGCGTCGTTCAGCGTACGACCACGCa                                 <  1:1378749/38‑1 (MQ=255)
                                                                   agaATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTg                  <  1:2459290/52‑1 (MQ=255)
                                                                   agaATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTg                  <  1:4286083/52‑1 (MQ=255)
                                                                     |                                                                 
TGAGTTAAAACCGATACGGGTCAGGAACATCTTCATCTGTTCGATGGTGGTGTTCTGGCTCTCATCGAGAATGATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTT  >  minE/415891‑416025

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: