Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,099 G→A 8.3% D448D (GAC→GAT) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,0990GA 8.3% 300.6 / 12.6 96D448D (GAC→GAT) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (50/38);  new base A (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 = 9.70e-01

GGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCG  >  minE/1933031‑1933168
                                                                    |                                                                     
ggcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                                                     >  1:3454360/1‑71 (MQ=255)
ggcggTAAGCCATGTTACAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                                                     >  1:3760169/1‑71 (MQ=255)
    gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCg                                                                       <  1:3195856/65‑1 (MQ=255)
    gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCg                                                                       <  1:1282271/65‑1 (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:152277/71‑1 (MQ=255)
          cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa                                                           <  1:3625449/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                                                <  1:1349057/64‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                                                <  1:2522804/64‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg                                                                <  1:994455/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:1174749/1‑71 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                         <  1:3954621/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                         <  1:1771479/71‑1 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                         <  1:363002/71‑1 (MQ=255)
              ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                         >  1:1757711/1‑69 (MQ=255)
              ttCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                         >  1:1157642/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:3797359/1‑71 (MQ=255)
                           aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                          <  1:515665/71‑1 (MQ=255)
                           aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                          >  1:3104385/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:3215066/1‑71 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                      >  1:2004059/1‑71 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                      >  1:1938148/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:578901/1‑70 (MQ=255)
                                                              aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa        <  1:3380628/70‑1 (MQ=255)
                                                                cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTg                                 <  1:1579981/43‑1 (MQ=255)
                                                                cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg     >  1:3917514/1‑71 (MQ=255)
                                                                cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg     >  1:3774672/1‑71 (MQ=255)
                                                                 tGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc    >  1:4093441/1‑71 (MQ=255)
                                                                  gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc   >  1:2384922/1‑71 (MQ=255)
                                                                  gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc   >  1:1709237/1‑71 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc   <  1:1423393/70‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc   <  1:3474159/70‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc   <  1:295992/70‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc    <  1:2446426/69‑1 (MQ=255)
                                                                   cGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCg  >  1:2379679/1‑71 (MQ=255)
                                                                    |                                                                     
GGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCG  >  minE/1933031‑1933168

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: