Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,690 C→T 18.6% A47A (GCC→GCT yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,6900CT18.6% 124.3 / 17.1 80A47A (GCC→GCTyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/39);  new base T (7/8);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGA  >  minE/1781636‑1781751
                                                      |                                                             
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:3330570/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1926195/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1233766/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:73585/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:2515442/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:2527170/69‑1 (MQ=255)
aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1447130/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:3415152/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1915364/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:62715/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:790044/68‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:1425186/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:141617/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:1568876/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:720225/69‑1 (MQ=255)
 gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:1150781/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGTCACTAttgtt                                          <  1:2981319/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:142868/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:2577819/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:3250214/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:795239/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:1908303/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:2844175/69‑1 (MQ=255)
       cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:1426249/69‑1 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGc                                                   >  1:2427475/1‑56 (MQ=37)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCggg                                       >  1:2574209/1‑68 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:1077078/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:358299/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:2081009/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:562729/1‑69 (MQ=255)
           tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:1910621/1‑69 (MQ=255)
            aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                     >  1:206445/1‑69 (MQ=255)
            aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                     >  1:1487992/1‑69 (MQ=255)
            aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                     >  1:2650515/1‑69 (MQ=255)
                      ttAATCGGCAATAACCTGGTCAATATTCTTGCCTc                                                             >  1:615845/1‑35 (MQ=255)
                      ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc                           >  1:2387500/1‑69 (MQ=255)
                      ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc                           >  1:1378778/1‑69 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:3033260/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:635200/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:675052/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:964196/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:2835426/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:194329/68‑1 (MQ=255)
                          tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:1895089/68‑1 (MQ=255)
                           cGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:3029205/67‑1 (MQ=255)
                            ggCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGAt                     >  1:2041821/1‑69 (MQ=255)
                                  aaCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGc               >  1:1583131/1‑69 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:2692847/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:552721/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:505960/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:419064/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2233189/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:187804/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:12834/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:3083514/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2505300/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:2963043/1‑67 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2883012/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2846702/67‑1 (MQ=255)
                                   aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:1066511/67‑1 (MQ=255)
                                       ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc          <  1:2381755/69‑1 (MQ=255)
                                       ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc          <  1:2593498/69‑1 (MQ=255)
                                          cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:2935050/69‑1 (MQ=255)
                                          cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:763862/69‑1 (MQ=255)
                                          cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:1239142/69‑1 (MQ=255)
                                            atatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:2522312/67‑1 (MQ=255)
                                             tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAtt      >  1:2963841/1‑67 (MQ=255)
                                             tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTg     >  1:2494517/1‑68 (MQ=255)
                                             tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATCGc    >  1:539002/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCg              >  1:2816003/1‑58 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc    >  1:2087856/1‑68 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc    >  1:2399167/1‑68 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc    >  1:2820023/1‑68 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:477606/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:3341696/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:187557/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:2841937/1‑69 (MQ=255)
                                              atTCTTGCCTCCGCGCTCGGCACTATTGGTGGGATGCGttt                               >  1:2824633/1‑41 (MQ=255)
                                               ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa  >  1:3246407/1‑69 (MQ=255)
                                               ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa  >  1:1280030/1‑69 (MQ=255)
                                                      |                                                             
AGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGA  >  minE/1781636‑1781751

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: