Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,993 C→T 12.1% L150L (CTG→CTA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9930CT12.1% 186.9 / 14.0 66L150L (CTG→CTAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base T (4/4);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGC  >  minE/1933926‑1934062
                                                                   |                                                                     
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:2707588/70‑1 (MQ=255)
       aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                                              >  1:3009342/1‑70 (MQ=255)
       aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                                              >  1:4118009/1‑70 (MQ=255)
       aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTg                                                              <  1:3746688/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)
            gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                         <  1:2757850/70‑1 (MQ=255)
             aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                                                         <  1:1107875/69‑1 (MQ=255)
             aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                                                         <  1:89432/69‑1 (MQ=255)
             aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                                                         <  1:1797788/69‑1 (MQ=255)
             aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                                                         <  1:1855939/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)
                   gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc                                                  <  1:2648389/70‑1 (MQ=255)
                   gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc                                                  <  1:3533304/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:1296557/1‑70 (MQ=255)
                           ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                         >  1:2914750/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:1743731/1‑51 (MQ=255)
                                        gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCaa                                                >  1:2685907/1‑51 (MQ=255)
                                        gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa                            <  1:530084/71‑1 (MQ=255)
                                            tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                                   <  1:3344496/60‑1 (MQ=255)
                                            tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTc                        <  1:3086584/71‑1 (MQ=255)
                                            tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTc                        >  1:2591501/1‑71 (MQ=255)
                                             ggAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAAc                                               >  1:162055/1‑47 (MQ=255)
                                              gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                              >  1:2257452/1‑63 (MQ=255)
                                              gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                              <  1:1069585/63‑1 (MQ=255)
                                              gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCg                       >  1:460989/1‑70 (MQ=255)
                                                aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc                                                  <  1:4027176/41‑1 (MQ=255)
                                                aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                                  >  1:1534944/1‑57 (MQ=255)
                                                aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                                  >  1:2629058/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:3563845/1‑35 (MQ=255)
                                                   tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGa                                                     >  1:1983854/1‑35 (MQ=255)
                                                   tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtt                 >  1:3514143/1‑71 (MQ=255)
                                                   tATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtt                 >  1:1964796/1‑71 (MQ=255)
                                                         gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGAGCAGTGGATAGACAg                           >  1:3058997/1‑55 (MQ=255)
                                                           aGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggttt          <  1:3387530/70‑1 (MQ=255)
                                                            gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg        <  1:1079027/71‑1 (MQ=255)
                                                            gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg        <  1:1151958/71‑1 (MQ=255)
                                                              ggACGCAGTGCTTTGGTCAGCAGACGCAACTCGGt                                          >  1:222158/1‑35 (MQ=255)
                                                                aCGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACa    <  1:1868732/71‑1 (MQ=255)
                                                                  gCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGc  >  1:1361091/1‑71 (MQ=255)
                                                                   cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                                    <  1:1479287/36‑1 (MQ=255)
                                                                   |                                                                     
CGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGC  >  minE/1933926‑1934062

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: