Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,258,429 A→G 9.5% intergenic (+100/+65) yecP → / ← torZ predicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,258,4290AG 9.5% 163.7 / 7.6 63intergenic (+100/+65)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (3/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAC  >  minE/1258360‑1258497
                                                                     |                                                                    
aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCa                                                                      <  1:1560218/70‑1 (MQ=255)
  cGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAgc                                                                  <  1:3361755/71‑1 (MQ=9)
  cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAg                                                                   >  1:3384363/1‑71 (MQ=255)
    ccTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                                                 >  1:2996392/1‑71 (MQ=255)
        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                                              <  1:2422172/69‑1 (MQ=9)
        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                                              >  1:2934480/1‑69 (MQ=9)
        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                                              >  1:3007148/1‑69 (MQ=9)
        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                                              <  1:1781149/69‑1 (MQ=9)
         tCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                                              >  1:714290/1‑69 (MQ=255)
            ggCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGAGTAGCGCATCAGGct                                                        >  1:3661923/1‑70 (MQ=9)
             gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAgg                                                          >  1:3797587/1‑69 (MQ=255)
              ccTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAg                                                       >  1:2191031/1‑71 (MQ=255)
                tACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                                                 <  1:1502178/59‑1 (MQ=255)
                tACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                                                 <  1:2401853/59‑1 (MQ=255)
                  cAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                                            >  1:2874553/1‑62 (MQ=255)
                   aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTa                                                  >  1:3594583/1‑71 (MQ=255)
                       ccTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGtt                                              >  1:2534001/1‑71 (MQ=255)
                           gCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                                             >  1:3749742/1‑68 (MQ=255)
                           gCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGc                                           <  1:2375998/70‑1 (MQ=255)
                            cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcgc                                                                >  1:1255422/1‑48 (MQ=255)
                            cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                                         <  1:3446304/71‑1 (MQ=255)
                               aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                                        <  1:54899/69‑1 (MQ=255)
                               aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                                        <  1:2988513/69‑1 (MQ=255)
                               aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                                      <  1:3152461/71‑1 (MQ=255)
                               aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                                      >  1:791146/1‑71 (MQ=255)
                                ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCAt                                                                     <  1:287292/39‑1 (MQ=9)
                                ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                                        >  1:594385/1‑68 (MQ=255)
                                ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAcc                                     <  1:1051038/71‑1 (MQ=255)
                                ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAcc                                     >  1:3091160/1‑71 (MQ=255)
                                 tCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                                                    >  1:3411357/1‑55 (MQ=255)
                                  cAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                                               >  1:1087146/1‑59 (MQ=9)
                                  cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                                    <  1:3874387/70‑1 (MQ=255)
                                  cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                                    <  1:3869756/70‑1 (MQ=255)
                                  cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                                    >  1:1664035/1‑70 (MQ=255)
                                    atatatTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                                                 <  1:2482192/39‑1 (MQ=255)
                                     tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCttt                                                   >  1:199104/1‑52 (MQ=255)
                                     tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCt                                                     >  1:1520544/1‑50 (MQ=255)
                                         ttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc                            >  1:874143/1‑71 (MQ=255)
                                         ttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc                            <  1:453880/71‑1 (MQ=255)
                                          tGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc                            >  1:2095191/1‑70 (MQ=255)
                                            cAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGt                          <  1:944924/70‑1 (MQ=255)
                                            cAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGt                          <  1:3854815/70‑1 (MQ=255)
                                             aGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGttt                        <  1:3293059/71‑1 (MQ=255)
                                               ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACg                                                <  1:3300766/45‑1 (MQ=255)
                                               ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTcc                      <  1:1931014/71‑1 (MQ=255)
                                               ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTcc                      <  1:3560809/71‑1 (MQ=255)
                                                    aTGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                                      >  1:2049194/1‑50 (MQ=255)
                                                    aTGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                                      >  1:1884088/1‑50 (MQ=255)
                                                     tGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgcgc                              >  1:791239/1‑57 (MQ=255)
                                                     tGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgcgc                              >  1:1514191/1‑57 (MQ=255)
                                                      gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGTTTgt                                           <  1:3617267/43‑2 (MQ=11)
                                                      gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGTTTgt                                           <  1:2320350/43‑2 (MQ=11)
                                                      gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGTTTgt                                           <  1:409011/43‑2 (MQ=11)
                                                         ggCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTAt             <  1:1309310/70‑1 (MQ=255)
                                                          accTGATAAGCATAGCGCATCAGGCTGCTTTGCGtt                                              <  1:21456/35‑1 (MQ=14)
                                                            cTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgc                                >  1:3824594/1‑48 (MQ=255)
                                                            cTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGttt                        <  1:2460936/56‑1 (MQ=255)
                                                            cTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCAc         >  1:2537592/1‑71 (MQ=255)
                                                             tGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCAc         >  1:1352657/1‑70 (MQ=255)
                                                              gATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTcc                      >  1:1925341/1‑56 (MQ=255)
                                                               aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTTAGCGCCCGTTTCCGGGCGCTATTCACGTc      >  1:2022117/1‑71 (MQ=255)
                                                                  aGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCtt    <  1:3014099/70‑1 (MQ=255)
                                                                   gCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGgcg                 <  1:3347696/56‑1 (MQ=255)
                                                                    cATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAc  <  1:52087/70‑1 (MQ=255)
                                                                     aTAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAg                                 >  1:1115070/1‑38 (MQ=255)
                                                                     |                                                                    
AACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTAC  >  minE/1258360‑1258497

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