Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,052 A→C 8.8% A13A (GCA→GCC yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0520AC 8.8% 337.4 / 30.0 137A13A (GCA→GCCyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (61/64);  new base C (8/4);  total (69/68)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 149 total aligned reads displayed.

GAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAG  >  minE/217983‑218122
                                                                     |                                                                      
 aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                                       >  1:4056465/1‑70 (MQ=255)
   tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                                       >  1:1863877/1‑68 (MQ=255)
   tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                                       >  1:484529/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)
       cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                                    <  1:247099/67‑1 (MQ=255)
          cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGa                                                                   >  1:970008/1‑65 (MQ=255)
             tggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAat                                                                 <  1:2760757/63‑1 (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)
                   cTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                           <  1:2803704/64‑1 (MQ=255)
                    tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                                                  >  1:1040555/1‑56 (MQ=255)
                    tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                                   >  1:1841463/1‑71 (MQ=255)
                    tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                                   >  1:422698/1‑71 (MQ=255)
                      gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattatt                                                             <  1:3994492/59‑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:3399054/71‑1 (MQ=255)
                            gTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                                           >  1:1807392/1‑71 (MQ=255)
                             tAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                           <  1:535521/54‑1 (MQ=255)
                             tAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                           <  1:1343104/54‑1 (MQ=255)
                               atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                        >  1:889177/1‑71 (MQ=255)
                               atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                        >  1:3557409/1‑71 (MQ=255)
                               atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                        >  1:1471576/1‑71 (MQ=255)
                               atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                        >  1:3641691/1‑71 (MQ=255)
                                tgaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                                                  <  1:2460343/60‑1 (MQ=255)
                                 gaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa                                      >  1:1154347/1‑71 (MQ=255)
                                  aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                           >  1:844124/1‑49 (MQ=39)
                                  aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa                                      >  1:2819341/1‑70 (MQ=255)
                                  aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                     >  1:1118129/1‑71 (MQ=255)
                                  aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGa                                                     >  1:1713938/1‑55 (MQ=255)
                                  aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                     >  1:3609852/1‑71 (MQ=255)
                                  aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                     >  1:2690485/1‑71 (MQ=255)
                                   atgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                     >  1:2960934/1‑70 (MQ=255)
                                    tgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                     >  1:884783/1‑69 (MQ=255)
                                      aaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                                  >  1:2244568/1‑70 (MQ=255)
                                      aaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                                            >  1:3590872/1‑60 (MQ=255)
                                      aaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                                  <  1:1125490/70‑1 (MQ=255)
                                       aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                                                >  1:1926202/1‑55 (MQ=255)
                                       aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAat                                                                 <  1:2990962/38‑1 (MQ=37)
                                       aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTAAGGGATACGGCACATCGCTAtt                                    <  1:971237/67‑1 (MQ=255)
                                        gacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa                                      <  1:469087/64‑1 (MQ=255)
                                        gacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                               <  1:3859721/71‑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:2617868/71‑1 (MQ=255)
                                            acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                                                    <  1:29086/46‑1 (MQ=255)
                                            acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                                              <  1:891890/52‑1 (MQ=255)
                                             cTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                                            >  1:3251094/1‑53 (MQ=255)
                                             cTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                                            >  1:160648/1‑53 (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)
                                              tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                                                    <  1:1523831/44‑1 (MQ=255)
                                              tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTa                                   <  1:616644/61‑1 (MQ=255)
                                              tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                         >  1:1151283/1‑71 (MQ=255)
                                                cAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                                           >  1:3559184/1‑35 (MQ=255)
                                                cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                        <  1:4025700/70‑1 (MQ=255)
                                                cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                        >  1:3335661/1‑70 (MQ=255)
                                                cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACgg                                                >  1:289175/1‑46 (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)
                                                     tAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                              <  1:1715851/43‑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:716141/71‑1 (MQ=255)
                                                         aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg              <  1:1545557/71‑1 (MQ=255)
                                                          agagGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAgg              >  1:1454612/1‑70 (MQ=255)
                                                          agagGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg              <  1:1554768/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:2711932/71‑1 (MQ=255)
                                                             ggAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTAt           <  1:1519628/70‑1 (MQ=255)
                                                             ggAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGat                      <  1:3280607/59‑2 (MQ=255)
                                                              gAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc         <  1:1069743/71‑1 (MQ=255)
                                                               aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc         <  1:2087874/70‑1 (MQ=255)
                                                               aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc         >  1:2414334/1‑70 (MQ=255)
                                                                aaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg              >  1:4120883/1‑64 (MQ=255)
                                                                aaCGCAGAAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg              <  1:335349/64‑1 (MQ=255)
                                                                     agcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                     >  1:3331619/1‑52 (MQ=255)
                                                                     agcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAg  <  1:3959747/71‑1 (MQ=255)
                                                                     |                                                                      
GAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAG  >  minE/217983‑218122

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: