Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,434,773 T→C 48.7% A465A (GCA→GCG yojI ← fused predicted multidrug transport subunits and membrane component and ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,434,7730TC48.7% ‑2.3 / 54.3 39A465A (GCA→GCGyojIfused predicted multidrug transport subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (11/9);  new base C (9/10);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACTCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCG  >  minE/1434703‑1434842
                                                                      |                                                                     
aaaaCTCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTtcttct                                                                       >  1:2341957/1‑71 (MQ=255)
aaaaCTCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTtcttct                                                                       >  1:1855574/1‑71 (MQ=255)
   aCTCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTgcc                                                                    <  1:3063575/71‑1 (MQ=255)
       acgacgGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCCgccagcg                                                                >  1:592247/1‑71 (MQ=255)
        cgacgGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTgccagcgc                                                               <  1:1473491/71‑1 (MQ=255)
          acgGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTgccagcgc                                                               >  1:626413/1‑69 (MQ=255)
          acgGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCCgccagcgcca                                                             <  1:3278529/71‑1 (MQ=255)
           cgGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCcgccagcgccag                                                            >  1:1926044/1‑71 (MQ=255)
                gtgtGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCaa                                                         >  1:2388169/1‑69 (MQ=255)
                gtgtGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACa                                                       <  1:78968/71‑1 (MQ=255)
                   tgGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAgcg                                                    >  1:3684295/1‑71 (MQ=255)
                   tgGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAgc                                                     <  1:1201628/70‑1 (MQ=255)
                        ccTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCAc                                                >  1:1883338/1‑70 (MQ=255)
                             tCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTgccagcgccagc                                                           <  1:3350898/54‑1 (MQ=255)
                             tCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGt                                            >  1:36391/1‑69 (MQ=255)
                             tCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGt                                            >  1:1396900/1‑69 (MQ=255)
                             tCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGttt                                          <  1:41755/71‑1 (MQ=255)
                                gccgccCATTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGtttttt                                       <  1:3125488/71‑1 (MQ=255)
                                 ccgccCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTc                                      <  1:3709886/71‑1 (MQ=255)
                                  cgccCATTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCt                                     >  1:1429297/1‑71 (MQ=255)
                                       aTTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCt                                >  1:709712/1‑71 (MQ=255)
                                       aTTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCt                                <  1:556907/71‑1 (MQ=255)
                                       aTTCATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCt                                <  1:1458767/71‑1 (MQ=255)
                                        ttCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCtt                               <  1:1759443/71‑1 (MQ=255)
                                          cATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCtttt                             >  1:3021820/1‑71 (MQ=255)
                                          cATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCtttt                             >  1:3518989/1‑71 (MQ=255)
                                          cATCCAGCAGGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCtttt                             >  1:1307678/1‑71 (MQ=255)
                                               agcagGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTg                            >  1:781621/1‑67 (MQ=255)
                                               agcagGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATaa                        >  1:1256479/1‑71 (MQ=255)
                                                  agGATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATaa                        >  1:3366865/1‑68 (MQ=255)
                                                   ggATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCttt                              <  1:2712021/61‑1 (MQ=255)
                                                    gATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACtt                     <  1:1426877/69‑1 (MQ=255)
                                                    gATAATATCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCa                   <  1:2882187/71‑1 (MQ=255)
                                                      taataTCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAgg                 >  1:2932208/1‑71 (MQ=255)
                                                      taataTCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAgg                 >  1:2179624/1‑71 (MQ=255)
                                                      taataTCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCacg                 >  1:124341/1‑69 (MQ=255)
                                                       aataTCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCtttt                             <  1:1750003/58‑1 (MQ=255)
                                                        atatCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGtt               <  1:3013724/71‑1 (MQ=255)
                                                         tatCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTa              <  1:276406/71‑1 (MQ=255)
                                                         tatCGCGTTCTTCCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTa              <  1:1306978/71‑1 (MQ=255)
                                                             gcgTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACaa          >  1:867669/1‑71 (MQ=255)
                                                                ttcttcCGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCt                                     >  1:2979089/1‑41 (MQ=255)
                                                                 tcttccGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACg      >  1:201973/1‑71 (MQ=255)
                                                                  cttccGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGc     >  1:929458/1‑71 (MQ=255)
                                                                   ttctGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTa              <  1:406964/61‑1 (MQ=255)
                                                                   ttccGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAAc            <  1:2239760/63‑1 (MQ=255)
                                                                    tctGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGc     <  1:1961052/69‑1 (MQ=255)
                                                                     ctGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCg  >  1:1637396/1‑71 (MQ=255)
                                                                      |                                                                     
AAAACTCACGACGGAAGTGTGGATCCTGATCCGCCGCCCATTCATCCAGCAGGATAATATCGCGTTCTTCTGCCAGCGCCAGCAACAGCGCCACGCGTTTTTTCTGCCCTTTTGATAACTTCAGGTTAACAATACGCCCG  >  minE/1434703‑1434842

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: