Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,100 T→G 8.3% D448A (GAC→GCC) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1000TG 8.3% 232.7 / 18.1 96D448A (GAC→GCC) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (50/38);  new base G (3/5);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGT  >  minE/1933031‑1933169
                                                                     |                                                                     
ggcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                                                      >  1:3454360/1‑71 (MQ=255)
ggcggTAAGCCATGTTACAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                                                      >  1:3760169/1‑71 (MQ=255)
      aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgc                                                                >  1:3583245/1‑71 (MQ=255)
       aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                                               >  1:311101/1‑71 (MQ=255)
         ccATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc                                                             <  1:2123592/71‑1 (MQ=255)
          cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                                                            <  1:3625449/71‑1 (MQ=255)
          cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                                                            <  1:152277/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                                                 <  1:994455/64‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                                                 <  1:2522804/64‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                                                 <  1:1349057/64‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCaa                                                           <  1:1067070/70‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                                                            <  1:1929935/69‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          <  1:3492382/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          <  1:363002/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          <  1:1771479/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          <  1:3954621/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          >  1:1174749/1‑71 (MQ=255)
              ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          >  1:1157642/1‑69 (MQ=255)
              ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          >  1:1757711/1‑69 (MQ=255)
                ccAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTg                                                       <  1:598391/70‑1 (MQ=255)
                  aGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                    >  1:782468/1‑71 (MQ=255)
                   gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTTGCGGCAACCTGATcc                                                   <  1:3507696/71‑1 (MQ=255)
                   gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                                                            <  1:1720469/62‑1 (MQ=255)
                   gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                                                   <  1:1430577/71‑1 (MQ=255)
                   gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                                                   <  1:470044/71‑1 (MQ=255)
                   gTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                    >  1:3010349/1‑70 (MQ=255)
                        gTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                                                            <  1:1574322/57‑1 (MQ=255)
                          gACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                    >  1:1316063/1‑63 (MQ=255)
                           aCGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCtttag                                                                 >  1:4024266/1‑47 (MQ=255)
                           aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                           <  1:515665/71‑1 (MQ=255)
                           aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                           >  1:3104385/1‑71 (MQ=255)
                           aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                           >  1:3797359/1‑71 (MQ=255)
                            cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGt                                          <  1:2941025/71‑1 (MQ=255)
                            cGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                          <  1:2302727/71‑1 (MQ=255)
                            cGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                          >  1:3679064/1‑71 (MQ=255)
                             gTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGtt                                         <  1:381791/71‑1 (MQ=255)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          <  1:2984727/54‑1 (MQ=255)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGtt                                         >  1:3037191/1‑71 (MQ=255)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                          <  1:3778205/70‑1 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGa                                                      >  1:699377/1‑56 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                    >  1:3720045/1‑58 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                    >  1:1490997/1‑58 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                       >  1:2004059/1‑71 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                       >  1:1938148/1‑71 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                       >  1:3215066/1‑71 (MQ=255)
                                aTCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGttcag                                      <  1:2250979/71‑4 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                                                        >  1:2074084/1‑53 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                          >  1:3357430/1‑67 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                      >  1:343493/1‑71 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                      >  1:2775405/1‑71 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                      >  1:1399802/1‑71 (MQ=255)
                                  cTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          >  1:1859216/1‑49 (MQ=255)
                                   ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                   >  1:1560276/1‑71 (MQ=255)
                                   ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                          <  1:1825225/48‑1 (MQ=255)
                                   ttCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAg                                                 >  1:194897/1‑57 (MQ=255)
                                    tCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc                                                             >  1:127367/1‑44 (MQ=255)
                                    tCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                   >  1:3326727/1‑70 (MQ=255)
                                    tCATCGTAGAACATCCCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                        >  1:1108282/1‑65 (MQ=255)
                                    tCATCGTAGAACATCCCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                        >  1:1112987/1‑65 (MQ=255)
                                     cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGcc                                    <  1:1098252/68‑1 (MQ=255)
                                      aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAt                                >  1:2998414/1‑71 (MQ=255)
                                       tcatAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCt                                   <  1:3612344/64‑1 (MQ=255)
                                       tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                    <  1:229421/50‑1 (MQ=255)
                                       tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                       >  1:3085872/1‑63 (MQ=255)
                                       tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                        >  1:878957/1‑62 (MQ=255)
                                       tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                               <  1:2753595/71‑1 (MQ=255)
                                        cGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGAGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCt                              >  1:2037954/1‑71 (MQ=255)
                                         atAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCtt                             <  1:1614705/70‑1 (MQ=255)
                                            gAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGAt                                >  1:507916/1‑65 (MQ=255)
                                               cATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGt                       <  1:1095778/71‑1 (MQ=255)
                                                   gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATcc                                                   <  1:2317202/39‑1 (MQ=255)
                                                   gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAtt                   >  1:133817/1‑71 (MQ=255)
                                                   gCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAt                    <  1:4038586/70‑1 (MQ=255)
                                                     ttCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTc                            >  1:511336/1‑60 (MQ=255)
                                                        gtcgtcACCTGCGTCTTTCGCGGCAACATGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCt              >  1:2566904/1‑71 (MQ=255)
                                                         tcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                           <  1:3866983/41‑1 (MQ=255)
                                                         tcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTc             >  1:1773734/1‑71 (MQ=255)
                                                          cgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGt                       <  1:1239521/60‑1 (MQ=255)
                                                           gtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAg                >  1:3465021/1‑66 (MQ=255)
                                                           gtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGc           >  1:2524913/1‑71 (MQ=255)
                                                             cACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCt                                   >  1:4053302/1‑45 (MQ=255)
                                                             cACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCt          >  1:113009/1‑70 (MQ=255)
                                                              aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTaa        <  1:144523/71‑1 (MQ=255)
                                                              aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa         <  1:3380628/70‑1 (MQ=255)
                                                              aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa         >  1:578901/1‑70 (MQ=255)
                                                                cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTg                                  <  1:1579981/43‑1 (MQ=255)
                                                                cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg      >  1:3774672/1‑71 (MQ=255)
                                                                cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg      >  1:3917514/1‑71 (MQ=255)
                                                                 tGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc     >  1:4093441/1‑71 (MQ=255)
                                                                  gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc    >  1:1709237/1‑71 (MQ=255)
                                                                  gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc    >  1:2384922/1‑71 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc    <  1:3474159/70‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc    <  1:295992/70‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc    <  1:1423393/70‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc     <  1:2446426/69‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCg   >  1:2379679/1‑71 (MQ=255)
                                                                     tCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGt  >  1:2540938/1‑70 (MQ=255)
                                                                     |                                                                     
GGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGT  >  minE/1933031‑1933169

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: