Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,864,336 G→A 12.2% intergenic (+83/+65) ygcP → / ← ygcQ predicted anti‑terminator regulatory protein/predicted flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,864,3360GA12.2% 230.1 / 18.9 82intergenic (+83/+65)ygcP/ygcQpredicted anti‑terminator regulatory protein/predicted flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (45/27);  new base A (6/4);  total (51/31)
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.47e-01

ACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTT  >  minE/1864267‑1864402
                                                                     |                                                                  
aCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGa                                                                   >  1:3655392/1‑71 (MQ=255)
  cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcg                                                                 >  1:1812785/1‑71 (MQ=255)
   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                >  1:889229/1‑71 (MQ=255)
   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                <  1:3915297/71‑1 (MQ=255)
   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                >  1:183237/1‑71 (MQ=255)
    aTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgcg                                                               >  1:3412795/1‑71 (MQ=255)
     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                <  1:1264782/69‑1 (MQ=255)
     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGt                                                              <  1:3197546/71‑1 (MQ=255)
     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGt                                                              >  1:2380230/1‑71 (MQ=255)
     gtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGt                                                              >  1:675455‑M1/48‑71 (MQ=255)
      tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTc                                                             >  1:244389/1‑71 (MQ=255)
       gattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTc                                                             >  1:2399343‑M1/46‑70 (MQ=255)
       gattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTc                                                             <  1:3374261‑M1/25‑1 (MQ=255)
       gAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTc                                                             <  1:2284437/70‑1 (MQ=255)
       gAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCt                                                            >  1:3083328/1‑71 (MQ=255)
        attttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCt                                                            >  1:1847499‑M1/45‑70 (MQ=255)
        attttgtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCtt                                                           >  1:2767252‑M1/45‑71 (MQ=255)
        aaTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCtt                                                           <  1:1233038/71‑1 (MQ=255)
         aTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                >  1:736332/1‑65 (MQ=255)
         aTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCtt                                                           <  1:3691480/70‑1 (MQ=255)
         aTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTa                                                          <  1:3628077/71‑1 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGt                                                              >  1:3517401/1‑66 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTa                                                          >  1:2067448/1‑70 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                         >  1:3749855/1‑71 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                         >  1:2289927/1‑71 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                         <  1:684465/71‑1 (MQ=255)
           ttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTAt                                                         <  1:3291702‑M1/29‑1 (MQ=255)
           ttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTc                                                             >  1:786558/1‑66 (MQ=255)
             ggTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGa                                                                   <  1:266688/58‑1 (MQ=255)
             ggTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                         >  1:921810/1‑68 (MQ=255)
              taggtcggataaggcgtttacgccgcatccgacatcaaTGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                      >  1:2123991‑M1/39‑70 (MQ=255)
              gTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                      <  1:783991/70‑1 (MQ=255)
               tttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                         <  1:1006433/66‑1 (MQ=255)
                ttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                   >  1:3014361/1‑71 (MQ=255)
                ggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                   <  1:3522756‑M1/35‑1 (MQ=255)
                 tGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                   >  1:3030767/1‑70 (MQ=255)
                   cggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                                <  1:1276829‑M1/38‑1 (MQ=255)
                     gataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCtacgc                                              >  1:3671455‑M1/32‑69 (MQ=255)
                      ataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCtacgc                                              <  1:1943384‑M1/40‑3 (MQ=255)
                       gTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                      >  1:1918647/1‑61 (MQ=255)
                              tttacgccgcatccgacatcaaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                  >  1:2150327‑M1/23‑58 (MQ=255)
                              tttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTacg                                               >  1:654835‑M1/23‑60 (MQ=255)
                                tacgccgcaaccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCt                                                            <  1:993478‑M1/26‑1 (MQ=255)
                                aCGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGa                                                                   >  1:3223308/1‑39 (MQ=255)
                                   ccACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgag                                                    >  1:4143059/1‑49 (MQ=255)
                                       aTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAAc                              >  1:998385/1‑69 (MQ=255)
                                       aTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg                            >  1:1853075/1‑71 (MQ=255)
                                       aTCCGACATTCAGCCCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTccc                                         >  1:2534839/1‑58 (MQ=255)
                                        tCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCtt                                                           <  1:1604680/39‑1 (MQ=255)
                                        tCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg                            >  1:2816287/1‑70 (MQ=255)
                                        tCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg                            >  1:2074387/1‑70 (MQ=255)
                                         tccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTa                                                          >  1:1285918‑M1/12‑39 (MQ=255)
                                         tccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTa                                                          >  1:18329‑M1/12‑39 (MQ=255)
                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                         >  1:3093080/1‑40 (MQ=255)
                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGa                                <  1:1149346/65‑1 (MQ=255)
                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg                            >  1:2861536/1‑69 (MQ=255)
                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg                            >  1:3902696/1‑69 (MQ=255)
                                          cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCtt                                                           <  1:2069277/37‑1 (MQ=255)
                                          cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAg                         >  1:2500803/1‑71 (MQ=255)
                                           gACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAgg                        >  1:3570257/1‑71 (MQ=255)
                                            aCATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAgg                        >  1:4106435/1‑70 (MQ=255)
                                              aTTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAg                     <  1:429319/71‑1 (MQ=255)
                                               ttCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAgg                    <  1:3127818/71‑1 (MQ=255)
                                                tCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGa                   >  1:647977/1‑71 (MQ=255)
                                                 cAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGa                   <  1:2363599/70‑1 (MQ=255)
                                                   gcgcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatca                <  1:2094197/71‑1 (MQ=255)
                                                     gcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCa                                       >  1:1542726/1‑46 (MQ=40)
                                                        tGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAg                     <  1:3818877/61‑1 (MQ=255)
                                                        tGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAg                     <  1:1783208/61‑1 (MQ=255)
                                                          aTGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGc          >  1:958563/1‑70 (MQ=255)
                                                          aTGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGc          >  1:344903/1‑70 (MQ=255)
                                                           tGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCGCAGGCCAGAACCGTAGGCAGGatca                <  1:4072608/63‑1 (MQ=255)
                                                            gCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCGGGATCATCCGGCaa        >  1:1316916/1‑70 (MQ=255)
                                                            gCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCaa        >  1:541131/1‑70 (MQ=255)
                                                            gCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCaa        >  1:1825670/1‑70 (MQ=255)
                                                              gACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAacac     >  1:2590887/1‑71 (MQ=255)
                                                               aCGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAAcc                             <  1:2212714/46‑1 (MQ=255)
                                                               aCGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAacac     <  1:1715271/70‑1 (MQ=255)
                                                                cGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGt   <  1:2476790/71‑1 (MQ=255)
                                                                cGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGt   >  1:281553/1‑71 (MQ=255)
                                                                 gCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGtt  <  1:1837918/71‑1 (MQ=255)
                                                                 gCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGt   >  1:27855/1‑70 (MQ=255)
                                                                     |                                                                  
ACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTT  >  minE/1864267‑1864402

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: