Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,298,201 A→T 100% intergenic (‑131/+198) katG ← / ← metF catalase/hydroperoxidase HPI(I)/5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,298,2010AT78.9% 31.9 / 6.9 19intergenic (‑131/+198)katG/metFcatalase/hydroperoxidase HPI(I)/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base T (0/15);  total (1/18)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCACGGATGCGATAGTTGAGATCTTACATATTGTTGGTTAAAGAGATGTAGATCAAATTGATCTTAATTAAGCTTTTTTGATTATGTCCGAATTCGGACAAATGCTTTATAAAAAAGG  >  minE/2298141‑2298258
                                                            |                                                         
ccACGGATGCGATAGTTGAGATCTTACATATTGTTGGTTAAAGAGATGTAGATCAAATTGAt                                                          <  1:64241/62‑1 (MQ=255)
                     tCTTACATATTGTTGGTTAAAGAGATGTAGATCAAATTGATCTTAATTAAGCTTTTTTGAtt                                     >  1:425056/1‑62 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:193609/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:64497/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:48978/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:476514/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:431585/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:425493/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:253394/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:249124/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:222329/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:214601/58‑1 (MQ=255)
                           atatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:175876/58‑1 (MQ=255)
                            tatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:458589/57‑1 (MQ=255)
                            tatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:531577/57‑1 (MQ=255)
                            tatTGTTGGTTAAAGAGATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:89020/57‑1 (MQ=255)
                                         agagATGTAGATCAAATTGTTCTTAATTAAGCTTTTTTGATTAt                                   <  1:510085/44‑1 (MQ=255)
                                           agaTGTAGATCAAATTGATCTTAATTAAGCTTTTTTGATTATGTCCGAATTCGGACAAATGc               <  1:364632/62‑1 (MQ=255)
                                                        aTTGATCTTAATTAAGCTTTTTTGATTATGTCCGAATTCGGACAAATGCTTTATAAAAAAgg  <  1:535672/62‑1 (MQ=255)
                                                            |                                                         
CCACGGATGCGATAGTTGAGATCTTACATATTGTTGGTTAAAGAGATGTAGATCAAATTGATCTTAATTAAGCTTTTTTGATTATGTCCGAATTCGGACAAATGCTTTATAAAAAAGG  >  minE/2298141‑2298258

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