Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,299,439 A→G 13.9% intergenic (‑364/+126) 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,4390AG13.9% 105.3 / 8.6 43intergenic (‑364/+126)cobU/yeeXbifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base G (2/4);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACA  >  minE/1299369‑1299504
                                                                      |                                                                 
aTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGaa                                                                   >  1:2986970/1‑71 (MQ=255)
 tCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAAt                                                                  <  1:1399187/71‑1 (MQ=255)
      gCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATcc                                                             <  1:2646719/71‑1 (MQ=255)
         aTTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCgg                                                           >  1:2517659/1‑70 (MQ=255)
          ttGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCg                                                            <  1:3277954/68‑1 (MQ=255)
           tGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAt                                                        >  1:1415213/1‑71 (MQ=255)
            ggttCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGAtg                                                       >  1:513489/3‑71 (MQ=255)
            gATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtg                                                       >  1:387079/1‑71 (MQ=255)
            gATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAt                                                        <  1:1958453/70‑1 (MQ=255)
              ttCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTAt                                                               <  1:2362783/61‑1 (MQ=255)
               tCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg                                                    >  1:1656376/1‑71 (MQ=255)
                cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTAt                                                               >  1:3089553/1‑59 (MQ=255)
                cAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg                                                    <  1:299404/70‑1 (MQ=255)
                 agcGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGc                                                  <  1:1515239/69‑1 (MQ=255)
                 aaCGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggt                                                     >  1:2573165/1‑68 (MQ=255)
                 aaCGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGc                                                  >  1:1999654/1‑71 (MQ=255)
                   cGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGc                                                  >  1:2213787/1‑69 (MQ=255)
                     tCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACt                                              <  1:176595/71‑1 (MQ=255)
                     tCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGATGGTGGCGagc                                              <  1:2394593/71‑3 (MQ=255)
                        gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatg                                           <  1:3785583/71‑2 (MQ=255)
                        gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatg                                           <  1:3822283/71‑2 (MQ=255)
                         cGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaa                                          >  1:1357433/1‑71 (MQ=255)
                          gACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACtatgca                                         >  1:1481928/1‑68 (MQ=255)
                              gATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGa                                                          <  1:3212400/50‑1 (MQ=255)
                               aTCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGa                                    >  1:1018209/1‑71 (MQ=255)
                                tCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATcc                                                             >  1:33409/1‑45 (MQ=255)
                                   tcgcTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCgg                                                           <  1:1036078/43‑1 (MQ=255)
                                   tcgcTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCgg                                                           <  1:1508678/43‑1 (MQ=255)
                                    cgcTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGTAATGGTGGCGACTATAATTTCGATTg                                 <  1:1071593/69‑1 (MQ=255)
                                    cgcTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGAt                                   <  1:3443648/67‑1 (MQ=255)
                                      cTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTtc                             >  1:321380/1‑71 (MQ=255)
                                          aaaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCg                          >  1:433460/1‑70 (MQ=255)
                                           aaaGAACTTGTAAAACAACGTAATAGAATTATCCGGAAtggtgg                                                   >  1:2604181/1‑44 (MQ=255)
                                              gAACTTGTAAAACAACGTAATAGAATTATCCGGATTGGTGGCGACTATAATTTCGATTGGTTCTCggat                       >  1:3381622/1‑69 (MQ=255)
                                               aaCTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTgg                                <  1:3322195/59‑1 (MQ=255)
                                               aaCTTGTAAAACAACGTAATAGAATTATCCGGAATGGAGGCGACTATAATTTCGATTGGTTCTCggatgag                    >  1:1299343/1‑71 (MQ=255)
                                                 cTTGTAAAACAACGTAATAGAATTATCCGGAAtggtg                                                    >  1:3760524/1‑37 (MQ=255)
                                                 cTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCggatgagga                  <  1:3250968/71‑1 (MQ=255)
                                                     tAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGgatgat              <  1:1762734/71‑1 (MQ=255)
                                                              gTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACgg         <  1:2348389/67‑1 (MQ=255)
                                                                 aTAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca  >  1:92975/1‑71 (MQ=255)
                                                                 aTAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATaca  >  1:3738676/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)
                                                                      |                                                                 
ATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACA  >  minE/1299369‑1299504

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: