Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,729	(G)5→4	coding (54/864 nt)	yicC →	UPF0701 family protein YicC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,725	0	G	.	100.0%	162.8 / NA	39	coding (50/864 nt)	yicC	UPF0701 family protein YicC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/19);  total (20/19)

TCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGC  >  NC_000913/3816645‑3816796                                                                                 |                                                                        tCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATg                                                                         >  2:750830/1‑81 (MQ=255)          ccTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAgcgc                                                                >  2:750826/1‑80 (MQ=255)              ttATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAAcc                                                            <  2:750774/80‑1 (MQ=255)              ttATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAAc                                                             <  2:750773/79‑1 (MQ=255)                 taAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTgg                                                         >  2:750844/1‑80 (MQ=255)                  aaaCAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTggg                                                        >  2:750808/1‑80 (MQ=255)                    aCAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGaaa                                                     >  2:750835/1‑81 (MQ=255)                      aGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATg                                                   <  1:750775/81‑1 (MQ=255)                      aGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATg                                                   >  2:750831/1‑81 (MQ=255)                            gTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCg                                             >  1:750822/1‑81 (MQ=255)                                  gATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAAc                                       >  1:750838/1‑81 (MQ=255)                                   aTCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAAcc                                      <  2:750776/81‑1 (MQ=255)                                      cGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAg                                    <  2:750778/80‑1 (MQ=255)                                      cGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGc                                   <  1:750777/81‑1 (MQ=255)                                        cAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGt                                 <  1:750779/81‑1 (MQ=255)                                         aGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGTAACCTGGGAAATGCGCTCGGTAAACCAGCGt                                 >  1:750853/1‑80 (MQ=255)                                         aGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTTGGTAAACCAGCGt                                 >  2:750852/1‑80 (MQ=255)                                          gTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTa                               >  1:750851/1‑81 (MQ=255)                                           tATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTAt                              <  1:750780/81‑1 (MQ=255)                                           tATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTAt                              >  2:750841/1‑81 (MQ=255)                                           tATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTAt                              <  2:750781/81‑1 (MQ=255)                                             tGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATc                             <  2:750783/80‑1 (MQ=255)                                             tGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATc                             <  1:750782/80‑1 (MQ=255)                                                  gCCTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGaaa                       >  1:750858/1‑81 (MQ=255)                                                    cTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAact                     >  2:750849/1‑81 (MQ=255)                                                    cTACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAact                     <  1:750784/81‑1 (MQ=255)                                                     tACGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAact                     >  1:750850/1‑80 (MQ=255)                                                       cGCCCGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAactt                    <  2:750785/79‑1 (MQ=255)                                                           cGGCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTc              <  2:750786/81‑1 (MQ=255)                                                            ggCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCg             >  2:750839/1‑81 (MQ=255)                                                             gCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGt            <  1:750787/81‑1 (MQ=255)                                                             gCGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGt            >  1:750847/1‑81 (MQ=255)                                                              cGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTc           <  1:750788/81‑1 (MQ=255)                                                              cGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTc           <  1:750789/81‑1 (MQ=255)                                                              cGTGAAATCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTc           >  1:750836/1‑81 (MQ=255)                                                                  aaaTCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGc        <  1:750790/80‑1 (MQ=255)                                                                   aaTCAAGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCg      >  2:750868/1‑81 (MQ=255)                                                                       aaGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAg   >  2:750854/1‑80 (MQ=255)                                                                       aaGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGc  <  2:750791/81‑1 (MQ=255)                                                                        aGGGTGAAT‑GGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAg   >  2:750866/1‑79 (MQ=255)                                                                                 |                                                                        TCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGC  >  NC_000913/3816645‑3816796

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