Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,555 C→T 9.0% S296S (TCG→TCA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,5550CT 9.0% 259.4 / 12.1 89S296S (TCG→TCAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (37/44);  new base T (4/4);  total (41/48)
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.80e-01

GTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTG  >  minE/1933486‑1933625
                                                                     |                                                                      
gTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCg                                                                       <  1:1305796/71‑1 (MQ=255)
 tCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGa                                                                      <  1:695720/71‑1 (MQ=255)
 tCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGa                                                                      <  1:3300165/71‑1 (MQ=255)
   cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                    <  1:951368/71‑1 (MQ=255)
   cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                    >  1:2614096/1‑71 (MQ=255)
   cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                    >  1:3461071/1‑71 (MQ=255)
   cAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAtt                                                                    <  1:3262513/71‑1 (MQ=255)
    aGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGAt                                                                     <  1:3573586/69‑1 (MQ=255)
     gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  <  1:940952/71‑1 (MQ=255)
     gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:1254485/1‑71 (MQ=255)
     gCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCg                                                                  >  1:2350707/1‑71 (MQ=255)
       gTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGt                                                                <  1:3242253/71‑1 (MQ=255)
            gTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                           >  1:1583329/1‑71 (MQ=255)
             tCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTc                                                          <  1:612638/71‑1 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGc                                                                        >  1:793209/1‑55 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCt                                                           <  1:1114222/68‑1 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                        >  1:1602227/1‑71 (MQ=255)
               gCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGa                                                        <  1:1798249/71‑1 (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:1711308/71‑1 (MQ=255)
                       tCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                <  1:2661728/71‑1 (MQ=255)
                       tCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCt                                                <  1:1821559/71‑1 (MQ=255)
                        cACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCtt                                               >  1:3284334/1‑71 (MQ=255)
                         aCTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCg                                                         <  1:182154/60‑1 (MQ=255)
                           ttCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGc                                                                        <  1:3466598/43‑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:2637526/1‑71 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          >  1:3661312/1‑71 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          >  1:544646/1‑71 (MQ=255)
                             ccGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTaa                                          >  1:1313704/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:1074683/1‑71 (MQ=255)
                                            tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                           >  1:768477/1‑71 (MQ=255)
                                            tatCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATg                           <  1:519088/71‑1 (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:1115578/1‑71 (MQ=255)
                                                       agagTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCcat                >  1:1726141/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:3179668/71‑1 (MQ=255)
                                                              ggAACAGTGACTCTGTCAGTTCAATCAAATCGTGGTAATCCGCATACGCCATGTAGAGTTCCATCATtgtg         <  1:1985606/71‑1 (MQ=255)
                                                              ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATgg           <  1:2105788/69‑1 (MQ=255)
                                                              ggAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGt          >  1:1987943/1‑70 (MQ=255)
                                                               gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCAt                            <  1:2437724/51‑1 (MQ=255)
                                                               gAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCAt                            <  1:3466422/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:2216259/71‑1 (MQ=255)
                                                                    gCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAActct   <  1:3855310/71‑1 (MQ=255)
                                                                     cGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGTGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTg  <  1:1743573/71‑1 (MQ=255)
                                                                     |                                                                      
GTCCAGCGTCACGTCGCCGTAGGTCACTTCCGTCTTACCGAGAATATCCTGTGCCAGAGTACGGAACAGCGATTCGGTCAGCTCGATCAGATCTTTGTAATCTGCGTAAGCCATGTAGAGTTCCATCATGGTGAACTCTG  >  minE/1933486‑1933625

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: