Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,707,042 G→A 15.0% V111V (GTC→GTT yphB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,707,0420GA15.0% 161.2 / 13.7 92V111V (GTC→GTTyphBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/47);  new base A (8/6);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCG  >  minE/1706976‑1707108
                                                                  |                                                                  
gACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGac                                                                  <  1:587997/69‑1 (MQ=255)
gACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGac                                                                  <  1:3226802/69‑1 (MQ=255)
gACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGac                                                                  <  1:1435136/69‑1 (MQ=255)
      gagCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCt                                                            >  1:1244421/1‑69 (MQ=255)
         cGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCg                                                         <  1:166454/69‑1 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAACCGCCTGACTTACCCGATAGTGATAGCCACCGCTGCGATg                                                      >  1:2043402/1‑69 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGAt                                                       >  1:1524774/1‑68 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGAt                                                       >  1:2062902/1‑68 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                      >  1:2730696/1‑69 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                      >  1:1576185/1‑69 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                      >  1:261765/1‑69 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                      >  1:2631197/1‑69 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                      >  1:780706/1‑69 (MQ=255)
            cACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATg                                                      >  1:1018230/1‑69 (MQ=255)
               cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                    >  1:1007881/1‑68 (MQ=255)
               cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                    >  1:3136175/1‑68 (MQ=255)
               cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                    >  1:245638/1‑68 (MQ=255)
               cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                    >  1:2372764/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)
               cGTCAGCGTATCCGCCGTTAAATGAAATGCCTGGCTTACCCGATAGTGATAAACACCGCTGAGATGtt                                                    >  1:1469525/1‑68 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:1505206/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:1444378/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:2288855/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:733795/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:764147/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:250157/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:2576270/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:2636827/1‑69 (MQ=255)
                gTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCa                                                  >  1:2942727/1‑69 (MQ=255)
                        aTCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGac                                                                  <  1:2952141/45‑1 (MQ=255)
                         tCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:819556/1‑69 (MQ=255)
                         tCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:1380583/1‑69 (MQ=255)
                         tCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:2457254/1‑69 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:3374497/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:1319449/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:247306/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:529796/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:2451192/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:3365613/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:1905363/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:1800038/1‑68 (MQ=255)
                          ccgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGa                                         >  1:91387/1‑68 (MQ=255)
                           cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                        <  1:665284/68‑1 (MQ=255)
                           cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                        <  1:3051966/68‑1 (MQ=255)
                           cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                        <  1:33810/68‑1 (MQ=255)
                           cgccgTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGac                                        <  1:786330/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:3042404/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:1300741/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:3172951/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:1262358/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:46744/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:1209724/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:613446/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:2078093/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:1185779/68‑1 (MQ=255)
                                                  ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                 <  1:1855154/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:464605/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:437108/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:588679/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:641449/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:1139669/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:1668845/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:3293655/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:3002865/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:2335742/68‑1 (MQ=255)
                                                      ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac             <  1:2504919/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:1882088/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:700534/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:679054/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:674856/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:2525517/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:2538249/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:1836660/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:2644411/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:1003143/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:2958422/68‑1 (MQ=255)
                                                       cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca            <  1:3078880/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:1441623/69‑1 (MQ=255)
                                                          tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc        <  1:507764/69‑1 (MQ=255)
                                                          tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc        <  1:2666356/69‑1 (MQ=255)
                                                          tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc        <  1:401521/69‑1 (MQ=255)
                                                          tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc        <  1:2231454/69‑1 (MQ=255)
                                                          tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc        <  1:1636476/69‑1 (MQ=255)
                                                                 aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg  <  1:1664555/68‑1 (MQ=255)
                                                                 aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg  <  1:1244677/68‑1 (MQ=255)
                                                                 aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg  <  1:2865798/68‑1 (MQ=255)
                                                                 aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg  <  1:3345398/68‑1 (MQ=255)
                                                                  |                                                                  
GACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCG  >  minE/1706976‑1707108

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: