Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,356,630 C→T 100% T6M (ACG→ATG)  mobA → molybdopterin‑guanine dinucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,356,6300CT77.8% 13.5 / ‑0.1 9T6M (ACG→ATG) mobAmolybdopterin‑guanine dinucleotide synthase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (0/7);  total (0/9)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGC  >  minE/2356562‑2356663
                                                                    |                                 
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:987064/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:108794/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:1093194/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:1133658/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:141542/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:1821512/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:582850/1‑71 (MQ=255)
ttGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGAcgacga                                 >  1:484048/1‑71 (MQ=255)
   agacagacACGTTTGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATaa                              <  1:711549/71‑1 (MQ=255)
                   aGGGTCAATCCCACAATAAAAGAGGCGATATCGGGGAATCTGATGACGACGATaa                              <  1:889603/55‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:1990658/71‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:724084/71‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:2751111/71‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:2296328/71‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:1002177/71‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:178835/71‑1 (MQ=255)
                               acaATAAAAGAGGCGATATCGGTGAATCTGATGACGATGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGc  <  1:1076903/71‑1 (MQ=255)
                                                                    |                                 
TTGAGACAGACACGTTAGCAGGGTCAATCCCACAATAAAAGAGGCGATATCGGTGAATCTGATGACGACGATAACAGGCGTTGTGCTGGCAGGCGGTAAAGC  >  minE/2356562‑2356663

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: ‑