Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,684 T→C 19.1% I45I (ATT→ATC yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,6840TC19.1% 117.6 / 17.2 78I45I (ATT→ATCyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (24/39);  new base C (7/8);  total (31/47)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGA  >  minE/1781616‑1781751
                                                                    |                                                                   
gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                                     <  1:2413440/69‑1 (MQ=255)
gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                                     <  1:1298126/69‑1 (MQ=255)
gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                                     <  1:1446769/69‑1 (MQ=255)
gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                                     <  1:1403515/69‑1 (MQ=255)
gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                                     <  1:328320/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:2515442/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1233766/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1447130/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:73585/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:2527170/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:1926195/69‑1 (MQ=255)
                    aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                                 <  1:3330570/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:1568876/69‑1 (MQ=255)
                     gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:1150781/69‑1 (MQ=255)
                     gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:720225/69‑1 (MQ=255)
                     gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:141617/69‑1 (MQ=255)
                     gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                                <  1:1425186/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGTCACTAttgtt                                          <  1:2981319/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:795239/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:2577819/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:2844175/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:3250214/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:1908303/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:1426249/69‑1 (MQ=255)
                           cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                                          <  1:142868/69‑1 (MQ=255)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGc                                                   >  1:2427475/1‑56 (MQ=37)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCggg                                       >  1:2574209/1‑68 (MQ=255)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:562729/1‑69 (MQ=255)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:358299/1‑69 (MQ=255)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:1077078/1‑69 (MQ=255)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:2081009/1‑69 (MQ=255)
                               tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                                      >  1:1910621/1‑69 (MQ=255)
                                aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                     >  1:1487992/1‑69 (MQ=255)
                                aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                                     >  1:206445/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:964196/68‑1 (MQ=255)
                                              tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:675052/68‑1 (MQ=255)
                                              tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:635200/68‑1 (MQ=255)
                                              tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:1895089/68‑1 (MQ=255)
                                              tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:3033260/68‑1 (MQ=255)
                                              tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:194329/68‑1 (MQ=255)
                                              tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                        <  1:2835426/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:505960/1‑67 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:12834/1‑67 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:187804/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:419064/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:3083514/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:2692847/1‑67 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:552721/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                >  1:2963043/1‑67 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2233189/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2883012/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2846702/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:1066511/67‑1 (MQ=255)
                                                       aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                <  1:2505300/67‑1 (MQ=255)
                                                           ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc          <  1:2593498/69‑1 (MQ=255)
                                                           ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc          <  1:2381755/69‑1 (MQ=255)
                                                              cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:763862/69‑1 (MQ=255)
                                                              cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:1239142/69‑1 (MQ=255)
                                                              cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt       <  1:2935050/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:2820023/1‑68 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc    >  1:2399167/1‑68 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc    >  1:2087856/1‑68 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:477606/1‑69 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:187557/1‑69 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:3341696/1‑69 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg   >  1:2841937/1‑69 (MQ=255)
                                                                  atTCTTGCCTCCGCGCTCGGCACTATTGGTGGGATGCGttt                               >  1:2824633/1‑41 (MQ=255)
                                                                   ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa  >  1:1280030/1‑69 (MQ=255)
                                                                   ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa  >  1:3246407/1‑69 (MQ=255)
                                                                    |                                                                   
GAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGA  >  minE/1781616‑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: