Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,839 A→C 42.0% intergenic (+221/‑154) gyrB → / → ibpA DNA gyrase, subunit B/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,493,8390AC42.0% ‑3.1 / 39.2 31intergenic (+221/‑154)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base C (8/5);  total (18/13)
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

CGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAT  >  minE/2493771‑2493905
                                                                    |                                                                  
cGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaa                                                                   >  1:2567523/1‑70 (MQ=255)
       tttATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaaa                                                              <  1:2347191/68‑1 (MQ=255)
                 cAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTg                                                    <  1:1244286/68‑1 (MQ=255)
                   aTCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGaaaaaata                                                             >  1:898010/1‑55 (MQ=255)
                     ccGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTcc                                                        <  1:3851612/60‑1 (MQ=255)
                      cGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCt                                            >  1:1795689/1‑71 (MQ=255)
                        cAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAAGTGAAAAAAAAAGTCCGCTg                                                    <  1:3527744/61‑1 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCGGCTGATAAGGCTTGAAAAg                                     >  1:1188538/1‑71 (MQ=255)
                             atcatcCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                     >  1:138067/1‑71 (MQ=255)
                             atcatcCCCCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAg                                     >  1:1885073/1‑71 (MQ=255)
                                atcCCGCAAATCTACTATGCAATAGACCATAAACTGa                                                                    <  1:2395583/37‑1 (MQ=255)
                                   ccGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAt                                >  1:2458817/1‑70 (MQ=255)
                                              ctaTGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGAc                        <  1:1224153/67‑1 (MQ=255)
                                                 tGCAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaaa                                       >  1:2284503/1‑49 (MQ=255)
                                                  gCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCAt                                <  1:2447549/55‑1 (MQ=255)
                                                   cAATAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGaa                                        >  1:3222247/1‑46 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATaa                                                >  1:428101/1‑38 (MQ=255)
                                                   cAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATa                                                 >  1:4577/1‑37 (MQ=255)
                                                    aaTAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTAc              <  1:314213/71‑1 (MQ=255)
                                                      tAGACCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACAt            >  1:3174523/1‑71 (MQ=255)
                                                       aGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATc           >  1:1653941/1‑71 (MQ=255)
                                                         aCCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCg          >  1:2827357/1‑70 (MQ=255)
                                                         aCCATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGt         <  1:1974917/71‑1 (MQ=255)
                                                          ccATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCg          >  1:3567961/1‑69 (MQ=255)
                                                           cATAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGa                         >  1:853118/1‑53 (MQ=255)
                                                           cATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCa                           >  1:2191080/1‑51 (MQ=255)
                                                           cATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAg       >  1:2900946/1‑71 (MQ=255)
                                                            aTAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGc      <  1:2815515/71‑1 (MQ=255)
                                                             tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTAc              <  1:1235998/62‑1 (MQ=255)
                                                             tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGcc     <  1:3748904/71‑1 (MQ=255)
                                                             tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGc      >  1:1507795/1‑70 (MQ=255)
                                                             tAAACTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGc      >  1:648804/1‑70 (MQ=255)
                                                               aaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCg    <  1:425543/70‑1 (MQ=255)
                                                               aaCTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCg    <  1:514266/70‑1 (MQ=255)
                                                                aCTGCAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt  <  1:915132/71‑1 (MQ=255)
                                                                    aaaaaaaaaGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAt  <  1:428783/67‑1 (MQ=255)
                                                                    |                                                                  
CGCCCTTTTTATTTCTGCAATCCGCAATTATCATCCCGCAAATCTACTATGCAATAGACCATAAACTGAAAAAAAAAGTCCGCTGATAAGGCTTGAAAAGTTCATTTCCAGACCCATTTTTACATCGTAGCCGAT  >  minE/2493771‑2493905

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: