Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,949 G→A 16.2% A165V (GCG→GTG) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9490GA16.2% 170.8 / 22.2 68A165V (GCG→GTG) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base A (4/7);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGC  >  minE/1933885‑1934014
                                                                |                                                                 
cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                             <  1:1039876/71‑1 (MQ=255)
   tttAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgcgc                                                                  >  1:1055040/1‑63 (MQ=255)
     tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                             >  1:656168/1‑66 (MQ=255)
        aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                                     <  1:3104870/71‑1 (MQ=255)
        aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                                     <  1:735565/71‑1 (MQ=255)
        aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                                     <  1:3453558/71‑1 (MQ=255)
             ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                             >  1:2727446/1‑58 (MQ=255)
                  ggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                                            <  1:3509917/70‑1 (MQ=255)
                   gATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGaa                                          <  1:1842658/71‑1 (MQ=255)
                    aTTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAt                                         >  1:352920/1‑71 (MQ=255)
                     tttaTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAtt                                        >  1:395179/4‑71 (MQ=255)
                     ttCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                             >  1:3150550/1‑50 (MQ=255)
                      tCATCGTTCGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAtt                                        >  1:3280268/1‑70 (MQ=255)
                       taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                       <  1:1860825/69‑1 (MQ=255)
                       taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                       <  1:2011318/69‑1 (MQ=255)
                       cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                                       >  1:3001961/1‑70 (MQ=255)
                       cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                                       <  1:2690654/70‑1 (MQ=255)
                        aTCATTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                             >  1:2223781/1‑47 (MQ=255)
                          cGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATcc                                   >  1:571348/1‑71 (MQ=255)
                            ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           <  1:878876/61‑1 (MQ=255)
                            ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg                                 <  1:3810504/71‑1 (MQ=255)
                              ggAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGa                                                           >  1:307841/1‑43 (MQ=255)
                               gagatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATc                                                         <  1:1448303/44‑1 (MQ=255)
                                agatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGc                             >  1:2199005/1‑71 (MQ=255)
                                 gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                            <  1:2532192/71‑1 (MQ=255)
                                 gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                            <  1:616547/71‑1 (MQ=255)
                                   tgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                          >  1:3823945/1‑71 (MQ=255)
                                    gaggtCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCa                                                    >  1:365271/1‑44 (MQ=255)
                                    gaggtCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCa                                                    >  1:370217/1‑44 (MQ=255)
                                    gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                          >  1:3702366/1‑70 (MQ=255)
                                    gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                         >  1:3210298/1‑71 (MQ=255)
                                     agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTa                                      <  1:3750292/57‑1 (MQ=255)
                                     agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                         <  1:385157/70‑1 (MQ=255)
                                     agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                         <  1:1913520/70‑1 (MQ=255)
                                       atCGAGATAACGCTGACGATTGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGc                                <  1:656722/61‑1 (MQ=255)
                                       atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACg                        >  1:3037702/1‑69 (MQ=255)
                                         cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                    >  1:1552525/1‑71 (MQ=255)
                                         cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                    >  1:2960645/1‑71 (MQ=255)
                                          cagaTAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTg                   <  1:230532/70‑1 (MQ=255)
                                            gaTAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGa                                           <  1:669919/45‑1 (MQ=255)
                                             aTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAtt                                        <  1:2560908/47‑1 (MQ=255)
                                              tAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                                       <  1:2460660/47‑1 (MQ=255)
                                              tAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCtt                <  1:1059104/70‑1 (MQ=255)
                                               aaCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              >  1:796552/1‑71 (MQ=255)
                                                aCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtttag             >  1:3589246/1‑69 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGaa                                          <  1:227987/42‑1 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGaa                                          <  1:1133831/42‑1 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           >  1:3483104/1‑41 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAt                                         >  1:1691153/1‑43 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAt                                         >  1:3235777/1‑43 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTgg             >  1:314281/1‑71 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              >  1:4118009/1‑70 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              <  1:2707588/70‑1 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              <  1:3746688/70‑1 (MQ=255)
                                                aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              >  1:3009342/1‑70 (MQ=255)
                                                  gttGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTc           <  1:3419323/69‑1 (MQ=255)
                                                  gCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTc           <  1:894359/71‑1 (MQ=255)
                                                   gtGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTc           <  1:650374/69‑1 (MQ=255)
                                                   cTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTc           <  1:169050/70‑1 (MQ=255)
                                                    tGACGATAGCGCGCTTACTGATCCTGCAAGCCGTGGAATTTATGCgg                                 >  1:2703058/1‑47 (MQ=255)
                                                     gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag         <  1:2757850/70‑1 (MQ=255)
                                                      aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag         <  1:1107875/69‑1 (MQ=255)
                                                      aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag         <  1:1855939/69‑1 (MQ=255)
                                                      aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag         <  1:1797788/69‑1 (MQ=255)
                                                      aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag         <  1:89432/69‑1 (MQ=255)
                                                        gATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGc                       <  1:2411612/53‑1 (MQ=255)
                                                        gATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGa     <  1:2087235/71‑1 (MQ=255)
                                                          tAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACg   >  1:3971513/1‑71 (MQ=255)
                                                           aGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAg                              >  1:556061/1‑43 (MQ=255)
                                                            gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc  <  1:3533304/70‑1 (MQ=255)
                                                            gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc  <  1:2648389/70‑1 (MQ=255)
                                                                |                                                                 
CACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGC  >  minE/1933885‑1934014

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: