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% 116.3 33intergenic (+285/‑220)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (23/9):  ref base (1/0):  total (24/9)

GCCGAACCGCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCC.GGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGAA  >  NC_014034/2185922‑2186109
                                                                                             |                                                                                                 
ccgacccccccccTTCCCCCTTGGCTCAAATACCCCGGGGGTCC.GGGGGCGGAGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAgc                                                                                         <  1:4396627/92‑1 (MQ=2)
               tccctttGGCTCAAATACCCGGGGGGTGC.GGGGGCGGAGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACtt                                                                          <  2:4653221/96‑1 (MQ=2)
                             attccccGGGGGGCCGGGGGGCGGGGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTc                                                             <  2:1677032/97‑1 (MQ=9)
                                                   ggcgcccccGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCcaggcag                                       <  2:272140/95‑1 (MQ=21)
                                                          cggccgGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCt                                <  2:828167/99‑1 (MQ=37)
                                                                        ctcgggcaaaGCCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    >  1:2090371/9‑101 (MQ=255)
                                                                        ctcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGGAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    >  1:3243625/9‑101 (MQ=255)
                                                                        ctcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    >  2:1385014/9‑101 (MQ=255)
                                                                        ctcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    >  1:2134597/9‑101 (MQ=255)
                                                                        ctcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGAccc                    >  1:433633/9‑101 (MQ=255)
                                                                         tcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCt                   <  1:3474525/94‑1 (MQ=255)
                                                                         tcgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCt                   >  2:2922719/8‑101 (MQ=255)
                                                                          cgggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCtt                  <  2:888844/95‑1 (MQ=255)
                                                                           gggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTc                 >  1:4209696/6‑101 (MQ=255)
                                                                           gggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTc                 >  1:1695785/6‑101 (MQ=255)
                                                                            ggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCg                >  1:543746/5‑101 (MQ=255)
                                                                            ggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCg                >  2:2873211/5‑101 (MQ=255)
                                                                            ggcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCg                >  2:3148061/5‑101 (MQ=255)
                                                                             gcaaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCgg               <  1:1385109/98‑1 (MQ=255)
                                                                              caaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGa              >  2:3539157/3‑101 (MQ=255)
                                                                              caaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGa              >  1:2475029/3‑101 (MQ=255)
                                                                               aaaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGAc             >  2:1786249/2‑101 (MQ=255)
                                                                                aaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGAcc            >  1:841170/1‑101 (MQ=255)
                                                                                aaGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGAcc            >  2:907658/1‑101 (MQ=255)
                                                                                 agcccGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCg           >  1:331811/4‑101 (MQ=255)
                                                                                     cGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGtt       >  1:3014732/1‑101 (MQ=255)
                                                                                     cGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGCCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCgggct       >  1:1933861/1‑99 (MQ=255)
                                                                                      naGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACNNNNtnn      >  1:695064/2‑99 (MQ=255)
                                                                                       aGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTgg     >  2:3673929/1‑101 (MQ=255)
                                                                                       aGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTgg     <  2:838305/101‑1 (MQ=255)
                                                                                       aGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGGTgg     >  2:3916926/1‑101 (MQ=255)
                                                                                       aGGCGAAACACCCCAGCGACCCCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTgg     <  2:2229571/101‑1 (MQ=255)
                                                                                        ggCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTggg    >  2:3359377/1‑101 (MQ=255)
                                                                                          cGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGaa  >  2:3536019/1‑101 (MQ=255)
                                                                                             |                                                                                                 
GCCGAACCGCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCC.GGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAGGCAGGACGGCTATGCCTTGACCCTTCGGACCGGGTTGGGAA  >  NC_014034/2185922‑2186109

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG