Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,380,292 2 bp→A coding (225‑226/735 nt) frdB ← fumarate reductase iron‑sulfur protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,380,2920G.94.1% 64.9 / ‑0.9 17coding (226/735 nt)frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (12/4);  minor base A (1/0);  total (13/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
*NC_0009134,380,2930TA100.0% 51.0 / NA 16K75N (AAA→AAT) frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (12/4);  total (12/4)

GGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCT  >  NC_000913/4380047‑4380523
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
ggcgTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAAtttgg                                                                                                                                                                                                                                    >  1:110109/1‑249 (MQ=255)
    tAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACg                                                                                                                                                                                                                               >  1:114539/1‑251 (MQ=255)
        acaAACCACAGTTGATGCAACCGGCGAACTGGTGATCCTTCGCCATCTGCGCCGCGGTCTGGATGCTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCATTCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACgttgt                                                                                                                                                                                                                           >  2:46152/1‑251 (MQ=255)
                                         gATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGcc                                                                                                                                                                                           >  1:166499/1‑250 (MQ=255)
                                                                              aCCCTGATCCGCGGTGGGGGAGTTGCCGATGATGTACGGTTTGATCTCTTCCAGCCTTTCGATGAAGTGGGTCATATCGCCCACCAGATCGCGTTCAATCGGGAAGTTAGCGAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCa                                                                                                                                                      <  1:182578/250‑1 (MQ=255)
                                                                                                        cGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCAt                                                                                                                           >  1:143308/1‑251 (MQ=255)
                                                                                                                  cGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGTGTGTCTTTGATGTAGCCCAGCGCATCCAGTAATg                                                                                                                  >  2:61547/1‑250 (MQ=255)
                                                                                                                   ggTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCAACCAGTAATg                                                                                                                  >  2:214588/1‑249 (MQ=255)
                                                                                                                                     cTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGc                                                                                              >  1:126117/1‑251 (MQ=255)
                                                                                                                                                        tatCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTAt                                                                           >  1:210597/1‑251 (MQ=255)
                                                                                                                                                                        gcgTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGAc                                                           <  1:124734/251‑1 (MQ=255)
                                                                                                                                                                                             tAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGcaccac                                      >  1:225363/1‑251 (MQ=255)
                                                                                                                                                                                               aCGCTTCAACCTTCATACCGTCGGTGTAATCACTCGGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATTCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATTTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCCCTTCATAGTATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCt                                    >  2:34051/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCCGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAACGCACCACCTCAATTTTCAGGTTTTTCAtct               >  2:138656/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaTCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCTGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCGTCTCAGCCa         <  1:46841/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGcct  >  1:229965/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCGCAATTTTCAGGTTTTTCATCTCAGCCATTCGcct  <  2:231778/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
GGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCT  >  NC_000913/4380047‑4380523

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑