Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,981 C→T 10.9% P154P (CCG→CCA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9810CT10.9% 214.2 / 13.7 73P154P (CCG→CCAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (38/27);  new base T (4/4);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

CGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTT  >  minE/1933911‑1934047
                                                                      |                                                                  
cGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATcc                                                                    >  1:571348/1‑71 (MQ=255)
  ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg                                                                  <  1:3810504/71‑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)
          gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                                                           >  1:3702366/1‑70 (MQ=255)
          gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                                                          >  1:3210298/1‑71 (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:2960645/1‑71 (MQ=255)
               cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                     >  1:1552525/1‑71 (MQ=255)
                cagaTAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTg                                                    <  1:230532/70‑1 (MQ=255)
                    tAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCtt                                                 <  1:1059104/70‑1 (MQ=255)
                     aaCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                               >  1:796552/1‑71 (MQ=255)
                      aCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtttag                                              >  1:3589246/1‑69 (MQ=255)
                      aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTgg                                              >  1:314281/1‑71 (MQ=255)
                      aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                               <  1:3746688/70‑1 (MQ=255)
                      aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                               <  1:2707588/70‑1 (MQ=255)
                      aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                               >  1:4118009/1‑70 (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:1855939/69‑1 (MQ=255)
                            aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                                          <  1:1107875/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)
                                         tCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCa                                                       >  1:3720702/1‑43 (MQ=255)
                                         tCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCa                                                       >  1:1409043/1‑43 (MQ=255)
                                          ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTca                             >  1:3455490/1‑67 (MQ=255)
                                          ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                          >  1:2914750/1‑70 (MQ=255)
                                          ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                          >  1:1296557/1‑70 (MQ=255)
                                          ccTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTg                          >  1:1651608/1‑71 (MQ=255)
                                            tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTg                          >  1:1354753/1‑69 (MQ=255)
                                                ccTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGCCGCAACTCGGTGCAGTg                     >  1:427977/1‑70 (MQ=255)
                                                ccTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                     >  1:2293852/1‑70 (MQ=255)
                                                ccTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTCCAGTg                     >  1:2777200/1‑70 (MQ=255)
                                                   gCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGt                           >  1:914479/1‑61 (MQ=255)
                                                   gCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                     <  1:101635/67‑1 (MQ=255)
                                                   gCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                   >  1:1457159/1‑69 (MQ=255)
                                                    cAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGa                                      >  1:1398104/1‑49 (MQ=255)
                                                    cAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                   >  1:1697902/1‑68 (MQ=255)
                                                       gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCaa                                 >  1:2685907/1‑51 (MQ=255)
                                                       gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCaa                                 >  1:1743731/1‑51 (MQ=255)
                                                       gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa             <  1:530084/71‑1 (MQ=255)
                                                           tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                    <  1:3344496/60‑1 (MQ=255)
                                                           tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTc         >  1:2591501/1‑71 (MQ=255)
                                                           tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTc         <  1:3086584/71‑1 (MQ=255)
                                                            ggAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAAc                                >  1:162055/1‑47 (MQ=255)
                                                             gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa               <  1:1069585/63‑1 (MQ=255)
                                                             gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa               >  1:2257452/1‑63 (MQ=255)
                                                             gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCg        >  1:460989/1‑70 (MQ=255)
                                                               aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc                                   <  1:4027176/41‑1 (MQ=255)
                                                               aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                   >  1:2629058/1‑57 (MQ=255)
                                                               aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                   >  1:1534944/1‑57 (MQ=255)
                                                                tttATCCGGCAGCGGACGCAGTGCTTTGGTcagca                                        <  1:2718200/35‑1 (MQ=255)
                                                                tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGc                         <  1:2526827/50‑1 (MQ=255)
                                                                 ttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCaa                                 >  1:579549/1‑41 (MQ=255)
                                                                  tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGa                                      >  1:1983854/1‑35 (MQ=255)
                                                                  tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGa                                      >  1:3563845/1‑35 (MQ=255)
                                                                  tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtt  >  1:1964796/1‑71 (MQ=255)
                                                                  tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtt  >  1:3514143/1‑71 (MQ=255)
                                                                      |                                                                  
CGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTT  >  minE/1933911‑1934047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: