Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 974,990 A→T 70.7% intergenic (+742/+250) lsrB → / ← ydeH AI2 transporter/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE974,9900AT70.7% 111.1 / 59.7 75intergenic (+742/+250)lsrB/ydeHAI2 transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/12);  new base T (24/29);  total (34/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 = 9.28e-01

ACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCA  >  minE/974922‑975053
                                                                    |                                                               
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCt                                                               >  1:454050/1‑71 (MQ=255)
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCt                                                               >  1:3118925/1‑71 (MQ=255)
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCt                                                               >  1:1597239/1‑71 (MQ=255)
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCt                                                               <  1:3168915/71‑1 (MQ=255)
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATAc                                                                <  1:151066/70‑1 (MQ=255)
   ggCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCg                                                            >  1:568763/1‑71 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAAttcttc                                                             <  1:5250062/69‑1 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGc                                                           >  1:3808385/1‑71 (MQ=255)
     cAAATTCTGACACCTGCACGACATGCTTCTTCATCTTTAGCCGCTTTGACAATAATGATAAATACt                                                               <  1:2651586/66‑1 (MQ=255)
     cAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACt                                                               <  1:4625247/66‑1 (MQ=255)
      aaaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGccc                                                         >  1:1569809/1‑71 (MQ=255)
      aaaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGccc                                                         >  1:1222844/1‑71 (MQ=255)
      aaaTTCTGACACCTGCACGACATGCTTCTTCATAATTAGCCGCTTTGACAATAATGATAAATACTTCGccc                                                         >  1:720883/1‑71 (MQ=255)
       aaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGcccc                                                        <  1:4624775/71‑1 (MQ=255)
        aTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGccccc                                                       >  1:5488270/1‑71 (MQ=255)
           cTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCg                                                      >  1:3686445/1‑69 (MQ=255)
           cTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGt                                                     <  1:3146865/70‑1 (MQ=255)
              acacCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGc                                                  <  1:318360/70‑1 (MQ=255)
                    gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATaaa                                            >  1:3162256/1‑70 (MQ=255)
                    gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATa                                              >  1:1129987/1‑68 (MQ=255)
                    gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAAc                                           >  1:5487342/1‑71 (MQ=255)
                       cGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGt                                        <  1:4149173/71‑1 (MQ=255)
                            tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGc                                                  <  1:2858456/56‑1 (MQ=255)
                            tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   <  1:3887291/71‑1 (MQ=255)
                            tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   <  1:5277042/71‑1 (MQ=255)
                            tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   <  1:2073717/71‑1 (MQ=255)
                             gCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   >  1:4779493/1‑70 (MQ=255)
                             gCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   >  1:4605926/1‑70 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGc                                                           >  1:4189037/1‑45 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAAc                                           >  1:5428693/1‑61 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTAGTaa                                 >  1:3637422/1‑71 (MQ=255)
                                   tcatcaTTAGCCGCTTTGACAATAATGATAAATTCTTCg                                                            >  1:2670062/1‑39 (MQ=255)
                                    catcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                             <  1:3052720/69‑1 (MQ=255)
                                    catcatTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                             <  1:3758576/69‑1 (MQ=255)
                                     atcatTAGCCGCTTTGACAATAATGATAAAttcttc                                                             >  1:1277513/1‑36 (MQ=255)
                                      tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATa                                              <  1:1442424/50‑1 (MQ=255)
                                      tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                <  1:2058439/64‑1 (MQ=255)
                                       catTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                        >  1:485284/1‑71 (MQ=255)
                                       catTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                        >  1:5678461/1‑71 (MQ=255)
                                       catTAGCAGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                        >  1:4622707/1‑71 (MQ=255)
                                        atTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                        >  1:1568312/1‑70 (MQ=255)
                                            gCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                >  1:2234746/1‑58 (MQ=255)
                                            gCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAAc                    <  1:4493710/70‑1 (MQ=255)
                                            gCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                   >  1:5143582/1‑71 (MQ=255)
                                            gCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                   <  1:5494066/71‑1 (MQ=255)
                                               gCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAAc                                           >  1:936965/1‑44 (MQ=255)
                                               gCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                <  1:1672963/55‑1 (MQ=255)
                                                cTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATATACCGtt                                       <  1:5868/47‑1 (MQ=25)
                                                cTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGtt                                       <  1:1411235/47‑1 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCa                              >  1:1482512/1‑54 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                      <  1:673418/62‑1 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                      <  1:59323/62‑1 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCt               >  1:2465626/1‑69 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag            >  1:2195049/1‑71 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa             <  1:2199451/70‑1 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa             <  1:1701111/70‑1 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa             <  1:3104732/70‑1 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCta              <  1:1848892/69‑1 (MQ=255)
                                                    gACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTc                                     >  1:5165567/1‑45 (MQ=255)
                                                     aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagt           <  1:2545339/70‑1 (MQ=255)
                                                     aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagt           <  1:4878970/70‑1 (MQ=255)
                                                      cAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagta          <  1:2822703/70‑1 (MQ=255)
                                                      cAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                   <  1:2154759/61‑1 (MQ=255)
                                                      cAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                   <  1:4045593/61‑1 (MQ=255)
                                                      cAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                   <  1:3427124/61‑1 (MQ=255)
                                                       aataatGATAAATTCTTCGCCCCCGTAGCGATAAACCg                                         >  1:3568427/1‑38 (MQ=255)
                                                        ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   <  1:3645971/43‑1 (MQ=255)
                                                        ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                 <  1:5562736/45‑1 (MQ=255)
                                                        ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCTTCCAACTGGCtaagtaag        <  1:1123286/70‑1 (MQ=255)
                                                           atGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGc    <  1:2569442/71‑1 (MQ=255)
                                                            tGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                   <  1:3190935/39‑1 (MQ=255)
                                                            tGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                <  1:2826843/42‑1 (MQ=255)
                                                            tGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGcc   <  1:2755606/71‑1 (MQ=255)
                                                              aTAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCa  <  1:1998514/70‑1 (MQ=255)
                                                               tAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag            >  1:5060459/1‑59 (MQ=255)
                                                               tAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag            >  1:4617259/1‑59 (MQ=255)
                                                               tAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag            >  1:3051362/1‑59 (MQ=255)
                                                                    |                                                               
ACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCA  >  minE/974922‑975053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: