Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,245,573 T→C 17.0% intergenic (‑89/‑6) nudC ← / → rsd NADH pyrophosphatase/stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,245,5730TC17.0% 114.3 / 17.7 53intergenic (‑89/‑6)nudC/rsdNADH pyrophosphatase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (4/5);  total (26/27)
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.25e-01

TCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCGCT  >  minE/2245506‑2245642
                                                                   |                                                                     
tCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGCCaa                                                                    <  1:2379956/71‑1 (MQ=255)
 cACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                                                                   >  1:48828/1‑71 (MQ=255)
 cACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                                                                   >  1:2496126/1‑71 (MQ=255)
  aCTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                                                                   >  1:350957/1‑70 (MQ=255)
  aCTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAAt                                                                   <  1:2942731/70‑1 (MQ=255)
  aCTGGCATGTTTAACATGGTTTTCACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATc                                                                  <  1:2384315/71‑1 (MQ=255)
    tGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCa                                                                 <  1:418021/70‑1 (MQ=255)
         tGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                                                           <  1:938289/71‑1 (MQ=255)
          gTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTaa                                                          >  1:2457878/1‑71 (MQ=255)
           tttAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCt                                                             <  1:1542080/67‑1 (MQ=255)
             tAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAAc                                                         >  1:61373/1‑69 (MQ=255)
               aCATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGt                                                                       >  1:1432585/1‑53 (MQ=255)
               aCATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGt                                                                       >  1:3475554/1‑53 (MQ=255)
                    gTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGAt                                                <  1:2914386/71‑1 (MQ=255)
                        ttACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAAcc                                            >  1:481483/1‑71 (MQ=255)
                        ttACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAAc                                             <  1:3724701/70‑1 (MQ=255)
                          aCATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt                                           >  1:827236/1‑70 (MQ=255)
                          aCATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCt                                           >  1:2220278/1‑70 (MQ=255)
                            aTTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGa                                         <  1:3001960/70‑1 (MQ=255)
                             tttCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACg                                       <  1:578601/71‑1 (MQ=255)
                               tctcACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCGCGATAACCTGACgg                                      <  1:1905517/70‑1 (MQ=255)
                                cacACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGa                                     <  1:1807284/68‑1 (MQ=255)
                                   aCTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTa                                                           <  1:2943362/45‑1 (MQ=255)
                                   aCTGAGCAGTTTTTGAATACAAACTTGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGAc                                        >  1:2292706/1‑64 (MQ=255)
                                   aCTGAGCAGTTTTTGAATACAAACTTGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGACGGAg                                    <  1:660035/68‑1 (MQ=255)
                                   aCTGAGCAGTATTTGAATACAAACTTGCGGAGTCAATCAAGCTAAACCAGCTCGATAACCTGACGGAgcgc                                 >  1:1529093/1‑71 (MQ=255)
                                       aGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAAc                                                         >  1:2123245/1‑43 (MQ=255)
                                       aGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAAc                                                         >  1:1760548/1‑43 (MQ=255)
                                       aGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGt                               >  1:2879229/1‑69 (MQ=255)
                                       aGCAGTTTTTGAATACAAACTTGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGt                               >  1:1051486/1‑69 (MQ=255)
                                       aGCAGTTTTTGAATACAAACTTGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTc                              <  1:757182/70‑1 (MQ=255)
                                            tttttGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACg                                       <  1:1938659/56‑1 (MQ=255)
                                            tttttGAATACAAACTTGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAgg                         >  1:486991/1‑70 (MQ=255)
                                             ttttGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATa                                               >  1:2161487/1‑47 (MQ=255)
                                             ttttGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCaga                           >  1:3087617/1‑67 (MQ=255)
                                              tttGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAg                      >  1:3184741/1‑71 (MQ=255)
                                                tGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAg                      <  1:1429551/69‑1 (MQ=255)
                                                  aaTACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCa                             <  1:3385285/60‑1 (MQ=255)
                                                      cAAACTTGCGGAGTCAATCATGCTTAACCAGCACGATAACCTGACgg                                      >  1:3611803/1‑47 (MQ=255)
                                                       aaaCTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAgc                                   >  1:2355838/1‑49 (MQ=255)
                                                       aaaCTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACt             >  1:1519888/1‑71 (MQ=255)
                                                        aaCTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTg            >  1:3908832/1‑71 (MQ=255)
                                                         aCTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTgg           <  1:2154754/71‑1 (MQ=255)
                                                         aCTTGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTgg           >  1:3776596/1‑71 (MQ=255)
                                                           ttGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGtt         >  1:2602607/1‑71 (MQ=255)
                                                           ttGCGGAGCCAATCATGCTTAACCAGCTCGATAACCTGa                                         <  1:2260759/39‑1 (MQ=255)
                                                              cGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGa       <  1:1376199/70‑1 (MQ=255)
                                                                gAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATc     >  1:439155/1‑70 (MQ=255)
                                                                gAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATc     <  1:1654841/70‑1 (MQ=255)
                                                                gAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCg    <  1:1152579/71‑1 (MQ=255)
                                                                gAGCCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCg    <  1:2466398/71‑1 (MQ=255)
                                                                 aGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCGc   <  1:2684914/71‑1 (MQ=255)
                                                                 aGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCGc   >  1:2536967/1‑71 (MQ=255)
                                                                 aGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCGc   <  1:1362204/71‑1 (MQ=255)
                                                                  gTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCGCt  <  1:231531/71‑1 (MQ=255)
                                                                   tCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCg    <  1:1896448/68‑1 (MQ=255)
                                                                   |                                                                     
TCACTGGCATGTTTAACATGGTTTTTACATTTCTCACTGAGCAGTTTTTGAATACAAACTTGCGGAGTCAATCATGCTTAACCAGCTCGATAACCTGACGGAGCGCGTCAGAGGAAGTAACAAACTGGTTGATCGCT  >  minE/2245506‑2245642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: