Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,969 G→A 15.8% H158H (CAC→CAT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9690GA15.8% 194.8 / 24.6 76H158H (CAC→CATlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base A (5/7);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACA  >  minE/1933903‑1934036
                                                                  |                                                                   
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)
    tCATCGTTCGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAtt                                                              >  1:3280268/1‑70 (MQ=255)
     taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                                             <  1:2011318/69‑1 (MQ=255)
     taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                                             <  1:1860825/69‑1 (MQ=255)
     cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                                                             <  1:2690654/70‑1 (MQ=255)
     cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                                                             >  1:3001961/1‑70 (MQ=255)
        cGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATcc                                                         >  1:571348/1‑71 (MQ=255)
          ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                                                 <  1:878876/61‑1 (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)
                   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:2960645/1‑71 (MQ=255)
                       cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                          >  1:1552525/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:1133831/42‑1 (MQ=255)
                              aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGaa                                                                <  1:227987/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:1855939/69‑1 (MQ=255)
                                    aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                               <  1:89432/69‑1 (MQ=255)
                                    aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcag                               <  1:1107875/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: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: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:1743731/1‑51 (MQ=255)
                                                               gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCaa                      >  1:2685907/1‑51 (MQ=255)
                                                               gCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa  <  1:530084/71‑1 (MQ=255)
                                                                  |                                                                   
GGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACA  >  minE/1933903‑1934036

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: