Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,948 C→G 16.1% A165A (GCG→GCC) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9480CG16.1% 159.9 / 21.2 68A165A (GCG→GCC) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base G (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 = 9.81e-01

GCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGC  >  minE/1933878‑1934014
                                                                      |                                                                  
gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgc                                                                    >  1:4098853/1‑71 (MQ=255)
gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgc                                                                    >  1:2401052/1‑71 (MQ=255)
       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:2690654/70‑1 (MQ=255)
                              cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                                       >  1:3001961/1‑70 (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:3746688/70‑1 (MQ=255)
                                                       aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              >  1:3009342/1‑70 (MQ=255)
                                                       aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg              <  1:2707588/70‑1 (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:89432/69‑1 (MQ=255)
                                                             aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag         <  1:1797788/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)
                                                                      |                                                                  
GCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGC  >  minE/1933878‑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: