Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,864,337 A→C 7.6% intergenic (+84/+64) 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,3370AC 7.6% 179.9 / 12.9 79intergenic (+84/+64)ygcP/ygcQpredicted anti‑terminator regulatory protein/predicted flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (45/28);  new base C (4/2);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

ACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGG  >  minE/1864267‑1864407
                                                                      |                                                                      
aCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGa                                                                        >  1:3655392/1‑71 (MQ=255)
  cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcg                                                                      >  1:1812785/1‑71 (MQ=255)
   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                     >  1:183237/1‑71 (MQ=255)
   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                     <  1:3915297/71‑1 (MQ=255)
   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                     >  1:889229/1‑71 (MQ=255)
    aTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgcg                                                                    >  1:3412795/1‑71 (MQ=255)
     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                                     <  1:1264782/69‑1 (MQ=255)
     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGt                                                                   >  1:2380230/1‑71 (MQ=255)
     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGt                                                                   <  1:3197546/71‑1 (MQ=255)
     gtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGt                                                                   >  1:675455‑M1/48‑71 (MQ=255)
      tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTc                                                                  >  1:244389/1‑71 (MQ=255)
       gattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTc                                                                  <  1:3374261‑M1/25‑1 (MQ=255)
       gattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTc                                                                  >  1:2399343‑M1/46‑70 (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:684465/71‑1 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                              >  1:3749855/1‑71 (MQ=255)
          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                                              >  1:2289927/1‑71 (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:18329‑M1/12‑39 (MQ=255)
                                         tccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTa                                                               >  1:1285918‑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:281553/1‑71 (MQ=255)
                                                                cGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGt        <  1:2476790/71‑1 (MQ=255)
                                                                 gCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGtt       <  1:1837918/71‑1 (MQ=255)
                                                                 gCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGt        >  1:27855/1‑70 (MQ=255)
                                                                      aCGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTgg  <  1:173742/71‑1 (MQ=255)
                                                                      |                                                                      
ACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGG  >  minE/1864267‑1864407

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: