Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,891,836	C→A	C467F (TGT→TTT)	yoaA ←	putative 5' to 3' DNA helicase implicated in DNA repair

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,891,836	0	C	A	100.0%	118.2 / NA	33	C467F (TGT→TTT)	yoaA	putative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/18);  total (15/18)

ACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTT  >  NC_000913/1891761‑1891910                                                                            |                                                                           aCAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGa                                                                       >  1:333638/1‑81 (MQ=255)  cAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGAc                                                                      <  1:333611/81‑1 (MQ=255)        gCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCgttgt                                                                >  1:333642/1‑81 (MQ=255)           tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAg                                                             <  1:333612/81‑1 (MQ=255)               ggTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTc                                                         >  2:333643/1‑81 (MQ=255)                gTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTc                                                         <  2:333614/80‑1 (MQ=255)                gTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTc                                                         <  2:333615/80‑1 (MQ=255)                gTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTc                                                         <  1:333613/80‑1 (MQ=255)                  aGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCgat                                                      >  2:333654/1‑81 (MQ=255)                  aGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCgat                                                      >  2:333651/1‑81 (MQ=255)                   gcgcGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCgatg                                                     <  1:333616/81‑1 (MQ=255)                         aaCTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGt                                               >  2:333655/1‑81 (MQ=255)                             gCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCg                                             >  2:333644/1‑79 (MQ=255)                              cTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCgcag                                          <  1:333618/81‑1 (MQ=255)                              cTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCgc                                            <  2:333617/79‑1 (MQ=255)                                 ggCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCAt                                       >  2:333649/1‑81 (MQ=255)                                  gCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCAtt                                      <  1:333619/81‑1 (MQ=255)                                              atcatGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTgg                          >  2:333666/1‑81 (MQ=255)                                              atcatGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTgg                          <  1:333620/81‑1 (MQ=255)                                                 atGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGc                         >  2:333647/1‑79 (MQ=255)                                                 atGGCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCg                        <  2:333621/80‑1 (MQ=255)                                                    gCGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAg                      >  2:333650/1‑79 (MQ=255)                                                     cGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAg                      >  1:333645/1‑78 (MQ=255)                                                     cGTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAg                   <  1:333622/81‑1 (MQ=255)                                                      gTGCGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGa                  <  2:333623/81‑1 (MQ=255)                                                         ggAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAcc               <  1:333625/80‑1 (MQ=255)                                                         cGAGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAc                <  1:333624/80‑1 (MQ=255)                                                           aGGTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctg             <  2:333626/81‑1 (MQ=255)                                                            ggTACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctgg            >  2:333656/1‑81 (MQ=255)                                                              tACAAAGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggct          <  1:333627/81‑1 (MQ=255)                                                                   aGCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttg     >  1:333659/1‑81 (MQ=255)                                                                    gCATAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttgg    >  1:333657/1‑81 (MQ=255)                                                                      aTAAAAAAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTggttggtt  <  1:333628/81‑1 (MQ=255)                                                                            |                                                                           ACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTT  >  NC_000913/1891761‑1891910

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