Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 974,990 A→T 54.6% 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,9900AT54.6% 24.8 / 103.5 77intergenic (+742/+250)lsrB/ydeHAI2 transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/12);  new base T (17/25);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 3.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCG  >  minE/974922‑975060
                                                                    |                                                                      
aCTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACt                                                                      >  1:3358226/1‑71 (MQ=255)
 cTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCtt                                                                     >  1:3711763/1‑71 (MQ=255)
 cTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCtt                                                                     >  1:324801/1‑71 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCtt                                                                     <  1:1190015/68‑1 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCtt                                                                     <  1:2817536/68‑1 (MQ=255)
    gCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGc                                                                  >  1:1965568/1‑71 (MQ=255)
      aaaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGccc                                                                >  1:1077307/1‑71 (MQ=255)
      aaaTTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGccc                                                                >  1:3193967/1‑71 (MQ=255)
         ttCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCg                                                             <  1:1858115/71‑1 (MQ=255)
           cTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCg                                                             >  1:394305/1‑69 (MQ=255)
             gACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGc                                                                  >  1:233961/1‑62 (MQ=255)
             gACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAAAAATGATAAATACTTCGc                                                                  >  1:419501/1‑62 (MQ=255)
              acacCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGc                                                         <  1:2477036/70‑1 (MQ=255)
              acacCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGc                                                         <  1:1860882/70‑1 (MQ=255)
                  cTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATaa                                                    <  1:1699199/71‑1 (MQ=255)
                    gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAg                                                          >  1:1820006/1‑63 (MQ=255)
                    gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCg                                                        <  1:3706898/65‑1 (MQ=255)
                    gCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAAc                                                  >  1:3427146/1‑71 (MQ=255)
                     cACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAg                                                          <  1:2172868/62‑1 (MQ=255)
                        gACATGCTTCTTCATCATTAGCCGCTTTGACAATAATTATAAATTCTTCGCCCCCGTAGCGATAAACCGtt                                              >  1:1182288/1‑71 (MQ=255)
                        gACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGtt                                              <  1:2253071/71‑1 (MQ=255)
                          cATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTc                                            <  1:3830649/71‑1 (MQ=255)
                           aTGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGtt                                              <  1:175568/68‑1 (MQ=255)
                             gCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGt                                          >  1:169892/1‑70 (MQ=255)
                              cttcttCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGc                                                                  >  1:3675220/1‑45 (MQ=255)
                               ttcttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTaa                                        >  1:2630398/1‑70 (MQ=255)
                                 cttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATa                                                     <  1:3923167/55‑1 (MQ=255)
                                 cttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAAc                                                  >  1:3165401/1‑58 (MQ=255)
                                 cttcATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATc                                      <  1:503374/70‑1 (MQ=255)
                                  ttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                                    >  1:2607807/1‑71 (MQ=255)
                                  ttcATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                                    >  1:2392169/1‑71 (MQ=255)
                                    catcatTAGCCGCTTTGACAATAATGATAACTTCTTCGccc                                                                <  1:1613492/41‑1 (MQ=38)
                                     atcatTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAcgc                                  <  1:31204/70‑1 (MQ=255)
                                        atTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTc                               >  1:2270606/1‑70 (MQ=255)
                                           aGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCAc                                    >  1:1768521/1‑62 (MQ=255)
                                           aGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGtt                                              >  1:84144/1‑52 (MQ=255)
                                            gCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                          <  1:3215864/71‑1 (MQ=255)
                                              cGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTgg                        <  1:2119203/71‑1 (MQ=255)
                                               gCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGc                       >  1:2858394/1‑71 (MQ=255)
                                                cgatGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCt                      <  1:3045598/68‑1 (MQ=255)
                                                 tttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGt                                >  1:3944282/1‑60 (MQ=255)
                                                 tttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGt                                >  1:453000/1‑60 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGc                       >  1:1330057/1‑68 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCt                      >  1:1115448/1‑69 (MQ=255)
                                                  ttGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCt                      >  1:565841/1‑69 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag                   >  1:1425173/1‑71 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaag                   >  1:2420000/1‑71 (MQ=255)
                                                   tGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa                    <  1:1665392/70‑1 (MQ=255)
                                                   tGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa                    <  1:2164395/70‑1 (MQ=255)
                                                    gACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa                    >  1:3298989/1‑69 (MQ=255)
                                                     aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagt                  <  1:459616/70‑1 (MQ=255)
                                                     aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagt                  <  1:2245837/70‑1 (MQ=255)
                                                     aCAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaa                    <  1:878866/68‑1 (MQ=255)
                                                     aCAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagta                 <  1:3771611/71‑1 (MQ=255)
                                                      cAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaa                >  1:1665382/1‑71 (MQ=255)
                                                       aataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaag               <  1:334920/71‑1 (MQ=255)
                                                       aataatGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaag               <  1:3439713/71‑1 (MQ=255)
                                                        ataatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTaagtaagt              <  1:3779200/71‑1 (MQ=255)
                                                         taatGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTaagtaagt              <  1:3377167/70‑1 (MQ=255)
                                                          aatGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg            >  1:1443756/1‑71 (MQ=255)
                                                          aatGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTg            >  1:3909025/1‑71 (MQ=255)
                                                            tGATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGc           <  1:185613/70‑1 (MQ=255)
                                                             gATAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCta                     >  1:3763870/1‑59 (MQ=255)
                                                             gATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagtaa                >  1:3557492/1‑64 (MQ=255)
                                                              aTAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACt                          >  1:3546535/1‑53 (MQ=255)
                                                              aTAAATTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCa         <  1:1064024/70‑1 (MQ=255)
                                                              aTAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCa         <  1:2412748/70‑1 (MQ=255)
                                                              aTAAATACTTCGCCACCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCa         <  1:3907389/70‑1 (MQ=255)
                                                                aaaTTCTTCGCCGCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCtaagt                  <  1:2241205/59‑1 (MQ=255)
                                                                aaaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                          <  1:727750/35‑1 (MQ=255)
                                                                aaaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                          <  1:1340235/35‑1 (MQ=255)
                                                                aaaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGt                                          <  1:1344741/35‑1 (MQ=255)
                                                                aaaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTaagtaagt              >  1:3577168/1‑63 (MQ=255)
                                                                 aaTTCTTCGCCCCCGTAGCGATAAACCGTTTCGTAAt                                       <  1:1813627/37‑1 (MQ=255)
                                                                  aTACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTg    >  1:1369654/1‑71 (MQ=255)
                                                                  aTACTTCGCCCCCGTAGCGATAAACCGTTTAGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTg    >  1:3050969/1‑71 (MQ=255)
                                                                   tACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGc   >  1:2223315/1‑71 (MQ=255)
                                                                    aCTTCGCCCCCGTAGCGATAAACCGTTTCGTAATAACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCg  >  1:914351/1‑71 (MQ=255)
                                                                    |                                                                      
ACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATACTTCGCCCCCGTAGCGATAAACCGTTTCGTAATCACGCGTCCAACTGGCTAAGTAAGTTGCCAGGGTGCG  >  minE/974922‑975060

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: