Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,052 A→C 14.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,0520AC14.8% 202.3 / 37.5 101A13A (GCA→GCCyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/40);  new base C (9/6);  total (55/46)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 110 total aligned reads displayed.

TAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAG  >  minE/217986‑218122
                                                                  |                                                                      
tAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcag                                                                       >  1:2739297/1‑68 (MQ=255)
   gCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                                                 <  1:2540746/71‑1 (MQ=255)
    cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                                                 >  1:1681841/1‑70 (MQ=255)
     tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                                               >  1:3952469/1‑71 (MQ=255)
     tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatt                                                                >  1:1169153/1‑70 (MQ=255)
       cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattat                                                              <  1:616730/70‑1 (MQ=255)
         ctggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                                                           >  1:299572/3‑71 (MQ=255)
         caGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                           <  1:1367906/71‑1 (MQ=255)
          aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                                    <  1:3086400/61‑1 (MQ=255)
          aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                           >  1:2538903/1‑70 (MQ=255)
            gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                           >  1:309878/1‑68 (MQ=255)
            gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                           >  1:3619811/1‑68 (MQ=255)
              gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                                                  >  1:3437745/1‑59 (MQ=255)
                cTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                                                     >  1:334325/1‑70 (MQ=255)
                 tGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                                          >  1:2199733/1‑64 (MQ=255)
                 tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                                   >  1:3495409/1‑71 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcagc                                                                      >  1:3611270/1‑50 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACg                                                 >  1:535588/1‑71 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                                                 <  1:2514967/55‑1 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                                   <  1:1624961/69‑1 (MQ=255)
                    aCGTGTTAATTAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                                                <  1:750625/71‑1 (MQ=255)
                     cGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                                                         >  1:1290340/1‑61 (MQ=255)
                       tGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                                              >  1:3189581/1‑70 (MQ=255)
                        ggTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                                            >  1:907077/1‑71 (MQ=255)
                         gTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                                           <  1:41141/71‑1 (MQ=255)
                          tAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAtt                                          >  1:2290569/1‑71 (MQ=255)
                          aaatgaatgAAGATGCCAGCCGGATAAGAAGAGGAAACgca                                                                        <  1:1021342/40‑1 (MQ=255)
                           aatgaatgAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                         >  1:3231337/1‑71 (MQ=255)
                            atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCa                                                         <  1:757557/54‑1 (MQ=255)
                            atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                        >  1:3676027/1‑55 (MQ=255)
                            atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                        <  1:1563068/71‑1 (MQ=255)
                             tgaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                                       <  1:1235152/71‑1 (MQ=255)
                               aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                        >  1:1358576/1‑52 (MQ=255)
                                atgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatta                                                               <  1:2802961/44‑1 (MQ=38)
                                 tgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                        <  1:24334/50‑1 (MQ=39)
                                  gaaGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                                    <  1:470995/37‑1 (MQ=255)
                                   aaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                                  >  1:133170/1‑70 (MQ=255)
                                    cgacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGAt                                                    <  1:567610/50‑1 (MQ=38)
                                    cgacgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTc                                                          <  1:2962830/44‑1 (MQ=38)
                                    cgacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                                                  <  1:1144971/52‑1 (MQ=255)
                                    aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAgg                                                       <  1:814196/48‑1 (MQ=39)
                                    aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACg                                                 <  1:3830853/54‑1 (MQ=255)
                                    aGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattatt                                                             <  1:1573422/42‑1 (MQ=255)
                                     gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                               >  1:121127/1‑71 (MQ=255)
                                     gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                               >  1:3281020/1‑71 (MQ=255)
                                     gacgacTACAAAATAAGAAGAGGGAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                               >  1:492523/1‑71 (MQ=255)
                                      acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                                          <  1:2079541/43‑1 (MQ=38)
                                      acgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                               >  1:1630863/1‑70 (MQ=255)
                                       cgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                                                         <  1:1548989/43‑1 (MQ=38)
                                        gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                            >  1:1274865/1‑71 (MQ=255)
                                        gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                            <  1:3392271/71‑1 (MQ=255)
                                         acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                            <  1:2097155/70‑1 (MQ=255)
                                           tgcAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                              <  1:1160048/48‑1 (MQ=39)
                                           tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCa                                                         >  1:863718/1‑39 (MQ=255)
                                             cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGGACATTGCTATTAATATTTTGACg                        <  1:1249133/70‑1 (MQ=255)
                                             cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                        >  1:2679737/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                        >  1:3446679/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGAt                                                    <  1:1649140/42‑1 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                        >  1:1025623/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                        <  1:26962/70‑1 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                        >  1:1587024/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                        >  1:1962411/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGa                       <  1:3178782/71‑1 (MQ=255)
                                              aaaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGaa                      >  1:2943575/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCa                                                         <  1:2566579/35‑1 (MQ=255)
                                               aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                     >  1:1162456/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                     >  1:1282403/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt                     >  1:2635144/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt                     >  1:2080023/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                                       >  1:3244053/1‑53 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                                       >  1:3248726/1‑53 (MQ=255)
                                                aaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                     >  1:714974/1‑70 (MQ=255)
                                                 aTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATga                   >  1:746633/1‑71 (MQ=255)
                                                   aagaagAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag                  >  1:1398727/1‑70 (MQ=255)
                                                    agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata                 <  1:3797051/70‑3 (MQ=255)
                                                    agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata                 <  1:629809/70‑3 (MQ=255)
                                                     gaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAat                                 <  1:3589853/53‑1 (MQ=255)
                                                     gaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAg               <  1:3688200/71‑1 (MQ=255)
                                                      aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAtt                                          <  1:2948928/43‑1 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                         <  1:3535227/43‑1 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                         <  1:83172/43‑1 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg              >  1:1311194/1‑70 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg              >  1:3544892/1‑70 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg              >  1:3544019/1‑70 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGt             <  1:2237575/71‑1 (MQ=255)
                                                        gagGAAATGCAGAAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTa            >  1:468475/1‑71 (MQ=255)
                                                         agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                         <  1:2126842/41‑1 (MQ=255)
                                                         agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAt           <  1:148946/71‑1 (MQ=255)
                                                         agGAAATGCAGCAGAATTATTATCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAt           >  1:2042261/1‑71 (MQ=255)
                                                         agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt           >  1:2218775/1‑71 (MQ=255)
                                                          ggAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTAtt          >  1:590098/1‑71 (MQ=255)
                                                           gAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTa            >  1:2488924/1‑68 (MQ=255)
                                                            aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                                    <  1:1810703/43‑1 (MQ=255)
                                                            aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                                    <  1:1234284/43‑1 (MQ=255)
                                                            aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATaa                >  1:3004111/1‑63 (MQ=255)
                                                            aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                     <  1:3268993/58‑1 (MQ=255)
                                                            aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGt             <  1:2092416/66‑1 (MQ=255)
                                                             aatgcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc         >  1:1713405/4‑69 (MQ=255)
                                                               tgcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAgg     <  1:397642/70‑1 (MQ=255)
                                                                gcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGc    <  1:357132/71‑1 (MQ=255)
                                                                  |                                                                      
TAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAG  >  minE/217986‑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: