Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,065,625 C→T 11.7% intergenic (+24/+55) rpoD → / ← ygjF RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,065,6250CT11.7% 87.4 / 3.2 34intergenic (+24/+55)rpoD/ygjFRNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (10/20);  new base T (2/2);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCA  >  minE/2065555‑2065693
                                                                      |                                                                    
gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAc                                                                      <  1:1513807/71‑1 (MQ=255)
gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAc                                                                      <  1:1708332/71‑1 (MQ=255)
gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGACCAGGCGTTAc                                                                      <  1:3806964/71‑1 (MQ=255)
 tCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcg                                                                     <  1:814064/71‑1 (MQ=255)
   aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgcc                                                                   >  1:1999746/1‑71 (MQ=255)
   aCCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc                                                                    <  1:734280/70‑1 (MQ=255)
     ccGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgccg                                                                  >  1:2733620/1‑70 (MQ=255)
      cGAGCCTTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa                                                                <  1:2153378/71‑1 (MQ=255)
       gAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCa                                                                >  1:2945663/1‑70 (MQ=255)
          ccGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCAGGATCAGGCGTTACGCCGCAcc                                                              <  1:520334/69‑1 (MQ=255)
            gTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAc                                                                      <  1:156026/59‑1 (MQ=255)
               cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCAc                                                        <  1:3844102/70‑1 (MQ=255)
               cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCAc                                                        <  1:3643103/70‑1 (MQ=255)
                  aaGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACt                                                       <  1:945741/68‑1 (MQ=255)
                  aaGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACt                                                       <  1:1825606/68‑1 (MQ=255)
                       ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCct                                                <  1:2598069/70‑1 (MQ=255)
                       ctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc                                               >  1:1226603/1‑71 (MQ=255)
                         gcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctct                                              >  1:3203436/1‑70 (MQ=255)
                              gCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaa                                        <  1:3166237/71‑1 (MQ=255)
                               cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcg                                                                     <  1:1345343/41‑1 (MQ=255)
                               cTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaa                                        <  1:1229329/70‑1 (MQ=255)
                                   cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa                                                                >  1:3582902/1‑42 (MQ=255)
                                   cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCa                                                                >  1:483004/1‑42 (MQ=255)
                                   cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGca                                           >  1:1682844/1‑63 (MQ=255)
                                    tGGACGATTAATCGGTAGGCCGGATCAGGCGTTAcgc                                                                    <  1:1477831/37‑1 (MQ=255)
                                      gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCAcc                                <  1:1116410/71‑1 (MQ=255)
                                      gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACCCCAcc                                <  1:1591364/71‑1 (MQ=255)
                                        cGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtt                              <  1:851146/71‑1 (MQ=255)
                                         gATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt                             <  1:623467/71‑1 (MQ=255)
                                           ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtttt                            <  1:3051385/70‑1 (MQ=255)
                                           ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCtttt                            <  1:264951/70‑1 (MQ=255)
                                             aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCAc                                 >  1:1776422/1‑63 (MQ=255)
                                                cGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc                                               <  1:3110541/46‑1 (MQ=255)
                                                   tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTg                                             <  1:504840/45‑1 (MQ=255)
                                                    aGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTg                       >  1:614746/1‑66 (MQ=255)
                                                     ggCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt                             <  1:1505723/59‑1 (MQ=255)
                                                     ggCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtt                  <  1:107689/70‑1 (MQ=255)
                                                          gATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGttttt               <  1:2218672/68‑1 (MQ=255)
                                                              aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGcc        >  1:3742418/1‑71 (MQ=255)
                                                               ggCGTTACGCCGCACCCGGCACTAGGCGCTCTGCACAAACGCCACCTTTTCgg                         <  1:1837046/53‑1 (MQ=255)
                                                                 cGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt                             <  1:2461272/47‑1 (MQ=255)
                                                                   ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACg    >  1:473691/1‑70 (MQ=255)
                                                                   ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACg    >  1:2059866/1‑70 (MQ=255)
                                                                    tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCCTTTCGGTGGCGTTTTTTATCGCCCACGCa  >  1:3598320/1‑71 (MQ=255)
                                                                      |                                                                    
GTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCA  >  minE/2065555‑2065693

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: