Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,299,487 G→A 30.8% intergenic (‑412/+78) cobU ← / ← yeeX bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,299,4870GA30.8% 48.2 / 28.3 39intergenic (‑412/+78)cobU/yeeXbifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/7);  new base A (5/7);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 7.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCA  >  minE/1299418‑1299548
                                                                     |                                                             
cTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgagga                                                              <  1:3250968/71‑1 (MQ=255)
    tAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGgatgat                                                          <  1:1762734/71‑1 (MQ=255)
             gTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACgg                                                     <  1:2348389/67‑1 (MQ=255)
                aTAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca                                              >  1:3738676/1‑71 (MQ=255)
                aTAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca                                              >  1:92975/1‑71 (MQ=255)
                 tAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca                                              >  1:3796255/1‑70 (MQ=255)
                  aGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca                                              >  1:2246264/1‑69 (MQ=255)
                    aaTTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggatg                                                            <  1:1116850/53‑1 (MQ=255)
                        aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggatg                                                            <  1:3350800/49‑1 (MQ=255)
                        aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgaggatg                                                            <  1:2820589/49‑1 (MQ=255)
                        aTCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGAt                                         >  1:3014685/1‑68 (MQ=255)
                             gAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATca                                 <  1:2027807/71‑1 (MQ=255)
                               aTGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                               >  1:604844/1‑71 (MQ=255)
                               aTGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                               >  1:600241/1‑71 (MQ=255)
                                tggtggCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                               >  1:910351/1‑70 (MQ=255)
                                 ggtggCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                               >  1:2527057/1‑69 (MQ=255)
                                  gtggCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgag                            >  1:1744817/1‑71 (MQ=255)
                                     gCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagca                          >  1:1586721/1‑70 (MQ=255)
                                     gCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgaaca                          >  1:1886354/1‑70 (MQ=255)
                                       gACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACg                                                      >  1:3731390/1‑40 (MQ=255)
                                       gACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGAt                                                 >  1:3283717/1‑45 (MQ=255)
                                        aCTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGt                                                    <  1:593159/41‑1 (MQ=255)
                                         cTATAATTTCGATTGGTTCTCGGATGAGGATTATGACGGTGATACACAGATAAGAATcaatgagcaatgag                     <  1:3739297/71‑1 (MQ=255)
                                           ttaATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAgg                    <  1:2018437/69‑1 (MQ=255)
                                            taATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                               >  1:2751144/1‑58 (MQ=255)
                                            taATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaat                               >  1:2055319/1‑58 (MQ=255)
                                             aaTTTCTATTGGTTCTCGGATGAGGATGATGACGGt                                                    >  1:504936/1‑36 (MQ=255)
                                             aaTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGtt                 >  1:2856136/1‑71 (MQ=255)
                                             aaTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGtt                 >  1:3332354/1‑71 (MQ=255)
                                             aaTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGt                  <  1:3812496/70‑1 (MQ=255)
                                              aTTTCGATTGGTTCTCGGATGAGAATGATGACTGTGATACACAGATAAGAATCAATGAGCAATGAGGGtt                 <  1:1358751/70‑1 (MQ=255)
                                              aTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATcaatgagcaatga                      <  1:3431210/65‑1 (MQ=255)
                                               tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGc               <  1:2733198/71‑1 (MQ=255)
                                                 tCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc              <  1:3546122/70‑1 (MQ=255)
                                                 tCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc              <  1:1949429/70‑1 (MQ=255)
                                                 tCGATTGGTTCTCGGATGAGAATGATGACCGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc              <  1:122121/70‑1 (MQ=255)
                                                  cGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTg                >  1:98555/1‑67 (MQ=255)
                                                  cGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg            >  1:799681/1‑71 (MQ=255)
                                                  cGATTGGTTCTCGGATGAGAATGATGACGGTGATa                                                >  1:1583546/1‑35 (MQ=255)
                                                   gATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagcaatga                      <  1:3513277/60‑1 (MQ=255)
                                                   gATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagcaatga                      <  1:3089664/60‑1 (MQ=255)
                                                   gATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCggg           >  1:1976288/1‑71 (MQ=255)
                                                      tGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagc                           >  1:1376066/1‑52 (MQ=255)
                                                          tctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCt    >  1:1686008/1‑71 (MQ=255)
                                                              ggatgaggatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCa  <  1:826620/69‑1 (MQ=255)
                                                              ggatgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTc   <  1:1155072/68‑1 (MQ=255)
                                                               gatgagaatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCa  >  1:2661686/1‑68 (MQ=255)
                                                                  gaggatgaTGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGtt                 >  1:3904441/1‑50 (MQ=255)
                                                                     gatgatGACGGTGATACACAGATAAGAATcaatgagcaatgag                     <  1:3556131/43‑1 (MQ=255)
                                                                     |                                                             
CTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCA  >  minE/1299418‑1299548

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: