Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,166 C→T 18.7% Y51Y (TAC→TAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,1660CT18.7% 215.0 / 39.3 91Y51Y (TAC→TATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (5/12);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

AATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACA  >  minE/218101‑218234
                                                                 |                                                                    
aaTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGc                                                                 <  1:1917859/71‑1 (MQ=255)
 aTGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCg                                                                >  1:127160/1‑71 (MQ=255)
 aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCg                                                                <  1:1117020/71‑1 (MQ=255)
 aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCg                                                                >  1:3207457/1‑71 (MQ=255)
  tGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGt                                                               >  1:2223288/1‑71 (MQ=255)
       aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCt                                                                    >  1:465618/1‑61 (MQ=255)
         gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCt                                                        >  1:2293620/1‑71 (MQ=255)
         gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGAACAGAAACTACCTGGCGTCAGTCCt                                                        >  1:3230561/1‑71 (MQ=255)
          tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCtt                                                       >  1:2291435/1‑71 (MQ=255)
           aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTg                                                      >  1:3813829/1‑71 (MQ=255)
           aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                      >  1:807176/1‑71 (MQ=255)
           aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                      >  1:1203054/1‑71 (MQ=255)
           aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                      >  1:4119509/1‑71 (MQ=255)
              cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCt                                                                    <  1:3415655/54‑1 (MQ=255)
              cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                                    >  1:187658/1‑70 (MQ=255)
              cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCTTTGCgg                                                   >  1:675992/1‑71 (MQ=255)
                 ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                  >  1:1901107/1‑69 (MQ=255)
                   cAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGc                                              <  1:833421/71‑1 (MQ=255)
                    aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGt                                                               <  1:2840579/53‑1 (MQ=255)
                    aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCg                                             <  1:4057331/71‑1 (MQ=255)
                    aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCg                                             <  1:3230083/71‑1 (MQ=255)
                     gggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgggg                                                 <  1:4075973/66‑1 (MQ=255)
                      ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGc                                              <  1:3726520/68‑1 (MQ=255)
                       gTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGa                                                <  1:600666/65‑1 (MQ=255)
                       gTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGa                                                <  1:1968890/65‑1 (MQ=255)
                       gTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTAc                                                                      >  1:4132895/1‑43 (MQ=255)
                        ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                         <  1:196569/71‑1 (MQ=255)
                        ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                                   >  1:1851928/1‑61 (MQ=255)
                         tAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCtt                                        >  1:1841244/1‑71 (MQ=255)
                          aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                                       >  1:2480533/1‑71 (MQ=255)
                           aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                  <  1:1726861/59‑1 (MQ=255)
                           aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                   >  1:1368511/1‑58 (MQ=255)
                           aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                   >  1:377940/1‑58 (MQ=255)
                            gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     <  1:3654419/71‑1 (MQ=255)
                              cGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCt                                                                    >  1:2426795/1‑38 (MQ=38)
                              cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGTAGCGGGCTTTCgta                                   <  1:3869391/71‑1 (MQ=255)
                              cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    <  1:1015421/70‑1 (MQ=255)
                              cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    <  1:2551760/70‑1 (MQ=255)
                              cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     >  1:1410887/1‑69 (MQ=255)
                               gTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                                    <  1:3732799/53‑1 (MQ=255)
                               gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                            <  1:3817410/45‑1 (MQ=255)
                               gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                            <  1:1199974/45‑1 (MQ=255)
                               gTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAg                                                            <  1:1195059/45‑1 (MQ=255)
                                tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGa                                                >  1:753863/1‑56 (MQ=255)
                                tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:1559773/71‑1 (MQ=255)
                                tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:2904104/71‑1 (MQ=255)
                                tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:308399/71‑1 (MQ=255)
                                tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtag                                  >  1:3508145/1‑70 (MQ=255)
                                 aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagt                                 >  1:851487/1‑70 (MQ=255)
                                 aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     >  1:665057/1‑66 (MQ=255)
                                  aGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGc                                              <  1:849033/56‑1 (MQ=255)
                                  aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTc                                                              >  1:675423/1‑40 (MQ=255)
                                  aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtta                                                      <  1:3263540/48‑2 (MQ=255)
                                  aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtta                                                      <  1:3404136/48‑2 (MQ=255)
                                  aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGGGGGAGCGGGCTTTCGTAGTaa                               >  1:457266/1‑70 (MQ=255)
                                  aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                                   >  1:3920696/1‑51 (MQ=255)
                                   gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCggg                                                  >  1:434918/1‑51 (MQ=255)
                                   gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     <  1:3550700/64‑1 (MQ=255)
                                    aTGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTc                                                          <  1:3601167/42‑1 (MQ=255)
                                      gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                   >  1:3530399/1‑63 (MQ=255)
                                        gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCtta                                                      <  1:2863869/42‑2 (MQ=255)
                                        gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaatc                                <  1:3379625/64‑4 (MQ=255)
                                         aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                                     <  1:2598712/42‑1 (MQ=255)
                                         aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                                     <  1:2527603/42‑1 (MQ=255)
                                         aaaaGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAg                                               <  1:772944/48‑1 (MQ=255)
                                          aaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    >  1:2063674/1‑58 (MQ=255)
                                           aaGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCg                                                    >  1:2681255/1‑41 (MQ=255)
                                           aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGa                       >  1:1611346/1‑70 (MQ=255)
                                              cagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                    <  1:2399621/54‑1 (MQ=255)
                                              cagcCATGGACAGAAACTACCTGGCGTCAGTACTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGActct                   >  1:1312792/1‑71 (MQ=255)
                                                gcCATGGACAGAAACTACCTGGCGTCAGTGCTTGCGGGGAGCGGGCTTTCgtagt                                 <  1:3970781/55‑1 (MQ=255)
                                                gcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcc                 >  1:1153224/1‑71 (MQ=255)
                                                 ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAcccc                         >  1:1669971/1‑62 (MQ=255)
                                                 ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc                <  1:4069344/71‑1 (MQ=255)
                                                 ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc                >  1:264221/1‑71 (MQ=255)
                                                   aTGGATATAAACTATCTCGCGTCAGTCCTTGCGGGGa                                                <  1:2497878/37‑1 (MQ=25)
                                                     ggACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAg                                               <  1:3851405/36‑1 (MQ=255)
                                                     ggACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                     <  1:811416/46‑1 (MQ=255)
                                                      gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc            <  1:1583840/70‑1 (MQ=255)
                                                      gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGGGGCTTTCGTAGTAATACCCCGa                       >  1:567215/1‑58 (MQ=255)
                                                         aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtag                                  >  1:1670170/1‑45 (MQ=255)
                                                         aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTccc                >  1:4109688/1‑63 (MQ=255)
                                                         aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCtt         >  1:2025505/1‑70 (MQ=255)
                                                         aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa        <  1:1098499/71‑1 (MQ=255)
                                                          gAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAc                            <  1:3293331/50‑1 (MQ=255)
                                                          gAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa        <  1:1166993/70‑1 (MQ=255)
                                                            aaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAccc                          >  1:3564837/1‑50 (MQ=255)
                                                             aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                               >  1:3507135/1‑44 (MQ=255)
                                                             aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTcccc               <  1:24954/60‑1 (MQ=255)
                                                             aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc            >  1:2463527/1‑63 (MQ=255)
                                                             aCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAaca    <  1:1809961/71‑1 (MQ=255)
                                                              cTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAacac   <  1:42574/71‑1 (MQ=255)
                                                                 ccTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAAcaca  <  1:473885/69‑1 (MQ=255)
                                                                 |                                                                    
AATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACA  >  minE/218101‑218234

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: