Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 5,593 A→T 8.7% intergenic (+63/+90) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE5,5930AT 8.7% 186.3 / 9.2 69intergenic (+63/+90)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base T (2/4);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

AAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGC  >  minE/5529‑5659
                                                                |                                                                  
aaaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                              >  1:158555/1‑71 (MQ=255)
 aaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                              <  1:1861934/70‑1 (MQ=255)
 aaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAc                                                             <  1:1659446/71‑1 (MQ=255)
     aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGtt                                                          >  1:3717708/1‑70 (MQ=255)
       aaaTGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                                        <  1:979416/70‑1 (MQ=255)
         aTGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTc                                                     >  1:1329133/1‑71 (MQ=255)
                 tAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatat                                             <  1:739098/71‑1 (MQ=255)
                   cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg                                           <  1:2440899‑M1/68‑19 (MQ=255)
                   cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg                                           <  1:2281548‑M1/68‑19 (MQ=255)
                   cgcatccGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtagg                                           <  1:3492425‑M1/68‑19 (MQ=255)
                   aCAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatata                                            >  1:3280752/1‑70 (MQ=255)
                    gcatcCGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggc                                          <  1:759492‑M1/69‑20 (MQ=255)
                    cAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCa                                                 >  1:1480437/1‑64 (MQ=255)
                     catCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggcc                                         <  1:719446‑M1/70‑21 (MQ=255)
                      aTCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgatt                                                  <  1:315886‑M1/61‑12 (MQ=255)
                      aTCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattt                                                 <  1:392510‑M1/62‑13 (MQ=255)
                       tCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAc                                                             >  1:3737090/1‑49 (MQ=255)
                       tCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTTATTCTGCAATATATTGaa                                       >  1:1590653/1‑71 (MQ=255)
                       tCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatata                                            <  1:3868539/66‑1 (MQ=255)
                        ccGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccgga                                      >  1:2010629‑M1/1‑48 (MQ=255)
                          ggCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt                                    <  1:2009738/71‑1 (MQ=255)
                          ggCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                               <  1:167075/44‑1 (MQ=255)
                           gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                   >  1:1148074/1‑71 (MQ=255)
                           gCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggat                                     >  1:579630‑M1/1‑45 (MQ=255)
                             aTTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagc                                 <  1:746257‑M1/71‑29 (MQ=255)
                               tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATg                               >  1:2321021/1‑71 (MQ=255)
                                caacgcctgATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagc                                 >  1:317321‑M1/4‑40 (MQ=255)
                                cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTTATTCTGCAATATATTGAATCTGCATGc                              >  1:772863/1‑71 (MQ=255)
                                cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGc                              >  1:1961496/1‑71 (MQ=255)
                                 aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt                             >  1:3159980/1‑71 (MQ=255)
                                  gcgcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGc                                  <  1:2370522/65‑1 (MQ=255)
                                    gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt                          >  1:1413923/1‑71 (MQ=255)
                                     ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt                          <  1:1291823/70‑1 (MQ=255)
                                         aTGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatc                       <  1:11235‑M1/69‑39 (MQ=255)
                                         aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTc                                                     <  1:1888829/39‑1 (MQ=255)
                                         aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                     <  1:2561009/71‑1 (MQ=255)
                                         aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg                      <  1:376284/70‑1 (MQ=255)
                                         aTGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccg                     >  1:1996752‑M1/1‑31 (MQ=255)
                                          tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                     >  1:3251649/1‑70 (MQ=255)
                                             gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                     <  1:2569631/67‑1 (MQ=255)
                                             gACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagc                                 <  1:1042922‑M1/55‑29 (MQ=255)
                                             gACGCTGACGCGTCTTATCTGGCCTACacgctgcgatt                                                  >  1:1632913‑M1/1‑27 (MQ=255)
                                               cGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacg                >  1:1249894‑M1/1‑25 (MQ=255)
                                               cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa                >  1:3668316/1‑70 (MQ=255)
                                               cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt               <  1:2723348/71‑1 (MQ=255)
                                                gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt               >  1:3482613/1‑70 (MQ=255)
                                                gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa              <  1:1389208/71‑1 (MQ=255)
                                                   ggCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa                <  1:841959/66‑1 (MQ=255)
                                                    gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGc          >  1:3199098/1‑71 (MQ=255)
                                                      gcgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGt        >  1:3160887/1‑71 (MQ=255)
                                                       cgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt       <  1:1265795/71‑1 (MQ=255)
                                                         tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                   <  1:1395336/41‑1 (MQ=255)
                                                         tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt       <  1:2891438/69‑1 (MQ=255)
                                                          cTTATCAGGCCTACGTTAATTCTGCAATATATTGAAt                                      <  1:785535/37‑1 (MQ=255)
                                                          cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa     <  1:1661665/70‑1 (MQ=255)
                                                          cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc    >  1:2892450/1‑71 (MQ=255)
                                                           ttATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcat                        <  1:2102570‑M1/50‑38 (MQ=255)
                                                           ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                   <  1:3158292/39‑1 (MQ=255)
                                                           ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                   <  1:1174162/39‑1 (MQ=255)
                                                            tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc  >  1:2386255/1‑71 (MQ=255)
                                                            tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc  >  1:2381545/1‑71 (MQ=255)
                                                             aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg           <  1:3352438/61‑1 (MQ=255)
                                                             aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg           <  1:3764262/61‑1 (MQ=255)
                                                             aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg   <  1:731706/69‑1 (MQ=255)
                                                             aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg   <  1:2396046/69‑1 (MQ=255)
                                                             aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg   <  1:1736853/69‑1 (MQ=255)
                                                             aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc    <  1:3514893/68‑1 (MQ=255)
                                                              tCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgca                         <  1:165201‑M1/46‑37 (MQ=255)
                                                               cAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg         >  1:479535/1‑61 (MQ=255)
                                                                |                                                                  
AAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGC  >  minE/5529‑5659

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: