Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,186,012 G→A 100% intergenic (+285/‑220) mtlK → / → RCAP_rcc02028 mannitol 2‑dehydrogenase/FeoA family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,186,0120G→A100.0% 113.6 29intergenic (+285/‑220)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (21/8):  ref base (0/0):  total (21/8)

TGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGAA  >  NC_014034/2185943‑2186109
                                                                       |                                                                                                 
gggCTCAAATTCCCGGGGGGTCCGGGGGCGGAGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACg                                                                    <  1:552240/100‑1 (MQ=9)
                  ggTCCGGGGGGGGAGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGcc                                                  <  1:529876/101‑1 (MQ=16)
                            gggcccccccGGCGGGTCCCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCcaggca                                        <  1:2896637/94‑1 (MQ=25)
                                                  ctcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    >  1:4046579/9‑101 (MQ=255)
                                                  ctcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    <  2:2521457/93‑1 (MQ=255)
                                                   tcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCt                   <  1:2452756/94‑1 (MQ=255)
                                                    cgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCtt                  <  2:3266904/95‑1 (MQ=255)
                                                     gggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTc                 >  1:4065203/6‑101 (MQ=255)
                                                      ggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCg                >  1:3323576/5‑101 (MQ=255)
                                                      ggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCg                <  2:98625/97‑1 (MQ=255)
                                                       gcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCgg               >  1:2497619/4‑101 (MQ=255)
                                                        caaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGa              <  2:59954/99‑1 (MQ=255)
                                                        caaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGa              >  2:2218041/3‑101 (MQ=255)
                                                         aaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGGCGGCCTGTCGCGCAAAAGCGCGCCCAGGCAGGAAGGCTATGCCTTGACCCTTCGGAc             >  1:3617218/2‑101 (MQ=255)
                                                         aaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGAc             >  1:2894322/2‑101 (MQ=255)
                                                         aaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGAc             <  1:1227066/100‑1 (MQ=255)
                                                          aaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGAcc            >  1:3690152/1‑101 (MQ=255)
                                                           agcccGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCg           >  2:195594/4‑101 (MQ=255)
                                                           agcccGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCg           >  2:989138/4‑101 (MQ=255)
                                                           agcccGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCg           >  1:2972313/4‑101 (MQ=255)
                                                              ccGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGt        >  1:4060593/1‑101 (MQ=255)
                                                                gAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTg      >  1:3237857/1‑101 (MQ=255)
                                                                gAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTg      >  2:588472/1‑101 (MQ=255)
                                                                gAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTg      >  2:871455/1‑101 (MQ=255)
                                                                gAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTg      >  1:1495749/1‑101 (MQ=255)
                                                                 aGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTgg     >  2:2188971/1‑101 (MQ=255)
                                                                 aGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCCCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGATCCTTCGTACCGGGtttg     >  1:2068836/1‑99 (MQ=255)
                                                                   gCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGa   >  2:2794282/1‑101 (MQ=255)
                                                                    cGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGaa  >  1:3392649/1‑101 (MQ=255)
                                                                    cGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGaa  >  1:1731583/1‑101 (MQ=255)
                                                                       |                                                                                                 
TGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGAA  >  NC_014034/2185943‑2186109

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG