Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,030 A→G 30.4% Y6C (TAC→TGC)  yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0300AG30.4% 157.2 / 98.3 115Y6C (TAC→TGC) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base G (21/14);  total (63/52)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 115 total aligned reads displayed.

CGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAC  >  minE/217963‑218099
                                                                   |                                                                     
cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa                                                                     >  1:2686778/1‑70 (MQ=255)
           gCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaag                                                           <  1:427443/69‑1 (MQ=255)
              cGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                                                       >  1:868611/1‑70 (MQ=255)
                tGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAATATaaga                                                             <  1:77/62‑1 (MQ=255)
                 ggAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACg                                                   >  1:3494015/1‑71 (MQ=255)
                  gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaatg                                                   >  1:1985133/1‑68 (MQ=255)
                    gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                <  1:1101612/71‑1 (MQ=255)
                     aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt                                                                 <  1:2377378/53‑1 (MQ=255)
                     aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                >  1:4056465/1‑70 (MQ=255)
                       tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaa                                                               >  1:233314/1‑53 (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)
                              cactGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaa                                                     >  1:3968423/1‑56 (MQ=255)
                              cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGa                                            >  1:970008/1‑65 (MQ=255)
                                 tggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAat                                          <  1:2760757/63‑1 (MQ=255)
                                 acgCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    <  1:1275039/68‑1 (MQ=255)
                                   gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaa                                                      <  1:2391668/50‑1 (MQ=255)
                                   gCGTCTGGACGTGGTAATGAATGAACACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                 >  1:3999352/1‑71 (MQ=255)
                                    cGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAgg                                <  1:3545393/71‑1 (MQ=255)
                                     gTCTGGACGTGGTAATTAATGAAGACGACTACAAAATAAGAAGAGGa                                                       <  1:940012/47‑1 (MQ=255)
                                     gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAgg                                >  1:1300484/1‑70 (MQ=255)
                                      tCTGGACGTTGTAATGAATTAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    <  1:224869/65‑1 (MQ=255)
                                      tCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAgg                                                        <  1:2561576/45‑1 (MQ=255)
                                      tCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    <  1:2602566/65‑1 (MQ=255)
                                      tCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGa                              <  1:541101/71‑1 (MQ=255)
                                       cTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    <  1:2803704/64‑1 (MQ=255)
                                        tGGACGTGGTAATGAATGAAGACGACTACAAAATaagaa                                                            >  1:1102258/1‑39 (MQ=255)
                                        tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                           >  1:1040555/1‑56 (MQ=255)
                                        tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                             >  1:3594168/1‑70 (MQ=255)
                                        tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                            >  1:422698/1‑71 (MQ=255)
                                        tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                            >  1:1841463/1‑71 (MQ=255)
                                          gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattatt                                      <  1:3994492/59‑1 (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)
                                                 tAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                    <  1:535521/54‑1 (MQ=255)
                                                 tAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    <  1:1343104/54‑1 (MQ=255)
                                                  aatgaatgAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAattatt                                      >  1:1194112/1‑51 (MQ=255)
                                                   ttgaattaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAgg                                <  1:3333265/55‑1 (MQ=255)
                                                   atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                 >  1:3557409/1‑71 (MQ=255)
                                                   atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                 >  1:3641691/1‑71 (MQ=255)
                                                   atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                 >  1:1471576/1‑71 (MQ=255)
                                                    tgaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGa                                            <  1:3383201/43‑1 (MQ=255)
                                                     gaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAt                    >  1:2789979/1‑66 (MQ=255)
                                                      aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                    >  1:844124/1‑49 (MQ=39)
                                                      aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTAtttt                                    >  1:3305832/1‑49 (MQ=39)
                                                      aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                    >  1:2725620/1‑65 (MQ=255)
                                                      aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa               >  1:891636/1‑70 (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)
                                                      aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACgg                         >  1:225047/1‑60 (MQ=255)
                                                      aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt              >  1:3609852/1‑71 (MQ=255)
                                                      aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt              >  1:2690485/1‑71 (MQ=255)
                                                      aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt              >  1:2222208/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)
                                                       atgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                       >  1:629782/1‑61 (MQ=255)
                                                        tgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt              >  1:884783/1‑69 (MQ=255)
                                                        tgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt             <  1:3743650/70‑1 (MQ=255)
                                                         gaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAggg                               <  1:1506492/51‑1 (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)
                                                            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:2962321/46‑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)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    >  1:2306284/1‑37 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAgg                                <  1:757129/41‑1 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGAt                             <  1:1523831/44‑1 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                           <  1:91473/46‑1 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                >  1:3523939/1‑57 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTa            <  1:616644/61‑1 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt     <  1:3558693/68‑1 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc  >  1:1151283/1‑71 (MQ=255)
                                                                  tACAAAATAAGAAGAGGAAACCCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc  >  1:3452223/1‑71 (MQ=255)
                                                                   |                                                                     
CGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAC  >  minE/217963‑218099

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: