Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 974,990 A→T 33.9% 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,9900AT33.9% 63.0 / 54.6 56intergenic (+742/+250)lsrB/ydeHAI2 transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/21);  new base T (9/10);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGT  >  minE/974922‑975057
                                                                    |                                                                   
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACt                                                                   <  1:1322411/71‑1 (MQ=255)
  tGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAAttcttc                                                                 >  1:92376/1‑71 (MQ=255)
  tGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTATGACAATAATGATAAATACTTc                                                                 >  1:4133488/1‑71 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGc                                                               <  1:680050/71‑1 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGc                                                               >  1:2949272/1‑71 (MQ=255)
       aaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGcc                                                              >  1:3121587/1‑69 (MQ=255)
              acacCTGCACGCCATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCg                                                          >  1:1465120/1‑66 (MQ=255)
              acacCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCg                                                          >  1:1469923/1‑66 (MQ=255)
                acCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGAt                                                   <  1:1126570/71‑1 (MQ=255)
                  cTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATa                                                  >  1:3605117/1‑70 (MQ=255)
                   tGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCg                                                                <  1:2271201/55‑1 (MQ=255)
                     cACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAAcc                                              >  1:3818603/1‑71 (MQ=255)
                     cACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAAcc                                              <  1:2791526/71‑1 (MQ=255)
                     cACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAAcc                                              >  1:905446/1‑71 (MQ=255)
                       cGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGa                                                    >  1:305697/1‑63 (MQ=255)
                       cGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGt                                            <  1:439845/71‑1 (MQ=255)
                       cGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGt                                            <  1:2717519/71‑1 (MQ=255)
                        gACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGtt                                           >  1:3497132/1‑71 (MQ=255)
                        gACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGtt                                           <  1:2681361/71‑1 (MQ=255)
                          cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGAt                                                   <  1:1740417/61‑1 (MQ=255)
                          cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTc                                         <  1:1795957/71‑1 (MQ=255)
                            tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCg                                                     <  1:3804669/57‑1 (MQ=255)
                            tGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGt                                       <  1:993138/71‑1 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                     >  1:714515/1‑71 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                     >  1:582925/1‑71 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                     >  1:3887948/1‑71 (MQ=255)
                               ttcttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTc                                         <  1:725386/66‑1 (MQ=255)
                               ttcttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTc                                         <  1:2371065/66‑1 (MQ=255)
                                 cttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                       <  1:2903284/66‑1 (MQ=255)
                                 cttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                    <  1:3862448/69‑1 (MQ=255)
                                 cttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATc                                   <  1:2062768/70‑1 (MQ=255)
                                   tcatcaTTAGCCGCTTTGACAATAATGATAAATTCt                                                                   <  1:88530/36‑1 (MQ=255)
                                   tcatcaTTAGCCGCATTGACAATAATGATAAATACTTCGc                                                               >  1:2767199/1‑40 (MQ=255)
                                    catcatTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAcgc                               >  1:62618/1‑71 (MQ=255)
                                     atcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAcgcg                              >  1:1664071/1‑71 (MQ=255)
                                     atcatTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAcgcg                              >  1:3092536/1‑71 (MQ=255)
                                      tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCg                                             <  1:1764931/55‑1 (MQ=255)
                                      tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                                 <  1:728316/67‑1 (MQ=255)
                                      tcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                                 <  1:3492683/67‑1 (MQ=255)
                                        atTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTcc                           <  1:3913614/71‑1 (MQ=255)
                                         ttAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTcc                           >  1:2361656/1‑70 (MQ=255)
                                         ttAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                          <  1:2184743/71‑1 (MQ=255)
                                           aGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                     >  1:3890376/1‑58 (MQ=255)
                                                    gACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                    <  1:1633997/50‑1 (MQ=255)
                                                       aataatGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaag            <  1:1184634/71‑1 (MQ=255)
                                                        ataatGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCa                          <  1:3706022/56‑1 (MQ=255)
                                                         taatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCt                   >  1:2084968/1‑62 (MQ=255)
                                                          aatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCaa                         >  1:3947266/1‑55 (MQ=255)
                                                          aatGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAAc                        >  1:761296/1‑56 (MQ=255)
                                                           atGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                     >  1:3318352/1‑42 (MQ=255)
                                                              aTAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAg     <  1:2689768/71‑1 (MQ=255)
                                                               tAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAgg    <  1:991727/71‑1 (MQ=255)
                                                                aaaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAcg                                >  1:1653425/1‑42 (MQ=255)
                                                                aaaTACTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                     >  1:3982209/1‑37 (MQ=255)
                                                                 aaTACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg         <  1:1745302/64‑1 (MQ=255)
                                                                 aaTACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGt  <  1:1084051/71‑1 (MQ=255)
                                                                 aaTACTTCGCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg         <  1:1690545/63‑1 (MQ=255)
                                                                    |                                                                   
ACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGT  >  minE/974922‑975057

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: