Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,133 A→C 10.0% L85L (CTT→CTG yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,1330AC10.0% 128.8 / 12.6 60L85L (CTT→CTGyajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/36);  new base C (2/4);  total (20/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

GCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAAC  >  minE/253066‑253196
                                                                   |                                                               
gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGtaaag                                                               <  1:719261/70‑1 (MQ=255)
gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                                               >  1:4614777/1‑70 (MQ=255)
gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                                               >  1:2865832/1‑70 (MQ=255)
gCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag                                                               >  1:3901398/1‑70 (MQ=255)
 cTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaaga                                                              <  1:2626579/70‑1 (MQ=255)
 cTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaaga                                                              <  1:3089030/70‑1 (MQ=255)
 cTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaaga                                                              <  1:4108210/70‑1 (MQ=255)
    gcgGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag                                                           <  1:866931/70‑1 (MQ=255)
     cgGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTg                                                         <  1:4687668/71‑1 (MQ=255)
           tAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaa                                                             <  1:2747707/61‑1 (MQ=255)
             aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                                                 >  1:2702795/1‑71 (MQ=255)
             aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                                                 >  1:212789/1‑71 (MQ=255)
             aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca                                                 >  1:430510/1‑71 (MQ=255)
             aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcac                                                  <  1:315668/70‑1 (MQ=255)
             aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcac                                                  <  1:1484874/70‑1 (MQ=255)
             aCAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcac                                                  <  1:1033690/70‑1 (MQ=255)
              cAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACAt                                                <  1:99663/71‑1 (MQ=255)
               aGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATc                                                      >  1:1689288/1‑64 (MQ=255)
                gACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca                                                   >  1:4552861/1‑66 (MQ=255)
                gACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGa                                              >  1:4604117/1‑71 (MQ=255)
                  cAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACa                                            >  1:4022855/1‑71 (MQ=255)
                    aaTTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaa                                                            <  1:1869637/53‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag                                                           <  1:3250545/46‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:4225871/69‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:2980141/69‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:1151928/69‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:1820171/69‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:1903508/69‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:1391312/69‑1 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  >  1:1838916/1‑71 (MQ=255)
                            tAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:3739492/71‑1 (MQ=255)
                              gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  >  1:25562/1‑69 (MQ=255)
                              gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATa                                <  1:3161041/71‑1 (MQ=255)
                              gCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATa                                <  1:1471233/71‑1 (MQ=255)
                                ggAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAgg                              >  1:3452352/1‑71 (MQ=255)
                                 gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg                                     <  1:840047/63‑1 (MQ=255)
                                 gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg                                     <  1:3590256/63‑1 (MQ=255)
                                 gAAAATCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCg                                     <  1:53627/63‑1 (MQ=255)
                                  aaaaTCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGgcg                            >  1:3221462/1‑71 (MQ=255)
                                  aaaaTCGGGCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGgcg                            >  1:1193590/1‑71 (MQ=255)
                                  aaaGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:2964105/63‑1 (MQ=255)
                                   aaGTCCGTCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGgcgc                           <  1:452907/71‑1 (MQ=255)
                                   aaGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGgcgc                           <  1:3661874/71‑1 (MQ=255)
                                     gTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:1387258/62‑1 (MQ=255)
                                      tCCGGCGACAGCTGCGTCCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAAc                        <  1:4210695/71‑1 (MQ=255)
                                      tCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:2923251/61‑1 (MQ=255)
                                       ccGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACAt                                                >  1:336821/1‑46 (MQ=255)
                                        gggCGACAGCTGCGACCAGTGGCGGGACAGAAAGTGATCCCACATCACATCCAGCGTAATAGGCGCGAcac                      <  1:538874/70‑3 (MQ=255)
                                        cGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATcc                                                     >  1:1792563/1‑40 (MQ=255)
                                         ggCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca                                                   <  1:2553607/41‑1 (MQ=255)
                                               aGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCaa                         <  1:4159342/61‑1 (MQ=255)
                                                gCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:4645886/49‑1 (MQ=255)
                                                 cTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGttt             <  1:4619550/71‑1 (MQ=255)
                                                 cTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGttt             <  1:1920091/71‑1 (MQ=255)
                                                   gCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTc            >  1:1516869/1‑70 (MQ=255)
                                                     gACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCAc          >  1:753007/1‑70 (MQ=255)
                                                     gACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACt         <  1:3199801/71‑1 (MQ=255)
                                                      aCCAGTGGCGGGAAAGAAAGTGATCCCACATGACATc                                          <  1:802504/37‑1 (MQ=255)
                                                        cAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt                                    <  1:3180278/41‑1 (MQ=255)
                                                         aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCTTAATAGGCGCAACGCGGCGCGTTTCACTACGa     <  1:2282900/71‑1 (MQ=255)
                                                         aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  >  1:818307/1‑42 (MQ=255)
                                                         aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCaa                         <  1:1233666/51‑1 (MQ=255)
                                                           tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:422384/40‑1 (MQ=255)
                                                           tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:4292858/40‑1 (MQ=255)
                                                           tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACt         <  1:3164160/65‑1 (MQ=255)
                                                             gCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  >  1:1575056/1‑38 (MQ=255)
                                                             gCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAAc  >  1:655327/1‑70 (MQ=255)
                                                               gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  >  1:4333926/1‑36 (MQ=255)
                                                                ggAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:606373/35‑1 (MQ=255)
                                                                ggAAAGAAAGTGATCCCACATGACATCCAGCGTaa                                  <  1:1415144/35‑1 (MQ=255)
                                                                ggAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAAc  <  1:452597/67‑1 (MQ=255)
                                                                  aaagaaagTGATCCCACATGACATCCAGCGTAATAGGCGCa                          >  1:2959334/1‑41 (MQ=255)
                                                                   aagaaagTGATCCCACATGACATCCAGCGTAATAGGCGCaa                         >  1:1396796/1‑41 (MQ=255)
                                                                   aagaaagTGATCCCACATGACATCCAGCGTAATAGGCGCaa                         >  1:1868362/1‑41 (MQ=255)
                                                                   |                                                               
GCTCGCGGGCATAACAGACAAATTCCTGTAGCGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAGGCGCAACGCGGCGCGTTTCACTACGAAAC  >  minE/253066‑253196

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: