Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,065,028 T→C 8.8% R423R (CGT→CGC rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,065,0280TC 8.8% 330.8 / 16.6 113R423R (CGT→CGCrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (64/39);  new base C (8/2);  total (72/41)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCA  >  minE/2064961‑2065096
                                                                   |                                                                    
tCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   >  1:2204377/1‑71 (MQ=255)
 ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   >  1:4955241/1‑70 (MQ=255)
 ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   >  1:425082/1‑70 (MQ=255)
    tGACCTGATTCAGGAAGGCAAAATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTAc                                                                >  1:1820771/1‑70 (MQ=255)
      aCCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             <  1:2093895/71‑1 (MQ=255)
       ccTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             <  1:2491365/70‑1 (MQ=255)
       ccTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                                                            <  1:1108804/71‑3 (MQ=12)
         tGATCCAGGAAGGCAATATCGGCCTGATGAAGGCTGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                                                            <  1:2569484/69‑3 (MQ=9)
          gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGt                                                                      >  1:2818843/1‑58 (MQ=255)
          gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:491083/1‑66 (MQ=255)
          gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:2829745/1‑66 (MQ=255)
            ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   <  1:4595197/59‑1 (MQ=255)
            ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   <  1:1495372/59‑1 (MQ=255)
            ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                                                       <  1:948051/71‑1 (MQ=255)
            ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                                                       <  1:1625803/71‑1 (MQ=255)
             ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                                                      >  1:3170143/2‑69 (MQ=255)
             ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                                                      >  1:3173542/2‑69 (MQ=255)
             ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                                                      >  1:2428536/2‑69 (MQ=255)
                 gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                                               <  1:2026326/58‑1 (MQ=255)
                  aaGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                                               <  1:3690460/57‑1 (MQ=255)
                  aaGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTAc                                                 >  1:2511348/1‑71 (MQ=255)
                   aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCtg                                                  >  1:1743491/1‑68 (MQ=255)
                   aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                                                  >  1:1291949/1‑69 (MQ=255)
                   aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                                                  >  1:5043968/1‑69 (MQ=255)
                   aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                                                  >  1:2254085/1‑69 (MQ=255)
                   aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTATCCACCTa                                                  >  1:1578098/1‑69 (MQ=255)
                        aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc                                           <  1:3107784/71‑1 (MQ=255)
                        aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc                                           >  1:4898665/1‑71 (MQ=255)
                        aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc                                           >  1:991692/1‑71 (MQ=255)
                         taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                                                           >  1:4970356/2‑54 (MQ=255)
                         taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                                                           >  1:75893/2‑54 (MQ=255)
                         taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                                                           >  1:2656642/2‑54 (MQ=255)
                         cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt                                          <  1:4866391/71‑1 (MQ=255)
                         cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt                                          <  1:483580/71‑1 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:742241/1‑50 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:3728995/1‑50 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:2881537/1‑50 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:2774363/1‑50 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                                              >  1:2294084/1‑50 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                                               >  1:72336/1‑49 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                                                        >  1:3423525/1‑56 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                                                        >  1:3295095/1‑56 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                                                        >  1:1387133/1‑56 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc                                               >  1:1390675/1‑65 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc                                               >  1:813768/1‑65 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc                                               >  1:2513981/1‑65 (MQ=255)
                          aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATAACGCCGTGGTTACAAGTtctc                                                        >  1:5133367/1‑56 (MQ=255)
                           tCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCtgg                                        <  1:3119972/71‑1 (MQ=255)
                               tCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGa                                    <  1:2469570/71‑1 (MQ=255)
                                cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                                                   >  1:4496882/1‑55 (MQ=255)
                                cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGAt                                   >  1:1189773/1‑71 (MQ=255)
                                cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTATCCACCTACGCAACCTGGTGGAt                                   >  1:5067196/1‑71 (MQ=255)
                                 tgatgaAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTAGACCTATGCGACGTGGTGCAtt                                  >  1:4832722/1‑70 (MQ=255)
                                 tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   <  1:3878127/38‑1 (MQ=255)
                                 tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   <  1:1532944/38‑1 (MQ=255)
                                 tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                                                   <  1:2316836/38‑1 (MQ=255)
                                 tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAc                                            >  1:2689457/1‑61 (MQ=255)
                                 tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAc                                            >  1:3162790/1‑61 (MQ=255)
                                 tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc                                  >  1:3301009/1‑71 (MQ=255)
                                   atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGtt                                                           >  1:4831013/1‑44 (MQ=255)
                                   atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGt                                                            >  1:4011548/1‑43 (MQ=255)
                                   atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATcc                                 <  1:2333211/70‑1 (MQ=255)
                                   atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATcc                                 >  1:3176972/1‑70 (MQ=255)
                                   atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATcc                                 <  1:2991830/70‑1 (MQ=255)
                                   atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCg                                >  1:523219/1‑71 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:4428889/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:5114657/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:4738799/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:2045447/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:2001162/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:2190888/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:4143497/1‑40 (MQ=255)
                                     gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                                             >  1:204568/1‑40 (MQ=255)
                                      aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGt                               >  1:1267602/1‑69 (MQ=255)
                                      aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGt                               >  1:3969095/1‑69 (MQ=255)
                                       aaGCGGTTGATAAATTCGAATACCGCCGTTGTTACAAGTTCTcc                                                       <  1:1189966/44‑1 (MQ=255)
                                         gCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGc                          <  1:875616/71‑1 (MQ=255)
                                          cGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCg                         <  1:1617499/71‑1 (MQ=255)
                                             ttGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                                                   <  1:4622198/42‑1 (MQ=255)
                                             ttGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                                                   <  1:877670/42‑1 (MQ=255)
                                             ttGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGt                               <  1:2898550/62‑1 (MQ=255)
                                                aTAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCa                             >  1:1089195/1‑61 (MQ=255)
                                                aTAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCa                             >  1:616585/1‑61 (MQ=255)
                                                aTAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCa                             >  1:4811615/1‑61 (MQ=255)
                                                aTAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCa                             >  1:3553761/1‑61 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCATCTACGCAACCTGGTGGATc                                  >  1:704683/1‑49 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc                                  >  1:887585/1‑49 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc                                  >  1:855779/1‑49 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc                                  >  1:1769515/1‑49 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc                                  >  1:3770414/1‑49 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc                                  >  1:1162379/1‑49 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGctct             <  1:3495342/70‑1 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGctct             <  1:4056713/70‑1 (MQ=255)
                                                       cGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGctct             <  1:4400454/70‑1 (MQ=255)
                                                                 cttGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCa  <  1:2055846/69‑1 (MQ=255)
                                                                 cgTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGTTCTATCGCGGATCa  <  1:3021800/71‑1 (MQ=255)
                                                                 cgTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGTTCTATCGCGGATCa  <  1:3015837/71‑1 (MQ=255)
                                                                 cgTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCa  <  1:2707594/71‑1 (MQ=255)
                                                                   tGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATc   <  1:480481/68‑1 (MQ=255)
                                                                   tGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATc   <  1:1108133/68‑1 (MQ=255)
                                                                   |                                                                    
TCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCA  >  minE/2064961‑2065096

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: