Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,719,768	G→A	*97* (TAG→TAA)	yiaG →	HTH_CROC1 family putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,719,768	0	G	A	100.0%	175.2 / NA	55	*97* (TAG→TAA)	yiaG	HTH_CROC1 family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/28);  total (27/28)

GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGA  >  NC_000913/3719639‑3719898                                                                                                                                  |                                                                                                                                   ggTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGCATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCTGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           >  2:292374/1‑139 (MQ=255)   ttttGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGGTTCAAGTCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                         <  1:965143/139‑1 (MQ=255)           gTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAGCTTTTATCCACTTTATTg                                                                                                                 <  2:1100834/139‑1 (MQ=255)            tATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                                >  2:1200077/1‑139 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                          <  2:24646/139‑1 (MQ=255)                       aTGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                 <  1:1060100/127‑1 (MQ=255)                       aTGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAACTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                 >  2:1060100/1‑127 (MQ=255)                       aTGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCC‑ACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTcc                                                                                                    >  1:595096/1‑139 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                    <  1:256289/111‑1 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                    >  2:256289/1‑111 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                          <  2:566931/121‑1 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                          >  1:566931/1‑121 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                        >  1:1339556/1‑139 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                 <  2:129665/113‑1 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                 >  1:129665/1‑113 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                        <  1:1488387/122‑1 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                        >  2:1488387/1‑122 (MQ=255)                                           cAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTc                                                                                 >  1:335353/1‑139 (MQ=255)                                             gACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCa                                                                               <  2:335353/139‑1 (MQ=255)                                                 cgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGttt                                                                                   >  1:783315/1‑131 (MQ=255)                                                 cgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGttt                                                                                   <  2:783315/131‑1 (MQ=255)                                                  gTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                          <  1:1200077/139‑1 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           <  1:1654872/82‑1 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           >  2:1654872/1‑82 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                          >  2:1319947/1‑97 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                          <  1:1319947/97‑1 (MQ=255)                                                                 gCCGAACTAAAATTGATGCGTTTGATTCAAGCC‑ACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGa                                                          <  2:595096/139‑1 (MQ=255)                                                                  ccGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGa                                                          <  2:1339556/139‑1 (MQ=255)                                                                  ccGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGa                                                          >  2:1098682/1‑139 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                 >  1:127175/1‑81 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                 <  2:127175/81‑1 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                                   <  2:995995/95‑1 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                                   >  1:995995/1‑95 (MQ=255)                                                                           aaTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAtt                                                 <  2:1562216/139‑1 (MQ=255)                                                                               gATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgttg                                             >  1:189439/1‑139 (MQ=255)                                                                                aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                                >  1:1417469/1‑71 (MQ=255)                                                                                aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                                <  2:1417469/71‑1 (MQ=255)                                                                                    gTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           >  2:1444633/1‑56 (MQ=255)                                                                                    gTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           <  1:1444633/56‑1 (MQ=255)                                                                                         aTTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGc                                   <  1:1098682/139‑1 (MQ=255)                                                                                               gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           <  2:509476/45‑1 (MQ=255)                                                                                               gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                           >  1:509476/1‑45 (MQ=255)                                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCg                            <  2:600830/135‑1 (MQ=255)                                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCg                            >  1:600830/1‑135 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAtt                                                                        >  1:267348/1‑72 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAtt                                                                        <  2:267348/72‑1 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                         <  1:1465582/71‑1 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                         >  2:1465582/1‑71 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                       >  2:70591/1‑121 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                       <  1:70591/121‑1 (MQ=255)                                                                                                                        tGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAtt                                                 >  1:1581049/1‑94 (MQ=255)                                                                                                                        tGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAtt                                                 <  2:1581049/94‑1 (MQ=255)                                                                                                                        tGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCa                                           <  1:605379/100‑1 (MQ=255)                                                                                                                        tGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCa                                           >  2:605379/1‑100 (MQ=255)                                                                                                                          aTGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGa  >  1:668803/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                   GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGA  >  NC_000913/3719639‑3719898

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