Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,558 T→C 9.8% E295E (GAA→GAG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5580TC 9.8% 356.7 / 22.2 132E295E (GAA→GAGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (59/60);  new base C (5/8);  total (64/68)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 132 total aligned reads displayed.

GCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGG  >  minE/1933491‑1933626
                                                                   |                                                                    
gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:616503/1‑71 (MQ=255)
  gTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:4477600/1‑69 (MQ=255)
   tCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTc                                                                <  1:3177385/71‑1 (MQ=255)
     aCGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:861474/1‑66 (MQ=255)
     aCGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:3219102/1‑66 (MQ=255)
     aCGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:1030118/1‑66 (MQ=255)
       gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTc                                                                    <  1:4816158/63‑1 (MQ=255)
       gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTc                                                                    <  1:5471403/63‑1 (MQ=255)
       gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTc                                                                    <  1:2364566/63‑1 (MQ=255)
       gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                            >  1:3284125/1‑71 (MQ=255)
        tCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTc                                                           >  1:5632444/1‑71 (MQ=255)
         cGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                            >  1:5144250/1‑69 (MQ=255)
          gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                         >  1:5327732/1‑71 (MQ=255)
          gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                         >  1:2898334/1‑71 (MQ=255)
          gCCGTAGGTCACTTCCGTCTTACCGAGAACATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                         <  1:5013603/71‑1 (MQ=255)
            cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                               >  1:2358104/1‑63 (MQ=255)
            cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                        <  1:25966/70‑1 (MQ=255)
            cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATc                                                       <  1:546588/71‑1 (MQ=255)
             gTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                        >  1:2143972/1‑69 (MQ=255)
              tAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAg                                                     <  1:2949391/71‑1 (MQ=255)
                 gTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                               <  1:529255/58‑1 (MQ=255)
                 gTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCa                                                               <  1:4394475/58‑1 (MQ=255)
                   cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:607336/1‑52 (MQ=255)
                   cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:494790/1‑52 (MQ=255)
                   cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                   >  1:283935/1‑52 (MQ=255)
                   cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                >  1:2054902/1‑71 (MQ=255)
                    aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                               >  1:2426161/1‑71 (MQ=255)
                    aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                               >  1:1174200/1‑71 (MQ=255)
                    aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                               <  1:5454919/71‑1 (MQ=255)
                      ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                             <  1:3616657/55‑1 (MQ=255)
                      ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                             <  1:3274494/55‑1 (MQ=255)
                       tCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTc                                                                    >  1:5166528/1‑47 (MQ=255)
                        ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATc                                                       >  1:208252/1‑59 (MQ=255)
                        ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                           >  1:2175117/1‑71 (MQ=255)
                        ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                           >  1:4152859/1‑71 (MQ=255)
                        ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                           >  1:3082879/1‑71 (MQ=255)
                        ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                            <  1:4825934/70‑1 (MQ=255)
                         cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                           <  1:2970151/70‑1 (MQ=255)
                         cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                           <  1:525005/70‑1 (MQ=255)
                         cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                          <  1:4032619/71‑1 (MQ=255)
                         cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                          <  1:5629324/71‑1 (MQ=255)
                         cGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                          <  1:1125447/71‑1 (MQ=255)
                           tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGt                                                                 <  1:3726233/46‑1 (MQ=255)
                           tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGt                                                                 <  1:1097762/46‑1 (MQ=255)
                           tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                <  1:1070019/63‑1 (MQ=255)
                           tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                <  1:1569337/63‑1 (MQ=255)
                           tCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                        <  1:4385361/71‑1 (MQ=255)
                           tCTTACCGAGAATATCCTGTGACAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                               <  1:1831827/64‑1 (MQ=255)
                             ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                             <  1:851229/48‑1 (MQ=255)
                             ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGc                                                             <  1:5087448/48‑1 (MQ=255)
                             ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGt                                             >  1:363106/1‑64 (MQ=255)
                               aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                              <  1:750842/61‑1 (MQ=255)
                                   agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                >  1:251775/1‑71 (MQ=255)
                                     aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                              >  1:5144597/1‑71 (MQ=255)
                                        atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                        <  1:1886964/58‑1 (MQ=255)
                                        atCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                            >  1:4320873/1‑70 (MQ=255)
                                         tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCg                                                          <  1:4707523/39‑1 (MQ=255)
                                         tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                            >  1:2809981/1‑69 (MQ=255)
                                         tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                            >  1:3359077/1‑69 (MQ=255)
                                         tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                          <  1:4120994/71‑1 (MQ=255)
                                           cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                        <  1:4287362/39‑1 (MQ=255)
                                           cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                         <  1:4458211/70‑1 (MQ=255)
                                           cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                         <  1:2918867/70‑1 (MQ=255)
                                           cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                          >  1:2636665/1‑69 (MQ=255)
                                            tgtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTagag                       >  1:489102/1‑71 (MQ=255)
                                            tgtgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                            <  1:4200930/50‑1 (MQ=255)
                                                  agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                              >  1:4668741/1‑58 (MQ=255)
                                                  agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                              >  1:4795297/1‑58 (MQ=255)
                                                  agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                 <  1:4582410/71‑1 (MQ=255)
                                                  agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                 <  1:4214173/71‑1 (MQ=255)
                                                   gagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                 <  1:3095157/70‑1 (MQ=255)
                                                    agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatc                >  1:1944328/1‑70 (MQ=255)
                                                    agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCAGCGTa                                   >  1:2250359/1‑51 (MQ=255)
                                                     gTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                           >  1:2766433/1‑42 (MQ=255)
                                                      tACGGAACAGCGATTCGGTCAGCTCGATCAGATCttt                                               <  1:3446813/37‑1 (MQ=255)
                                                      tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                              <  1:3404465/38‑1 (MQ=255)
                                                      tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                              <  1:5620026/38‑1 (MQ=255)
                                                      tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                              <  1:4564623/38‑1 (MQ=255)
                                                      tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg             <  1:3973897/71‑1 (MQ=255)
                                                       aCGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTg                                              <  1:3490659/37‑1 (MQ=255)
                                                        cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt           <  1:4610910/71‑1 (MQ=255)
                                                        cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt           >  1:3218272/1‑71 (MQ=255)
                                                        cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt           >  1:3125452/1‑71 (MQ=255)
                                                         ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg          <  1:1453187/71‑1 (MQ=255)
                                                         ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg          <  1:1628587/71‑1 (MQ=255)
                                                          gAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg          <  1:3540865/70‑1 (MQ=255)
                                                          gAACAGTGACTCTGTCAGTTCAATCAAATCGTAGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg          <  1:207200/70‑1 (MQ=255)
                                                          gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg          >  1:4119688/1‑70 (MQ=255)
                                                           aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa        >  1:4333181/1‑71 (MQ=255)
                                                             cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct      >  1:573568/1‑71 (MQ=255)
                                                              aGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc     >  1:3809577/1‑71 (MQ=255)
                                                               gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct    >  1:3805927/1‑71 (MQ=255)
                                                               gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct    >  1:2719605/1‑71 (MQ=255)
                                                                cGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTg   >  1:4254158/1‑71 (MQ=255)
                                                                  aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                       >  1:2781024/1‑49 (MQ=255)
                                                                  aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                       >  1:2863214/1‑49 (MQ=255)
                                                                  aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                       >  1:3804555/1‑49 (MQ=255)
                                                                  aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTagag                       >  1:4052027/1‑49 (MQ=255)
                                                                  aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                         <  1:1526961/47‑1 (MQ=255)
                                                                  aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg  >  1:1752748/1‑70 (MQ=255)
                                                                   |                                                                    
GCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGG  >  minE/1933491‑1933626

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: