Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,594 A→C 23.0% A283A (GCT→GCG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5940AC23.0% 161.6 / 72.5 113A283A (GCT→GCGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (47/40);  new base C (13/13);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 125 total aligned reads displayed.

AGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAC  >  minE/1933526‑1933664
                                                                    |                                                                      
agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                                                      >  1:251775/1‑71 (MQ=255)
  aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                                    >  1:5144597/1‑71 (MQ=255)
  aaTATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                                    <  1:3865103/71‑1 (MQ=255)
      tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                                  >  1:2809981/1‑69 (MQ=255)
      tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                                                                  >  1:3359077/1‑69 (MQ=255)
      tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                                <  1:4120994/71‑1 (MQ=255)
        cTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                                               <  1:1164551/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)
          gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                                            >  1:2978081/1‑71 (MQ=255)
          gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                                                            >  1:2505897/1‑71 (MQ=255)
             ccAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                                                          <  1:2002871/70‑1 (MQ=255)
               agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                                    >  1:4795297/1‑58 (MQ=255)
               agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCa                                                                    >  1:4668741/1‑58 (MQ=255)
               agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                       <  1:4214173/71‑1 (MQ=255)
               agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                       <  1:4582410/71‑1 (MQ=255)
                gagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                                                       <  1:3095157/70‑1 (MQ=255)
                   tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATg                                                   <  1:3973897/71‑1 (MQ=255)
                     cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt                                                 <  1:4610910/71‑1 (MQ=255)
                     cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGt                                                 <  1:5607190/71‑1 (MQ=255)
                     cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                                                 >  1:3125452/1‑71 (MQ=255)
                     cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                                                 >  1:3174490/1‑71 (MQ=255)
                     cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                                                 >  1:3218272/1‑71 (MQ=255)
                      ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                                                <  1:2298891/71‑1 (MQ=255)
                      ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                                                <  1:1453187/71‑1 (MQ=255)
                      ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                                <  1:4345748/55‑1 (MQ=255)
                      ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                                <  1:1628587/71‑1 (MQ=255)
                       gAACAGTGACTCTGTCAGTTCAATCAAATCGTAGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg                                                <  1:207200/70‑1 (MQ=255)
                       gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                                >  1:4119688/1‑70 (MQ=255)
                        aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa                                              >  1:4333181/1‑71 (MQ=255)
                        aaCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGaa                                              <  1:655701/71‑1 (MQ=255)
                         aCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt                                                 <  1:4098631/67‑1 (MQ=255)
                          cAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAAct                                            <  1:4462152/71‑1 (MQ=255)
                          cAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                                            >  1:573568/1‑71 (MQ=255)
                           aGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc                                           >  1:3809577/1‑71 (MQ=255)
                            gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                                          >  1:5368644/1‑71 (MQ=255)
                            gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc                                           >  1:4926756/1‑70 (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)
                               aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                                                               <  1:1526961/47‑1 (MQ=255)
                               aTTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg                                        >  1:1752748/1‑70 (MQ=255)
                                  cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATTGTGAAct                                            <  1:4907958/63‑1 (MQ=255)
                                  cGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAAct                                            <  1:2997913/63‑1 (MQ=255)
                                   tgTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTAt                                   >  1:4743265/2‑71 (MQ=255)
                                    gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                                  <  1:658033/71‑1 (MQ=255)
                                    gTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATg                                  <  1:3003014/71‑1 (MQ=255)
                                     tCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                                                          >  1:1794025/1‑46 (MQ=255)
                                     tCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcatt                                                   <  1:4831082/53‑2 (MQ=255)
                                     tCAGTTCAATCAAATAGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTa                                    <  1:3702415/68‑1 (MQ=255)
                                       aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg                               >  1:1631772/1‑71 (MQ=255)
                                       aGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg                               >  1:2144389/1‑71 (MQ=255)
                                          tCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                                                                >  1:168352/1‑35 (MQ=255)
                                           cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                                          <  1:249314/56‑1 (MQ=255)
                                           cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                                          <  1:3499983/56‑1 (MQ=255)
                                           cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct                                          <  1:290926/56‑1 (MQ=255)
                                           cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTggg                                       >  1:4682924/1‑59 (MQ=255)
                                           cGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg                           >  1:5580556/1‑71 (MQ=255)
                                            gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc                            <  1:2618209/69‑1 (MQ=255)
                                            gATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTAc                            <  1:4813837/69‑1 (MQ=255)
                                             aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                                  <  1:3696514/46‑1 (MQ=255)
                                             aTCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                                  <  1:3918940/46‑1 (MQ=255)
                                              taaGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGa                         <  1:2749264/68‑1 (MQ=255)
                                              tCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcat                                                    <  1:3926212/43‑1 (MQ=255)
                                                aGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGaaa                       >  1:4235832/1‑70 (MQ=255)
                                                    cTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACg                           >  1:3107546/1‑62 (MQ=255)
                                                     tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg                               <  1:5237807/57‑1 (MQ=255)
                                                     tttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACg                               <  1:1190534/57‑1 (MQ=255)
                                                     gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCtt                 >  1:4159908/4‑71 (MQ=255)
                                                     gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCtt                 >  1:5217119/4‑71 (MQ=255)
                                                     gtggTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCtt                 >  1:2696559/4‑71 (MQ=255)
                                                      ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa                                    >  1:320873/1‑51 (MQ=255)
                                                      ttGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTa                                    >  1:110065/1‑51 (MQ=255)
                                                       tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                                <  1:70909/38‑1 (MQ=255)
                                                       tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTg                                                <  1:1784391/38‑1 (MQ=255)
                                                       tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTa                             >  1:4402170/1‑57 (MQ=255)
                                                       tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCt                  <  1:4508666/68‑1 (MQ=255)
                                                       tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTc                <  1:281655/70‑1 (MQ=255)
                                                       tGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCg               >  1:4408375/1‑71 (MQ=255)
                                                        gTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                                  >  1:4041909/1‑35 (MQ=255)
                                                        gTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg                                                  >  1:1801738/1‑35 (MQ=255)
                                                         tAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGt              <  1:3745819/70‑1 (MQ=255)
                                                         tAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTca               <  1:4238562/69‑2 (MQ=255)
                                                         tAATCCGCATACGCCATGTAGAGTTCCATCATTGTGAACTCAGGATTATGGCGAACAGAAATACCTTca               <  1:2503967/69‑2 (MQ=255)
                                                           aTCCGCATACGCCATGTAGAGTTACATCATTGTGAACTCAGGATTATGGCGAACAGaa                        >  1:4041847/1‑58 (MQ=255)
                                                              tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg         >  1:3010245/1‑70 (MQ=255)
                                                              tGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACgg         >  1:46491/1‑70 (MQ=255)
                                                                cGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGa        <  1:2944428/69‑1 (MQ=255)
                                                                cGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAg      <  1:5230058/71‑1 (MQ=255)
                                                                   aaGCCATGTAGATTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa   >  1:3969084/1‑71 (MQ=255)
                                                                   aaGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGt     <  1:5298862/69‑1 (MQ=255)
                                                                   aaGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa   >  1:2642069/1‑71 (MQ=255)
                                                                    aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATa                     >  1:1549854/1‑52 (MQ=255)
                                                                    aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTa   <  1:125228/70‑1 (MQ=255)
                                                                    aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc  <  1:489138/71‑1 (MQ=255)
                                                                    aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc  >  1:4836025/1‑71 (MQ=255)
                                                                    aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc  <  1:166163/71‑1 (MQ=255)
                                                                    aGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAc  <  1:1531389/71‑1 (MQ=255)
                                                                    |                                                                      
AGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTAC  >  minE/1933526‑1933664

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: