Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,759 G→A 6.6% N228N (AAC→AAT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,7590GA 6.6% 247.3 / 5.4 75N228N (AAC→AATlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base A (3/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

AGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTT  >  minE/1933689‑1933829
                                                                      |                                                                      
aGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGcgcg                                                                        >  1:4067940/1‑71 (MQ=255)
 gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCTCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                       <  1:47521/71‑1 (MQ=255)
 gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                       <  1:3471807/71‑1 (MQ=255)
 gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                       >  1:2937929/1‑71 (MQ=255)
  ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGtt                                                                      <  1:2867274/71‑1 (MQ=255)
    accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                    <  1:940860/71‑1 (MQ=255)
    accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                    <  1:3777217/71‑1 (MQ=255)
     ccaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtg                                                                   <  1:1778782/71‑1 (MQ=255)
       aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                                                                 >  1:2263193/1‑71 (MQ=255)
       aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                                                                 <  1:4006010/71‑1 (MQ=255)
        accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg                                                                >  1:860261/1‑71 (MQ=255)
        accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                                                                 >  1:4068966/1‑70 (MQ=255)
         caaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtgatgg                                                               >  1:2820590/3‑71 (MQ=255)
         caaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtgatgg                                                               >  1:3355044/3‑71 (MQ=255)
          caccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtggg                                                              >  1:1790564/1‑71 (MQ=255)
            ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg                                                                <  1:3342585/67‑1 (MQ=255)
            ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggtg                                                            <  1:3637099/71‑1 (MQ=255)
            ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggt                                                             >  1:2704845/1‑70 (MQ=255)
              aGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggtg                                                            >  1:1325526/1‑69 (MQ=255)
               gACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                         <  1:146971/71‑1 (MQ=255)
                aCGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGtaataa                                                        <  1:708954/71‑1 (MQ=255)
                aCGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaa                                                        <  1:1254616/71‑1 (MQ=255)
                  gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                                                                 >  1:3728706/1‑60 (MQ=255)
                    ttGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg                                                                >  1:2975726/1‑59 (MQ=255)
                      gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTa                                                                     >  1:2845528/1‑52 (MQ=255)
                      gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtg                                                                   >  1:4072739/1‑54 (MQ=255)
                            cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                       >  1:2624653/1‑44 (MQ=255)
                            cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCa                                             <  1:1590013/70‑1 (MQ=255)
                            cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                            >  1:2046856/1‑71 (MQ=255)
                             aGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                                                       <  1:3268151/59‑1 (MQ=255)
                             aGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGc                                           <  1:2405991/71‑1 (MQ=255)
                              gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                                                       <  1:3847449/58‑1 (MQ=255)
                               ttCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg                                         >  1:662680/1‑71 (MQ=255)
                                tCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtg                                                                   >  1:4006402/1‑44 (MQ=16)
                                 ccGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                    >  1:799523/1‑42 (MQ=255)
                                 ccGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                    >  1:1645162/1‑42 (MQ=255)
                                 ccGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAAc                                                      >  1:3298229/1‑56 (MQ=255)
                                 ccGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCa                                       >  1:237568/1‑71 (MQ=255)
                                  cGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCAc                                      <  1:870545/71‑1 (MQ=255)
                                  cGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCAc                                      <  1:3474181/71‑1 (MQ=255)
                                  cGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAt                                                                    <  1:3633974/41‑1 (MQ=16)
                                   ggCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg                                         <  1:3369360/67‑1 (MQ=255)
                                   ggCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCAcc                                     <  1:3105635/71‑1 (MQ=255)
                                    gcgcGATACGCAGGTACATGTCGAGATCCAGCGCGTTa                                                                     <  1:2051441/38‑1 (MQ=255)
                                     cgcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGc                                   <  1:2667051/71‑1 (MQ=255)
                                      gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                                  <  1:3123512/71‑1 (MQ=255)
                                        gATACTCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACg                                                     <  1:131490/50‑1 (MQ=255)
                                        gATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACg                                                     <  1:3825922/50‑1 (MQ=255)
                                          tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGa                                                   >  1:2506259/1‑50 (MQ=255)
                                          tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGc                                   <  1:2999373/66‑1 (MQ=255)
                                          tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAg                                <  1:838136/69‑1 (MQ=255)
                                           aCGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                            >  1:1155076/1‑56 (MQ=255)
                                           aCGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCa                                 >  1:142193/1‑67 (MQ=255)
                                            cGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                            >  1:3541252/1‑55 (MQ=255)
                                            cGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                             <  1:1136795/70‑1 (MQ=255)
                                             gCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATc                           <  1:884411/71‑1 (MQ=255)
                                             gCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATc                           <  1:1747792/71‑1 (MQ=255)
                                               aGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAc                         <  1:2002269/71‑1 (MQ=255)
                                                 gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                         <  1:3805697/37‑1 (MQ=255)
                                                 gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAggg                              <  1:1973442/64‑1 (MQ=255)
                                                  tACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGc                                   <  1:93284/58‑1 (MQ=255)
                                                   aCATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                                                       <  1:1291034/37‑1 (MQ=255)
                                                   aCATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGa                                                   >  1:1211037/1‑41 (MQ=255)
                                                   aCATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCa                                             >  1:182780/1‑47 (MQ=255)
                                                   aCATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCa                                 >  1:1299086/1‑59 (MQ=255)
                                                     aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACg                                                     >  1:1346350/1‑37 (MQ=255)
                                                       gTCGAGATACAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATc                           >  1:2050277/1‑61 (MQ=255)
                                                           agaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCgg             <  1:3994758/71‑1 (MQ=255)
                                                            gaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                            >  1:353782/1‑39 (MQ=255)
                                                            gaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGAt                            >  1:1468868/1‑55 (MQ=255)
                                                            gaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAcc                        >  1:1717186/1‑59 (MQ=255)
                                                            gaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGc            >  1:467328/1‑71 (MQ=255)
                                                             aTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCg           >  1:2762374/1‑71 (MQ=255)
                                                              tCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGt          <  1:938572/71‑1 (MQ=255)
                                                               ccAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGc                                           >  1:937731/1‑37 (MQ=255)
                                                               ccAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGtt         >  1:1948045/1‑71 (MQ=255)
                                                                   cgcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCa                                 <  1:858353/43‑1 (MQ=255)
                                                                    gcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCaa     >  1:3396438/1‑70 (MQ=255)
                                                                    gcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAAc    >  1:1204740/1‑71 (MQ=255)
                                                                     cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAAc    >  1:2137886/1‑70 (MQ=255)
                                                                     cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAAc    >  1:250027/1‑70 (MQ=255)
                                                                      gTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACtt  <  1:594159/71‑1 (MQ=255)
                                                                      |                                                                      
AGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTT  >  minE/1933689‑1933829

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