Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,707,052 C→A 10.8% R108L (CGC→CTC)  yphB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,707,0520CA10.8% 297.7 / 18.0 129R108L (CGC→CTC) yphBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (70/45);  new base A (8/6);  total (78/51)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 129 total aligned reads displayed.

CGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCA  >  minE/1706985‑1707112
                                                                   |                                                            
cGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCg                                                             <  1:166454/69‑1 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAACCGCCTGACTTACCCGATAGTGATAGCCACCGCTGCGATg                                                          >  1:2043402/1‑69 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGAt                                                           >  1:1524774/1‑68 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                          >  1:1018230/1‑69 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                          >  1:261765/1‑69 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                          >  1:2631197/1‑69 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                          >  1:2730696/1‑69 (MQ=255)
   cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                          >  1:1576185/1‑69 (MQ=255)
      cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                        >  1:1007881/1‑68 (MQ=255)
      cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                        >  1:3136175/1‑68 (MQ=255)
      cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                        >  1:576145/1‑68 (MQ=255)
      cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                        >  1:83775/1‑68 (MQ=255)
      cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                        >  1:910088/1‑68 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:1444378/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:2942727/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:1505206/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:733795/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:764147/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:2576270/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:250157/1‑69 (MQ=255)
       gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                      >  1:2288855/1‑69 (MQ=255)
                tCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:819556/1‑69 (MQ=255)
                tCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:2457254/1‑69 (MQ=255)
                 ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:91387/1‑68 (MQ=255)
                 ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:3365613/1‑68 (MQ=255)
                 ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:2451192/1‑68 (MQ=255)
                 ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:247306/1‑68 (MQ=255)
                 ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:529796/1‑68 (MQ=255)
                 ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                             >  1:1319449/1‑68 (MQ=255)
                  cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                            <  1:33810/68‑1 (MQ=255)
                  cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                            <  1:665284/68‑1 (MQ=255)
                  cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                            <  1:786330/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:1262358/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:613446/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:1300741/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:46744/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:3172951/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:3042404/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:1209724/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:1185779/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:2078093/68‑1 (MQ=255)
                                         ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                     <  1:1855154/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:2335742/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:1139669/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:464605/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:437108/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:641449/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:3002865/68‑1 (MQ=255)
                                             ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:1668845/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:1003143/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:700534/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:2958422/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:2644411/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:3078880/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:679054/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:2525517/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:1836660/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:1882088/68‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATGTGCATAGACCAGACACAAACTATCATCGCTGTGCGAGacac                 <  1:2847347/67‑1 (MQ=255)
                                              cGATAGTGATAGACACCGCTGCGATCTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:210285/68‑1 (MQ=255)
                                               gATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                <  1:1919366/67‑1 (MQ=255)
                                                 tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc            <  1:2666356/69‑1 (MQ=255)
                                                 tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc            <  1:2231454/69‑1 (MQ=255)
                                                 tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc            <  1:507764/69‑1 (MQ=255)
                                                 tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc            <  1:1636476/69‑1 (MQ=255)
                                                 tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc            <  1:1441623/69‑1 (MQ=255)
                                                        aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg      <  1:3345398/68‑1 (MQ=255)
                                                        aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg      <  1:2865798/68‑1 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:3349733/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2351151/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:1112239/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2414976/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:50941/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:1785264/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2500704/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:368392/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:364370/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:1707165/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2091840/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:251454/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2331494/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2220828/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2996514/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca  >  1:2869830/1‑67 (MQ=255)
                                                             ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTCCGAGACACATTGCCACTCGCcca  >  1:1765849/1‑67 (MQ=255)
                                                              cGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGGGCGAGAcaca                <  1:160136/52‑1 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:885852/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:1466124/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:344522/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:1328951/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:3259605/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2929578/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2702372/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2631757/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:1638175/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2530992/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2454554/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2145725/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:2013886/1‑41 (MQ=255)
                                                               gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                          >  1:1965604/1‑41 (MQ=255)
                                                                   |                                                            
CGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCA  >  minE/1706985‑1707112

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: