Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,936 C→T 11.9% Q169Q (CAG→CAA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9360CT11.9% 189.8 / 13.0 67Q169Q (CAG→CAAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (4/4);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTC  >  minE/1933872‑1934005
                                                                |                                                                     
ggATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGa                                                                 <  1:3847390/71‑1 (MQ=255)
     tGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGa                                                                    <  1:719903/63‑1 (MQ=255)
      gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgc                                                           >  1:2401052/1‑71 (MQ=255)
      gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgc                                                           >  1:4098853/1‑71 (MQ=255)
      gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcg                                                            >  1:2299453/1‑70 (MQ=255)
             cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGAc                                                                   <  1:3576545/56‑1 (MQ=255)
             cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                    <  1:1039876/71‑1 (MQ=255)
                tttAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgcgc                                                         >  1:1055040/1‑63 (MQ=255)
                  tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGa                                                                    >  1:2789213/1‑50 (MQ=255)
                  tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                    >  1:656168/1‑66 (MQ=255)
                     aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                            <  1:3453558/71‑1 (MQ=255)
                     aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                            <  1:3104870/71‑1 (MQ=255)
                     aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                            <  1:735565/71‑1 (MQ=255)
                          ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACg                                                                  <  1:1639053/44‑1 (MQ=255)
                          ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACg                                                                  <  1:860012/44‑1 (MQ=255)
                          ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                                    >  1:2727446/1‑58 (MQ=255)
                            gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGAt                                                                <  1:1261431/44‑1 (MQ=255)
                            gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcg                                                            >  1:2007165/1‑48 (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:2011318/69‑1 (MQ=255)
                                    taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                              <  1:1860825/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: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:227987/42‑1 (MQ=255)
                                                             aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGaa                                 <  1:1133831/42‑1 (MQ=255)
                                                             aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                  >  1:3483104/1‑41 (MQ=255)
                                                             aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAt                                >  1:3235777/1‑43 (MQ=255)
                                                             aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAt                                >  1:1691153/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)
                                                               gCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTc  <  1:894359/71‑1 (MQ=255)
                                                                cTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTc  <  1:169050/70‑1 (MQ=255)
                                                                |                                                                     
GGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTC  >  minE/1933872‑1934005

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: