Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,370,824 C→T 36.9% intergenic (+81/‑49) yohK → / → cdd predicted inner membrane protein/cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,370,8240CT36.9% 89.9 / 82.8 135intergenic (+81/‑49)yohK/cddpredicted inner membrane protein/cytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base C (65/20);  new base T (22/28);  total (87/48)
Fisher's exact test for biased strand distribution p-value = 1.92e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCA  >  minE/1370765‑1370889
                                                           |                                                                 
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:732606/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:3281149/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:3198090/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:3368058/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:307231/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:534698/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:620941/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2884906/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2739555/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2695947/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2664306/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2389352/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2656256/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2469857/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:1229544/1‑68 (MQ=255)
ttttGATGTAAGTTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGtt                                                           >  1:2592319/1‑68 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGTGTTCGTAAACTGt                                                 <  1:978365/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:1043714/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:2136325/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:2130460/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:2040638/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:786912/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:3363772/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:3435146/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGtt                                                >  1:1450234/1‑69 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:3110690/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:3374554/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:2945197/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:499556/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:2646540/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:931527/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:1336069/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:283619/68‑1 (MQ=255)
          tgtTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 >  1:2449291/1‑68 (MQ=255)
          tgtTTAACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:2631052/68‑1 (MQ=255)
           gtTTCACGCGTTGCTTAATTAATGAGATTTAGATCTCATATAAAGCCATAACGGGTTCGt                                                        <  1:816283/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTcct                                                        <  1:2924418/60‑3 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:325369/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1175100/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:761436/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:2800249/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:724132/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:695698/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:3162733/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:320027/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:2555005/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1176714/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1586178/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1586402/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1792096/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1812654/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1902818/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1936244/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:1961323/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:2199141/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:227857/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:2528455/60‑1 (MQ=255)
           gtTTCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:2346236/60‑1 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTTAGATCACATATAAAGCCATAACGGGTTCGt                                                        <  1:873215/59‑1 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:2819855/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:1511335/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:1519563/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:381339/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:2396878/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:1818321/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:3206310/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:1948506/1‑69 (MQ=255)
            tttCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTAt                                              >  1:2053549/1‑69 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACaa                                                                 >  1:444228/1‑44 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAt                                         >  1:2312747/1‑68 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAt                                         >  1:1345906/1‑68 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAt                                         >  1:779426/1‑68 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAt                                         >  1:2633029/1‑68 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAt                                         >  1:426750/1‑68 (MQ=255)
                  gcgTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCACAt                                         >  1:844568/1‑68 (MQ=255)
                  gcgTTCCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAt                                         >  1:1583571/1‑68 (MQ=255)
                       gCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATg                                   >  1:2340138/1‑69 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTa                                                       >  1:2140573/1‑44 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:1649647/50‑1 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 >  1:462090/1‑50 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:475325/50‑1 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 >  1:1636589/1‑50 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 >  1:2665537/1‑50 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 >  1:2982313/1‑50 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGt                                                 <  1:348621/50‑1 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGAtt                                >  1:2853463/1‑67 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGAtt                                >  1:1394948/1‑67 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGAtt                                >  1:1430739/1‑67 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGAtt                                >  1:2743122/1‑67 (MQ=255)
                            attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGAtt                                >  1:2015552/1‑67 (MQ=255)
                            attaatGAGATTCAGATAACATATAAAGCCACAACGGGttg                                                          >  1:2230846/1‑40 (MQ=39)
                                 tGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGc                         >  1:60293/1‑69 (MQ=255)
                                  gagaTTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGc                         <  1:2567158/68‑1 (MQ=255)
                                     aTTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAAc                      >  1:1968277/1‑68 (MQ=255)
                                       tCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGc                    <  1:3280372/68‑1 (MQ=255)
                                                   taAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACTCCATGCATCCACGt       >  1:687340/1‑69 (MQ=255)
                                                   taAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGt       >  1:1387237/1‑69 (MQ=255)
                                                   taAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGt       >  1:855790/1‑69 (MQ=255)
                                                   taAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGt       >  1:1675033/1‑69 (MQ=255)
                                                        ccACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCa  <  1:992067/69‑1 (MQ=255)
                                                           |                                                                 
TTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCA  >  minE/1370765‑1370889

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: