Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,049 C→T 26.9% N12N (AAC→AAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0490CT26.9% 255.6 / 102.8 145N12N (AAC→AATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (53/53);  new base T (21/18);  total (74/71)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 146 total aligned reads displayed.

GGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAG  >  minE/217980‑218118
                                                                     |                                                                     
    aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                                   >  1:4056465/1‑70 (MQ=255)
     aTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAAc                                                                       <  1:3589637/65‑1 (MQ=255)
      tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                                   >  1:484529/1‑68 (MQ=255)
      tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                                   >  1:1863877/1‑68 (MQ=255)
      tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                                <  1:327962/71‑1 (MQ=255)
       aaGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                                >  1:3768469/1‑70 (MQ=255)
          cTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAat                                                             >  1:2275791/1‑70 (MQ=255)
             cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGa                                                               >  1:970008/1‑65 (MQ=255)
                  gCGTCTGGACGTGGTAATGAATGAACACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                    >  1:3999352/1‑71 (MQ=255)
                   cGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAgg                                                   <  1:3545393/71‑1 (MQ=255)
                    gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAgg                                                   >  1:1300484/1‑70 (MQ=255)
                     tCTGGACGTTGTAATGAATTAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                       <  1:224869/65‑1 (MQ=255)
                     tCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                       <  1:2602566/65‑1 (MQ=255)
                     tCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGa                                                 <  1:541101/71‑1 (MQ=255)
                       tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                                              >  1:1040555/1‑56 (MQ=255)
                       tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                                                >  1:3594168/1‑70 (MQ=255)
                       tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                               >  1:422698/1‑71 (MQ=255)
                         gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatt                                                            <  1:1584520/56‑1 (MQ=255)
                         gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                                              <  1:3889417/70‑1 (MQ=255)
                           cGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                                            <  1:3190175/70‑1 (MQ=255)
                            gTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                          <  1:2592846/71‑1 (MQ=255)
                            gTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                          >  1:3332144/1‑71 (MQ=255)
                              ggTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatt                                                            <  1:3949504/51‑1 (MQ=39)
                               gTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                                                     >  1:1325497/1‑57 (MQ=255)
                               gTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                                       >  1:1807392/1‑71 (MQ=255)
                               gTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                                       <  1:3399054/71‑1 (MQ=255)
                                  atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                    >  1:889177/1‑71 (MQ=255)
                                  atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                    >  1:3557409/1‑71 (MQ=255)
                                   tgaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                                              <  1:2460343/60‑1 (MQ=255)
                                    gaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                                       >  1:2789979/1‑66 (MQ=255)
                                    gaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa                                  >  1:1154347/1‑71 (MQ=255)
                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                       >  1:3305832/1‑49 (MQ=39)
                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                       >  1:844124/1‑49 (MQ=39)
                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                                       >  1:2725620/1‑65 (MQ=255)
                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa                                  >  1:2819341/1‑70 (MQ=255)
                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                 >  1:247012/1‑71 (MQ=255)
                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                 >  1:1118129/1‑71 (MQ=255)
                                     aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGa                                                 >  1:1713938/1‑55 (MQ=255)
                                     aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                 >  1:3609852/1‑71 (MQ=255)
                                      atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGa                                                               <  1:1689835/40‑1 (MQ=38)
                                      atgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                 >  1:2960934/1‑70 (MQ=255)
                                      atgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAg                                                    <  1:442812/51‑1 (MQ=255)
                                       tgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                 >  1:884783/1‑69 (MQ=255)
                                       tgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                                <  1:3743650/70‑1 (MQ=255)
                                        gaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAggg                                                  <  1:1506492/51‑1 (MQ=255)
                                         aaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                              >  1:2244568/1‑70 (MQ=255)
                                         aaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                              <  1:1125490/70‑1 (MQ=255)
                                          aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                                            >  1:1926202/1‑55 (MQ=255)
                                          aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTAAGGGATACGGCACATCGCTAtt                                <  1:971237/67‑1 (MQ=255)
                                          aGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                                          <  1:2834492/57‑1 (MQ=255)
                                            acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAattatt                                                         <  1:3100439/40‑1 (MQ=38)
                                             cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                       >  1:1403203/1‑41 (MQ=38)
                                             cgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATAttt                         <  1:1071762/71‑1 (MQ=255)
                                              gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                                                 >  1:910397/1‑46 (MQ=38)
                                              gacTACAAAATAAGTAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                                        >  1:1207762/1‑55 (MQ=255)
                                              gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                        >  1:4099478/1‑71 (MQ=255)
                                               acTACAAAATAAGAAGAGTAAACGCAGCAGAATTATTTTCAGGGAt                                                <  1:267106/46‑1 (MQ=255)
                                               acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                                                <  1:29086/46‑1 (MQ=255)
                                               acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                                                <  1:2962321/46‑1 (MQ=255)
                                               acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                                          <  1:891890/52‑1 (MQ=255)
                                                cTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGa                      >  1:3510357/1‑71 (MQ=255)
                                                 tgcAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                          <  1:3960874/48‑1 (MQ=39)
                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                       >  1:579730/1‑37 (MQ=255)
                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                       >  1:2306284/1‑37 (MQ=255)
                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAgg                                                   <  1:757129/41‑1 (MQ=255)
                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTa                               <  1:616644/61‑1 (MQ=255)
                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                     >  1:1151283/1‑71 (MQ=255)
                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                     >  1:1412998/1‑71 (MQ=255)
                                                 tACAAAATAAGAAGAGGAAACCCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                     >  1:3452223/1‑71 (MQ=255)
                                                   cAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                       >  1:3559184/1‑35 (MQ=255)
                                                   cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                    >  1:3335661/1‑70 (MQ=255)
                                                   cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                    >  1:1544234/1‑70 (MQ=255)
                                                   cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                    >  1:2184564/1‑70 (MQ=255)
                                                     aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                                                 <  1:3898359/39‑1 (MQ=255)
                                                     aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                 >  1:3252562/1‑71 (MQ=255)
                                                       aTAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTa                                  >  1:1226933/1‑52 (MQ=255)
                                                        tAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGTTATTaa                              <  1:3987488/55‑1 (MQ=255)
                                                        tAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag              <  1:843026/71‑1 (MQ=255)
                                                        tAAGAAGAGGAAACGCAGCAGAATAATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                 >  1:2261194/1‑68 (MQ=255)
                                                         aagaagAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga             <  1:2361289/71‑1 (MQ=255)
                                                          agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGcaca                                        >  1:2637555/1‑43 (MQ=255)
                                                          agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata             <  1:3397028/70‑3 (MQ=255)
                                                          agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgat              >  1:2761210/1‑68 (MQ=255)
                                                          agaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga             <  1:3523578/70‑1 (MQ=255)
                                                            aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg          <  1:1545557/71‑1 (MQ=255)
                                                            aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg          <  1:716141/71‑1 (MQ=255)
                                                             agagGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAgg          >  1:1454612/1‑70 (MQ=255)
                                                             agagGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg          <  1:4002204/70‑1 (MQ=255)
                                                              gagGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                                       >  1:2418975/1‑40 (MQ=255)
                                                              gagGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                 <  1:158657/46‑1 (MQ=255)
                                                               agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTCGACGAATGATAAGGTAt       >  1:530055/1‑71 (MQ=255)
                                                               agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGa                   <  1:2462396/59‑1 (MQ=255)
                                                               agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt       <  1:1766001/71‑1 (MQ=255)
                                                               agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt       <  1:2711932/71‑1 (MQ=255)
                                                                ggAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTAt       <  1:1519628/70‑1 (MQ=255)
                                                                 gAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc     <  1:1069743/71‑1 (MQ=255)
                                                                 gAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc     <  1:1858237/71‑1 (MQ=255)
                                                                  aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc     <  1:2087874/70‑1 (MQ=255)
                                                                  aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc     >  1:2414334/1‑70 (MQ=255)
                                                                   aaCGCAGAAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg          <  1:335349/64‑1 (MQ=255)
                                                                     |                                                                     
GGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAG  >  minE/217980‑218118

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: