Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,221,569 C→T 44.1% T48T (ACC→ACT yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,221,5690CT44.1% 16.7 / 51.2 68T48T (ACC→ACTyebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/25);  new base T (14/16);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GCTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGCG  >  minE/1221515‑1221636
                                                      |                                                                   
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:2145361/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:11343/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:686042/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:232208/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:1785961/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:2749440/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:1541644/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:343060/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:2888328/69‑1 (MQ=255)
gcTAATGATTGCCAGTTGCCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACtt                                                       <  1:655058/69‑1 (MQ=255)
   aaTGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTAt                                                    <  1:1635145/69‑1 (MQ=255)
   aaTGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTAt                                                    <  1:536672/69‑1 (MQ=255)
         tGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCTGGGTAATACCGTCACCATCGACTTTATGTcgg                                               <  1:906733/68‑1 (MQ=255)
         tGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTcgg                                               <  1:1484253/68‑1 (MQ=255)
         tGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTcgg                                               <  1:1381476/68‑1 (MQ=255)
         tGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTcgg                                               <  1:2621985/68‑1 (MQ=255)
         tGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTcgg                                               <  1:556834/68‑1 (MQ=255)
         tGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTcgg                                               <  1:1184318/68‑1 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGt                                                                   >  1:562028/3‑45 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:363911/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:106520/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:2543859/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:2500254/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:2462945/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:2615558/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:2078928/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:1167107/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:1462656/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:1517151/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:1617556/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:1201516/3‑69 (MQ=255)
            tggTTGGCTGATTTGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGa                                           >  1:1647306/3‑69 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:2822596/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:189609/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:3387558/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:1446983/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:2835460/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:122710/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:2682592/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:620326/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:572178/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:2194391/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:1689325/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:1747270/69‑1 (MQ=255)
                         tGGGACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:242004/69‑1 (MQ=255)
                          gggACTGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:411168/68‑1 (MQ=255)
                          gggACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:1887445/68‑1 (MQ=255)
                          gggACAGTTATCAGGACCGGGGCAATACTGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCa                              <  1:1919746/68‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:479463/69‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:3112778/69‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:283803/69‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:1075113/69‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:1891539/69‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:1563724/69‑1 (MQ=255)
                                                 aaTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAg      <  1:3080317/69‑1 (MQ=255)
                                                  aTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAgg     >  1:778263/1‑69 (MQ=255)
                                                  aTACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAgg     >  1:3304233/1‑69 (MQ=255)
                                                    aCCGTCACCTTCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:72790/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTa          >  1:136821/1‑62 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAggg  >  1:1896011/1‑68 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:36640/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:799184/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:3202548/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:516737/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:1619657/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:564146/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:314249/1‑69 (MQ=255)
                                                    aCCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:726661/1‑69 (MQ=255)
                                                     ccGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:3271301/1‑68 (MQ=255)
                                                     ccGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGc   >  1:396618/1‑68 (MQ=255)
                                                     ccGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGCg  >  1:116913/1‑69 (MQ=255)
                                                     ccGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGCg  >  1:3400993/1‑69 (MQ=255)
                                                     ccGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGCg  >  1:2889666/1‑69 (MQ=255)
                                                     ccGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGCg  >  1:1903233/1‑69 (MQ=255)
                                                      |                                                                   
GCTAATGATTGCCAGTTGGCTGATTTGGGACAGTTATCAGGACCGGGGTAATACCGTCACCATCGACTTTATGTCGGCGGATGGTATTGTTCCAGGCCGTACGCCTGTTCGTTATCAGGGCG  >  minE/1221515‑1221636

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: