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%	158.3 / NA	51	*97* (TAG→TAA)	yiaG	HTH_CROC1 family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/30);  total (21/30)

GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGT  >  NC_000913/3719639‑3719894                                                                                                                                  |                                                                                                                               ggTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                       <  2:279983/139‑1 (MQ=255)        ggCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACACGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                <  1:1318088/139‑1 (MQ=255)           gTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                             >  1:258572/1‑139 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      <  2:202544/139‑1 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      <  1:732733/139‑1 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      <  1:22328/139‑1 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      >  1:19274/1‑139 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGACGCGTTTGATTCAAGTCAACCCGACATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      <  2:822970/139‑1 (MQ=255)                    gCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                    <  2:1352703/139‑1 (MQ=255)                         ggTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                               <  2:410685/139‑1 (MQ=255)                         ggTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                               <  2:478518/139‑1 (MQ=255)                                gAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACa                                                                                        <  2:19274/139‑1 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGa                                                                                     >  2:500508/1‑139 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      <  1:1162921/121‑1 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                      >  2:1162921/1‑121 (MQ=255)                                           cAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTc                                                                                                 <  1:860049/119‑1 (MQ=255)                                           cAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTGATTGCTGTTTACGGTc                                                                                                 >  2:860049/1‑119 (MQ=255)                                               cgcgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                         >  2:58175/1‑139 (MQ=255)                                               cgcgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                         <  1:1235435/139‑1 (MQ=255)                                                gcgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                         >  1:1302111/1‑138 (MQ=255)                                                gcgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAAc                                                                         <  2:1302111/138‑1 (MQ=255)                                                  gTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                    <  2:51040/109‑1 (MQ=255)                                                  gTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                    >  1:51040/1‑109 (MQ=255)                                                  gTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                      <  2:258572/139‑1 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                     <  1:500508/139‑1 (MQ=255)                                                       gCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAAt                                                                 <  1:568531/139‑1 (MQ=255)                                                                                 tGCGTTTGATTCAAGCCAACCTGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGAGAGGACCGTTTTCCAACCGAt                                                                     <  1:804982/109‑1 (MQ=255)                                                                                 tGCGTTTGATTCAAGCCAACCTGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                     >  2:804982/1‑109 (MQ=255)                                                                                   cGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                      >  2:1277207/1‑106 (MQ=255)                                                                                   cGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                      <  1:1277207/106‑1 (MQ=255)                                                                                        gATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                    >  2:1204738/1‑87 (MQ=255)                                                                                        gATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                    <  1:1204738/87‑1 (MQ=255)                                                                                              aGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                     <  2:621675/48‑1 (MQ=255)                                                                                              aGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                     >  1:621675/1‑48 (MQ=255)                                                                                                ccAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCg                        <  2:554528/139‑1 (MQ=255)                                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGc                    >  2:502598/1‑139 (MQ=255)                                                                                                         cATTAAGTAAGCAGTTGATGGAATAAGCTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAAt               <  1:871950/139‑1 (MQ=255)                                                                                                            taagtaagCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGca            <  2:650256/139‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTCTCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCg                                                                                  >  2:661490/1‑62 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                    >  2:477133/1‑60 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                    <  1:477133/60‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCg                                                                                  <  1:661490/62‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaataat                                              >  2:542616/1‑98 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaataat                                              <  1:542616/98‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAAc     >  2:704746/1‑139 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAAc     >  1:1348345/1‑139 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAAc     >  2:1108543/1‑139 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGt  >  1:709766/1‑138 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGt  <  2:709766/138‑1 (MQ=255)                                                                                                                          aTGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGc                               >  1:657146/1‑106 (MQ=255)                                                                                                                          aTGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGc                               <  2:657146/106‑1 (MQ=255)                                                                                                                                  |                                                                                                                               GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGT  >  NC_000913/3719639‑3719894

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: ‑