Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,707,088 T→C 24.9% H96R (CAC→CGC)  yphB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,707,0880TC24.9% 148.2 / 46.0 125H96R (CAC→CGC) yphBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/59);  new base C (24/7);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 1.30e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 140 total aligned reads displayed.

TTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATC  >  minE/1707026‑1707143
                                                              |                                                       
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:46744/68‑1 (MQ=255)
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:3172951/68‑1 (MQ=255)
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:1185779/68‑1 (MQ=255)
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:1209724/68‑1 (MQ=255)
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:2078093/68‑1 (MQ=255)
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:1262358/68‑1 (MQ=255)
ttACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCgag                                                    <  1:1855154/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:641449/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:464605/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:588679/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:2504919/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:3002865/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:2335742/68‑1 (MQ=255)
    ccGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGacac                                                <  1:3293655/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:1836660/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:2538249/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:2644411/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:1003143/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:3078880/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:679054/68‑1 (MQ=255)
     cGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGAcaca                                               <  1:700534/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:507764/69‑1 (MQ=255)
        tAGTGATAAACACCGCTGAGATGTTCATACACCAGACACAAGCTATCATCGCTGCGCGAGACACATTGc                                           <  1:401521/69‑1 (MQ=255)
               aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCg                                     <  1:1664555/68‑1 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2869830/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2996514/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:3349733/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2500704/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:251454/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2414976/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2351151/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:364370/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2331494/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:272944/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:50941/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2220828/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:871035/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:2091840/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:1785264/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:966177/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:1112239/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:1119431/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:1707165/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                 >  1:1626135/1‑67 (MQ=255)
                    ccGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTCCGAGACACATTGCCACTCGCcca                                 >  1:1765849/1‑67 (MQ=255)
                     cGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGGGCGAGAcaca                                               <  1:160136/52‑1 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2631757/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:1328951/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2929578/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2883195/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:1466124/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2702372/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:1638175/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2530992/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:1965604/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2454554/1‑41 (MQ=255)
                      gCTGCGATGTTCATACACCAGACACAAACTATCATCGCtgt                                                         >  1:2145725/1‑41 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGcc                    >  1:1240862/1‑64 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGcc                    >  1:2128303/1‑64 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:3199790/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:479002/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:3405830/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:3408903/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:827053/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:650320/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:431721/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:589254/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:257763/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:2975125/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:1729313/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:1842155/1‑65 (MQ=255)
                                    acacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:2776716/1‑65 (MQ=255)
                                    aaacCAGACACAAGCTATCATCGCTGCGCGAGACACATTGCCACTCGCCCAACCAGCCATCGCCg                   >  1:358920/3‑65 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCATCCAGCCATCGCCGTGTAGATAg          <  1:864330/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:3153017/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:972974/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:1512909/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:2223151/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:2272111/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:2275691/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:262821/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:401098/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:2966859/69‑1 (MQ=255)
                                         aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg          <  1:3173154/69‑1 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCGTCGCCGTGTAGATAGTGCGCAt   >  1:2809766/1‑69 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCAt   >  1:1001229/1‑69 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCAt   >  1:2223367/1‑69 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCAt   >  1:615729/1‑69 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCAt   >  1:1702347/1‑69 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCAt   >  1:3323650/1‑69 (MQ=255)
                                                aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCAt   >  1:1009052/1‑69 (MQ=255)
                                                 aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATc  <  1:16213/69‑1 (MQ=255)
                                                 aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATc  <  1:3008822/69‑1 (MQ=255)
                                                 aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATc  <  1:1430342/69‑1 (MQ=255)
                                                 aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATc  <  1:1426023/69‑1 (MQ=255)
                                                 aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATc  <  1:1290987/69‑1 (MQ=255)
                                                              |                                                       
TTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATC  >  minE/1707026‑1707143

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: