Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,118,451 G→A 6.0% A162V (GCG→GTG)  ydiU ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,118,4510GA 6.0% 203.1 / 2.5 66A162V (GCG→GTG) ydiUconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/2);  total (34/32)
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.22e-01

AAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGAT  >  minE/1118383‑1118520
                                                                    |                                                                     
aaaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgcc                                                                     >  1:2341774/1‑71 (MQ=255)
 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                                    >  1:934468/1‑71 (MQ=255)
 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                                    >  1:2702688/1‑71 (MQ=255)
 aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                                    >  1:849264/1‑71 (MQ=255)
  aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGt                                                                   <  1:4438858/71‑1 (MQ=255)
   aaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGtt                                                                  >  1:4217380/1‑71 (MQ=255)
     tGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTc                                                                >  1:2794939/1‑71 (MQ=255)
       ttCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                              <  1:2121411/71‑1 (MQ=255)
         cGAAATTACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                             <  1:2178618/70‑1 (MQ=255)
         cGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                             <  1:1764701/70‑1 (MQ=255)
         cGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                            >  1:1689918/1‑71 (MQ=255)
         cGAAATGACCAAAGCGCAGATGACTTGGTGCAACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                            >  1:3379215/1‑71 (MQ=255)
           aaaTGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATa                                                          >  1:2361756/1‑71 (MQ=255)
                aCCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTg                                                     <  1:3155759/71‑1 (MQ=255)
                   aaaGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                              <  1:4083480/59‑1 (MQ=255)
                       cgcAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCTTTTCCCGATACACTGGCGAATc                                              <  1:2357817/71‑1 (MQ=255)
                         cAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCg                                             <  1:2794847/70‑1 (MQ=255)
                          aGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTcc                                                               >  1:1650615/1‑51 (MQ=255)
                             tGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATaca                                                        >  1:2616093/1‑55 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTcc                                                               >  1:3202578/1‑47 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                            >  1:4550111/1‑50 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                        >  1:3649497/1‑70 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                        >  1:562364/1‑70 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                        >  1:2817458/1‑70 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                        >  1:1691475/1‑70 (MQ=255)
                              gACTTGGTGCCACACGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGCTGGt                                        <  1:714957/70‑1 (MQ=255)
                               aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                       <  1:4536341/70‑1 (MQ=255)
                               aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                       <  1:1694466/70‑1 (MQ=255)
                               aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                      <  1:280709/71‑1 (MQ=255)
                               aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                      <  1:4274794/71‑1 (MQ=255)
                                cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                       <  1:4316446/69‑1 (MQ=255)
                                cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGAc                                     <  1:139018/71‑1 (MQ=255)
                                  tGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                                   >  1:3008774/1‑71 (MQ=255)
                                   ggTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                                   >  1:348860/1‑70 (MQ=255)
                                    gTGCCACACGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGCTGGTGACGAt                                  <  1:3725507/70‑1 (MQ=255)
                                       ccACACGCATCAGCATCTCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACtt                              <  1:3090433/71‑1 (MQ=255)
                                         acacGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGc                                            >  1:2334425/1‑55 (MQ=255)
                                         acacGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTa                             >  1:753958/1‑70 (MQ=255)
                                           acGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcg                          >  1:3109144/1‑71 (MQ=255)
                                            cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGc                                                   >  1:2168606/1‑45 (MQ=255)
                                            cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATAc                                >  1:1435511/1‑64 (MQ=255)
                                            cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                         <  1:2171750/71‑1 (MQ=255)
                                            cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                         <  1:1846979/71‑1 (MQ=255)
                                            cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                         >  1:4545402/1‑71 (MQ=255)
                                            cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                         <  1:1617045/71‑1 (MQ=255)
                                             gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                        >  1:709789/1‑71 (MQ=255)
                                             gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                        <  1:4379989/71‑1 (MQ=255)
                                             gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                        >  1:1716141/1‑71 (MQ=255)
                                              cATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcgc                       <  1:4565444/71‑1 (MQ=255)
                                              cATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                        <  1:3078624/70‑1 (MQ=255)
                                                 cAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGc                                            <  1:3693050/47‑1 (MQ=255)
                                                 cAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAc                           >  1:3690947/1‑64 (MQ=255)
                                                  aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGTTGGTGACGATACTTAACGCGCggg                    >  1:4446537/1‑70 (MQ=255)
                                                  aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                            <  1:306003/62‑1 (MQ=255)
                                                  aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgt                   >  1:2724991/1‑71 (MQ=255)
                                                      tCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgt                >  1:4296942/1‑70 (MQ=255)
                                                        gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTTACGATACTTAACGCGCGGGTCGTCgg             <  1:1426267/71‑1 (MQ=255)
                                                        gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg             <  1:933094/71‑1 (MQ=255)
                                                        gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg             <  1:3197870/71‑1 (MQ=255)
                                                        gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTCGTGACGATACTTAACGCGCGGGTCGTCgg             <  1:4412169/71‑1 (MQ=255)
                                                         cgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa            <  1:4305943/71‑1 (MQ=255)
                                                         cgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa            <  1:1705868/71‑1 (MQ=255)
                                                         cgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa            <  1:2947960/71‑1 (MQ=255)
                                                              ggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                                   >  1:3225689/1‑43 (MQ=255)
                                                              ggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc       >  1:3325683/1‑71 (MQ=255)
                                                               gCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAAt          >  1:153602/1‑67 (MQ=255)
                                                                cTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc       >  1:452392/1‑69 (MQ=255)
                                                                cTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc       >  1:961029/1‑69 (MQ=255)
                                                                cTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc       >  1:2859559/1‑69 (MQ=255)
                                                                 tCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc       >  1:4384031/1‑68 (MQ=255)
                                                                  ccgccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAg    <  1:3281071/70‑1 (MQ=255)
                                                                   cgccgTTTCCCGATACACTGGCGAATCGCTGGTGAc                                     <  1:4251742/36‑1 (MQ=255)
                                                                   cgccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg             <  1:65595/60‑1 (MQ=255)
                                                                   cgccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGAt  <  1:575348/71‑1 (MQ=255)
                                                                    |                                                                     
AAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGAT  >  minE/1118383‑1118520

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: