Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 785,203 C→T 41.2% A14A (GCG→GCA dsbB ← oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE785,2030CT41.2% 40.4 / 69.0 97A14A (GCG→GCAdsbBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base T (25/15);  total (57/40)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

TCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAG  >  minE/785139‑785271
                                                                |                                                                    
tCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACgcgc                                                                  >  1:2932938/1‑69 (MQ=255)
tCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACgcgc                                                                  >  1:2069693/1‑69 (MQ=255)
tCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACgcgc                                                                  >  1:187393/1‑69 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc                                                           >  1:1014482/1‑69 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc                                                           >  1:2602714/1‑69 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc                                                           >  1:2049371/1‑69 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc                                                           >  1:1506028/1‑69 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc                                                           >  1:663925/1‑69 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGc                                                            >  1:354652/1‑68 (MQ=255)
       cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGc                                                            >  1:1823071/1‑68 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:2985039/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:1317850/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:2545692/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:3203644/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:1784179/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:206865/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:2319687/69‑1 (MQ=255)
              cacaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa                                                    <  1:2221457/69‑1 (MQ=255)
                                   gCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCaa                                >  1:480189/1‑68 (MQ=255)
                                   gCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCaa                                >  1:664688/1‑68 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCTCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:2500989/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:3260902/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:2195888/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:2171464/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:1124187/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:3240237/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:1129500/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:870412/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:2010537/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:1934895/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:2674292/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:2815384/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:584688/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:1589182/69‑1 (MQ=255)
                                            gTAAACGCCATCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCAtt                      <  1:1858541/69‑1 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATaa                           >  1:1170382/1‑58 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTcc                   >  1:1075701/1‑66 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:1423677/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:1942145/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:2825615/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:573662/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:2669957/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:2125068/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:3310503/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:157899/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:60683/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:3289234/1‑68 (MQ=255)
                                                  gCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:1473640/1‑68 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCt                 >  1:174591/1‑67 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               <  1:499034/69‑1 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               <  1:2618031/69‑1 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               >  1:127142/1‑69 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               >  1:2733088/1‑69 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               >  1:563285/1‑69 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               >  1:908930/1‑69 (MQ=255)
                                                   ccATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGc               >  1:353430/1‑69 (MQ=255)
                                                   ccATCAACAGCCACGCGCCACGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTg                >  1:1412521/1‑68 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAAc                               >  1:3424036/1‑50 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:335133/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1037980/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:920079/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:701229/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:3141438/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:939650/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:3131648/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:3108847/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:798496/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1892861/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1048998/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1166514/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1310607/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1376217/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1580040/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:167295/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1741168/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:1866098/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:2991638/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:2947545/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:2780308/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:2144768/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCCCAGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:2983786/1‑51 (MQ=255)
                                                      tCAACAGCCATGCGCACCGGCCTTGTGAACATTGGTTCAAAAATCGCAACa                              >  1:655672/1‑51 (MQ=255)
                                                        aaCAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACAt                             <  1:1168374/50‑1 (MQ=255)
                                                        aaCAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACAt                             <  1:2275407/50‑1 (MQ=255)
                                                          cAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATg         <  1:2042343/68‑1 (MQ=255)
                                                           aGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCa       <  1:1536804/69‑1 (MQ=255)
                                                           aGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCa       <  1:2202012/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:2364639/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:1916523/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:133687/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:2878626/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:1272901/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:192061/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:2018692/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:925110/69‑1 (MQ=255)
                                                              cACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:341445/69‑1 (MQ=255)
                                                               aCGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATag    <  1:328724/68‑1 (MQ=255)
                                                                cgcgCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATagag  <  1:1427709/69‑1 (MQ=255)
                                                                |                                                                    
TCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAG  >  minE/785139‑785271

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: