Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,832,946	T→C	R173R (CGT→CGC)	xthA →	exonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,832,946	0	T	C	100.0%	131.7 / NA	37	R173R (CGT→CGC)	xthA	exonuclease III
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/16);  total (21/16)

AAAAGCGCAGTTTTATCAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGCTCAAAAG  >  NC_000913/1832799‑1833082                                                                                                                                                    |                                                                                                                                         aaaaGCGCAGTTTTATCAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAcc                                                                                                                                       >  2:542228/1‑151 (MQ=255)           ttttATCAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGct                                                                                                                             >  2:59270/1‑151 (MQ=255)            tttATCAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctc                                                                                                                            >  1:260017/1‑151 (MQ=255)             ttATCAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctc                                                                                                                            <  2:277834/150‑1 (MQ=255)                 cAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                                                                       <  1:579914/151‑1 (MQ=255)                    aagcTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGc                                                                                                                    <  1:360106/148‑1 (MQ=255)                         gCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCgg                                                                                                                 <  2:630107/149‑1 (MQ=255)                                    tGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatg                                                                                                    >  1:584466/1‑151 (MQ=255)                                                tCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAgg                                                                                          >  2:590753/1‑149 (MQ=255)                                                   aaCGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgat                                                                                     >  1:246172/1‑151 (MQ=255)                                                     cGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgatga                                                                                   <  2:117079/151‑1 (MQ=255)                                                       tGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAg                                                                                  >  2:392273/1‑150 (MQ=255)                                                              ccGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGt                                                                          <  2:14941/151‑1 (MQ=255)                                                               cGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGtt                                                                         >  1:10832/1‑151 (MQ=255)                                                                ggTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTg                                                                        >  2:283/1‑151 (MQ=255)                                                                   aCTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGt                                                                      <  1:369741/150‑1 (MQ=255)                                                                   aCTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTc                                                                     >  2:44006/1‑151 (MQ=255)                                                                     tGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCg                                                                    >  2:119692/1‑150 (MQ=255)                                                                          aTGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCt                                                              >  2:97586/1‑151 (MQ=255)                                                                          aTGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCt                                                              >  1:300424/1‑151 (MQ=255)                                                                                     gAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCg                                                   <  2:285249/151‑1 (MQ=255)                                                                                           cAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATcc                                              <  2:345200/150‑1 (MQ=255)                                                                                             gCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCa                                           <  2:381465/151‑1 (MQ=255)                                                                                                cTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAAc                                        >  1:276601/1‑151 (MQ=255)                                                                                                  aCAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAAca                                       <  1:7419/150‑1 (MQ=255)                                                                                                    aGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATAGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAAcagc                                     >  1:515796/1‑150 (MQ=255)                                                                                                           gATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCCCCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGAt                                 <  1:328631/147‑1 (MQ=255)                                                                                                              atCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTtc                           >  1:406485/1‑150 (MQ=255)                                                                                                                         cGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGAt               >  2:490469/1‑151 (MQ=255)                                                                                                                          gaagaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGAtt              >  1:272533/1‑151 (MQ=255)                                                                                                                              aaaaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCg          >  1:177048/1‑151 (MQ=255)                                                                                                                               aaaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGc         <  2:194968/151‑1 (MQ=255)                                                                                                                               aaaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGc         >  2:216018/1‑151 (MQ=255)                                                                                                                                  ccGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGCTc       <  1:60425/150‑1 (MQ=255)                                                                                                                                    gTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGCTCaa     >  1:500513/1‑150 (MQ=255)                                                                                                                                      aaGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGCTCAAAAg  <  2:330724/151‑1 (MQ=255)                                                                                                                                      aaGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGCTCAAAAg  <  2:72972/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                         AAAAGCGCAGTTTTATCAGAATCTGCAAAACTACCTGGAAACCGAACTCAAACGTGATAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAATCCGCAAACAGCAGATCGTTTCTCATGGTTTGATTACCGCTCAAAAG  >  NC_000913/1832799‑1833082

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