Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,716,832 A→C 18.3% intergenic (+90/+108) yphH → / ← glyA predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,716,8320AC18.3% 131.2 / 43.8 98intergenic (+90/+108)yphH/glyApredicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (49/31);  new base C (8/10);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAT  >  minE/1716763‑1716902
                                                                     |                                                                      
gCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTacac                                                                       >  1:79268/1‑71 (MQ=255)
    cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACg                                                                      >  1:725093/1‑68 (MQ=255)
    cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACg                                                                      >  1:5406542/1‑68 (MQ=255)
    cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:1278834/70‑1 (MQ=255)
       gaaggaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    >  1:4538483/1‑67 (MQ=255)
           gaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGc                                                              <  1:830322/69‑1 (MQ=255)
           gaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCaa                                                            <  1:1382750/71‑1 (MQ=255)
           gaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCaa                                                            <  1:4219521/71‑1 (MQ=255)
                  tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATAtt                                                      >  1:2158306/1‑70 (MQ=255)
                  tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATAtt                                                      >  1:5167321/1‑70 (MQ=255)
                       atCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGc                                                              <  1:4501282/57‑1 (MQ=255)
                       atCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGc                                                              <  1:1722259/57‑1 (MQ=255)
                        tCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGa                                                    >  1:5353944/1‑66 (MQ=255)
                        tCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                               >  1:2183963/1‑71 (MQ=255)
                         cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAttt                                                >  1:2953982/1‑69 (MQ=255)
                         cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAttt                                                >  1:2209273/1‑69 (MQ=255)
                         cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAttt                                                >  1:1408933/1‑69 (MQ=255)
                         cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                               >  1:4178114/1‑70 (MQ=255)
                           tAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAatatat                                                       <  1:1149807/60‑1 (MQ=255)
                            aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:5673330/1‑48 (MQ=255)
                             aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACgg                                                                     >  1:4679555/1‑44 (MQ=255)
                             aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTg                                                     >  1:2996130/1‑60 (MQ=255)
                             aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTg                                                     >  1:4525507/1‑60 (MQ=255)
                             aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAg                                          >  1:621064/1‑71 (MQ=255)
                                  caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                     <  1:1355129/70‑1 (MQ=11)
                                  caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                     <  1:1070969/70‑1 (MQ=11)
                                  caaaTTGCCTCATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                     <  1:4129553/70‑1 (MQ=11)
                                  aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAatatat                                                       >  1:4963288/1‑53 (MQ=255)
                                   agcttgcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTg                                               >  1:4318637‑M2/20‑60 (MQ=255)
                                   agcttgcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                     <  1:345139‑M2/51‑1 (MQ=255)
                                   agcttgcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                     <  1:334856‑M2/51‑1 (MQ=255)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:378560/1‑41 (MQ=40)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:3653164/1‑41 (MQ=40)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:4091429/1‑41 (MQ=40)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:3717694/1‑41 (MQ=40)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:1479918/1‑41 (MQ=40)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                  >  1:1042618/1‑41 (MQ=40)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                     >  1:1029988/1‑70 (MQ=11)
                                   aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGt                                    <  1:1405716/71‑1 (MQ=11)
                                    gcttgcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                     >  1:100149‑M2/19‑69 (MQ=255)
                                    aaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                    <  1:1965428/70‑1 (MQ=255)
                                    aaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGt                                    <  1:1123517/70‑1 (MQ=11)
                                     aTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAg                                  <  1:280107/71‑1 (MQ=255)
                                      ttgcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAata                                                          <  1:829910‑M2/30‑1 (MQ=255)
                                      ttGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAg                                  >  1:4964856/1‑70 (MQ=11)
                                        gcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAatata                                                        >  1:24882‑M2/15‑46 (MQ=255)
                                        gcctgatgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAatata                                                        >  1:3719605‑M2/15‑46 (MQ=255)
                                         ccTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAggac                               <  1:3626879/70‑3 (MQ=11)
                                         ccTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAgga                                <  1:5333843/69‑2 (MQ=11)
                                           tGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                            >  1:87508/1‑71 (MQ=255)
                                            gATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAt                                                  <  1:3524256/48‑1 (MQ=255)
                                             atgcgcgatGCTTATCAGGCCTACCAGAAGATTGCa                                                             >  1:4517591‑M2/10‑36 (MQ=255)
                                             atgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATa                          <  1:5482606‑M2/62‑1 (MQ=255)
                                             atgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATa                          <  1:276209‑M2/62‑1 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                    >  1:5241154/1‑61 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                    >  1:3472443/1‑61 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                    >  1:3465622/1‑61 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                            <  1:64705/69‑1 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                           <  1:1271959/70‑1 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                           <  1:5506823/70‑1 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGAt                           <  1:1233144/70‑1 (MQ=11)
                                             aTGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGAt                           <  1:4987178/70‑1 (MQ=11)
                                             aTGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATa                          <  1:2054670/71‑1 (MQ=11)
                                              tgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                     <  1:303370‑M2/51‑1 (MQ=255)
                                              tgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                     <  1:5484282‑M2/51‑1 (MQ=255)
                                                cgcgatGCTTATCAGGCCTACCAGAAGATTGCAata                                                          >  1:4175209‑M2/7‑36 (MQ=255)
                                                cgcTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                       >  1:1208984/1‑71 (MQ=255)
                                                cgcTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAg                        <  1:3703517/70‑1 (MQ=11)
                                                cgcTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAg                        <  1:4010410/70‑1 (MQ=11)
                                                  cgatGCTTATCAGGCCTACAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                    <  1:2714889‑M2/67‑1 (MQ=255)
                                                    atGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGtttt                                      <  1:5031911‑M2/50‑1 (MQ=255)
                                                    atGCTTATCAGGCCTACCAGAAGATAGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt                   <  1:4203362‑M2/69‑1 (MQ=255)
                                                    aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                            >  1:1349458/1‑62 (MQ=255)
                                                    aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                            >  1:4984271/1‑62 (MQ=255)
                                                    aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAgga                                >  1:4816340/1‑57 (MQ=2)
                                                    aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAgga                                >  1:5616760/1‑57 (MQ=2)
                                                    aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAgga                                >  1:5119812/1‑57 (MQ=2)
                                                    aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAgga                                >  1:293058/1‑57 (MQ=2)
                                                    aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgct                   >  1:4413369/1‑69 (MQ=2)
                                                     cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCaaa                                           >  1:2113532/1‑46 (MQ=255)
                                                     cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCaaa                                           >  1:1220582/1‑46 (MQ=255)
                                                     cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                   <  1:5436240/70‑1 (MQ=255)
                                                       cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                 >  1:2513460/1‑70 (MQ=255)
                                                       cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                 >  1:4071410/1‑70 (MQ=255)
                                                       cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc                >  1:3248722/1‑71 (MQ=255)
                                                       cTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATATGCGATGTTTTGTAGGACTGATAAGGCGCTcg                 >  1:5395870/1‑69 (MQ=2)
                                                         tATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAgga                                >  1:109984/1‑52 (MQ=2)
                                                           tCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTa                                   <  1:4110570/48‑1 (MQ=11)
                                                             aGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAgg                                 >  1:3164094/1‑48 (MQ=11)
                                                             aGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAgg                                 >  1:3956611/1‑48 (MQ=11)
                                                                ccTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                           <  1:5218794/51‑1 (MQ=255)
                                                                  tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                          <  1:1792834/50‑1 (MQ=255)
                                                                  tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                          <  1:529082/50‑1 (MQ=255)
                                                                    cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                    >  1:5162570/1‑38 (MQ=255)
                                                                    cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                         >  1:1550923/1‑49 (MQ=255)
                                                                    cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                        >  1:1556834/1‑50 (MQ=255)
                                                                    cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                        >  1:4125468/1‑50 (MQ=255)
                                                                    cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg     <  1:5050166/69‑1 (MQ=255)
                                                                     acGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                         >  1:5061555/1‑48 (MQ=255)
                                                                     acGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt  >  1:5566762/1‑71 (MQ=255)
                                                                     |                                                                      
GCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAT  >  minE/1716763‑1716902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: