Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,558 T→C 9.4% 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.4% 234.8 / 12.7 85E295E (GAA→GAGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (38/39);  new base C (4/4);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGT  >  minE/1933489‑1933628
                                                                     |                                                                      
cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                       <  1:3262513/71‑1 (MQ=255)
cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                       <  1:951368/71‑1 (MQ=255)
cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                       >  1:3461071/1‑71 (MQ=255)
cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                       >  1:2614096/1‑71 (MQ=255)
 aGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAt                                                                        <  1:3573586/69‑1 (MQ=255)
  gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                     >  1:2350707/1‑71 (MQ=255)
  gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                     >  1:1254485/1‑71 (MQ=255)
  gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                     <  1:940952/71‑1 (MQ=255)
    gTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGt                                                                   <  1:3242253/71‑1 (MQ=255)
         gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                              >  1:1583329/1‑71 (MQ=255)
          tCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTc                                                             <  1:612638/71‑1 (MQ=255)
            gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                              <  1:1114222/68‑1 (MQ=255)
            gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                           <  1:1798249/71‑1 (MQ=255)
            gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                           >  1:1602227/1‑71 (MQ=255)
              cGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                          >  1:2935580/1‑70 (MQ=255)
                  ggTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                          >  1:1501163/1‑66 (MQ=255)
                   gTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                     <  1:1243590/54‑1 (MQ=255)
                    tCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:2661728/71‑1 (MQ=255)
                    tCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:1821559/71‑1 (MQ=255)
                    tCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:1711308/71‑1 (MQ=255)
                     cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                                  >  1:3284334/1‑71 (MQ=255)
                      aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCg                                                            <  1:182154/60‑1 (MQ=255)
                        ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                           <  1:1496774/59‑1 (MQ=255)
                        ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:1979196/67‑1 (MQ=255)
                        ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:528870/67‑1 (MQ=255)
                          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTc                                                                  >  1:1264942/1‑50 (MQ=255)
                          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   >  1:1112316/1‑65 (MQ=255)
                          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                             >  1:3661312/1‑71 (MQ=255)
                          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                             >  1:544646/1‑71 (MQ=255)
                          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                             >  1:1313704/1‑71 (MQ=255)
                          ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                             >  1:2637526/1‑71 (MQ=255)
                            gTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATc                                           >  1:871269/1‑71 (MQ=255)
                             tCTTACCGAGAATATCCTGTGCCAGAGTACTGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                          <  1:1710787/71‑1 (MQ=255)
                               ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                           >  1:612403/1‑52 (MQ=255)
                               ttACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                        >  1:2458388/1‑71 (MQ=255)
                                 aCCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGt                                               <  1:2297364/62‑1 (MQ=255)
                                  ccGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTa                                     <  1:2822034/71‑1 (MQ=255)
                                   cGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGAt                                                          >  1:3484164/1‑49 (MQ=255)
                                     agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                  >  1:1486025/1‑71 (MQ=255)
                                     agaATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGc                                  >  1:1038735/1‑71 (MQ=255)
                                         tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGt                                               >  1:2745682/1‑54 (MQ=255)
                                         tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAAt                                            >  1:823861/1‑57 (MQ=255)
                                         tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAg                                   >  1:2768729/1‑66 (MQ=255)
                                         tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              <  1:519088/71‑1 (MQ=255)
                                         tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              >  1:768477/1‑71 (MQ=255)
                                         tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              >  1:1074683/1‑71 (MQ=255)
                                           tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              >  1:2619424/1‑69 (MQ=255)
                                           tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              >  1:425492/1‑69 (MQ=255)
                                           tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              <  1:174480/69‑1 (MQ=255)
                                           tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                            <  1:993183/71‑1 (MQ=255)
                                           tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                            <  1:1605697/71‑1 (MQ=255)
                                           tCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTa                            <  1:3829995/71‑1 (MQ=255)
                                            ccTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGc                                        >  1:730151/1‑58 (MQ=255)
                                              tgtgCCAGAGTACGGAACAGCGATTCGGTCAGCTCg                                                            <  1:2083726/36‑1 (MQ=255)
                                               gtgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGt                        >  1:2980744/1‑71 (MQ=255)
                                               gtgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGt                        >  1:1470473/1‑71 (MQ=255)
                                                tgCCAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGtt                       <  1:1240079/71‑1 (MQ=255)
                                                tgCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGtt                       >  1:1981415/1‑71 (MQ=255)
                                                  ccAGGGTGCGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTc                      >  1:1078007/1‑70 (MQ=255)
                                                  ccAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                              <  1:786839/62‑1 (MQ=255)
                                                   cAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGAt                                                     <  1:2258319/38‑1 (MQ=255)
                                                   cAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:1820228/40‑1 (MQ=255)
                                                    agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                   >  1:1726141/1‑71 (MQ=255)
                                                    agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                   >  1:1115578/1‑71 (MQ=255)
                                                     gagTACGGAACAGCGATTCGGTCAGCTCGATCAGATca                                                   >  1:909922/1‑37 (MQ=255)
                                                      agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                   <  1:738327/37‑1 (MQ=255)
                                                      agTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                              <  1:1207152/42‑1 (MQ=255)
                                                        tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTa                                              <  1:1514177/40‑1 (MQ=255)
                                                        tACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCt                                          <  1:3159745/44‑1 (MQ=255)
                                                         gcGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcattg              >  1:1729220/2‑69 (MQ=255)
                                                          cGGAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCcatcattg              <  1:90056/70‑3 (MQ=255)
                                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTag                           <  1:3330462/57‑1 (MQ=255)
                                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg              <  1:546070/70‑1 (MQ=255)
                                                          cGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt             >  1:2260033/1‑71 (MQ=255)
                                                           ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg            <  1:1985606/71‑1 (MQ=255)
                                                           ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg            <  1:3179668/71‑1 (MQ=255)
                                                           ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg              <  1:2105788/69‑1 (MQ=255)
                                                           ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt             >  1:1987943/1‑70 (MQ=255)
                                                            gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCAt                               <  1:3466422/51‑1 (MQ=255)
                                                            gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCAt                               <  1:2437724/51‑1 (MQ=255)
                                                            gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCAt                               <  1:3551587/51‑1 (MQ=255)
                                                             aaCAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATTGTGaa          >  1:3827449/1‑71 (MQ=255)
                                                             aaCAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGaa          >  1:1244421/1‑71 (MQ=255)
                                                                aGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCg                                       <  1:3917914/39‑1 (MQ=255)
                                                                 gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct      <  1:3855310/71‑1 (MQ=255)
                                                                 gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct      <  1:2216259/71‑1 (MQ=255)
                                                                  cGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGTGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTg     <  1:1743573/71‑1 (MQ=255)
                                                                   gATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTaa                                    >  1:1856602/1‑39 (MQ=255)
                                                                   gATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActc       >  1:2322600/1‑68 (MQ=255)
                                                                   gATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTgg    >  1:485489/1‑71 (MQ=255)
                                                                     ttCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatca                 <  1:2949102/56‑1 (MQ=255)
                                                                     ttCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcatca                 <  1:2876671/56‑1 (MQ=255)
                                                                     ttCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct      <  1:1983370/67‑1 (MQ=255)
                                                                     ttCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGt  >  1:3195258/1‑71 (MQ=255)
                                                                     |                                                                      
CAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTGGGT  >  minE/1933489‑1933628

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: