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

GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACG  >  NC_000913/3719639‑3719900                                                                                                                                  |                                                                                                                                     ggTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             >  2:285311/1‑139 (MQ=255)            tATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                                  >  2:288326/1‑139 (MQ=255)             aTCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCt                                                                                                                 <  2:915973/139‑1 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                            <  1:1137140/139‑1 (MQ=255)                    gCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                          <  1:1415678/139‑1 (MQ=255)                       aTGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                            >  2:186723/1‑134 (MQ=255)                       aTGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                            <  1:186723/134‑1 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCt                                                                                                                 <  2:1489827/116‑1 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCt                                                                                                                 >  1:1489827/1‑116 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                          >  1:1110563/1‑139 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAc                                                                                          >  2:604960/1‑139 (MQ=255)                                      gAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCg                                                                                        <  1:1192959/139‑1 (MQ=255)                                             gACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             <  1:34133/95‑1 (MQ=255)                                             gACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             >  2:34133/1‑95 (MQ=255)                                              aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATcc                                                                                                                            <  2:71932/95‑1 (MQ=255)                                              aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATcc                                                                                                                            >  1:71932/1‑95 (MQ=255)                                              aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             <  1:1430738/94‑1 (MQ=255)                                              aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             >  2:1430738/1‑94 (MQ=255)                                              aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCaa                                                                                <  2:87807/139‑1 (MQ=255)                                              aCGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCaa                                                                                <  2:1110563/139‑1 (MQ=255)                                                 cgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCg                                                                             >  1:1495815/1‑139 (MQ=255)                                                  gTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                            <  1:604960/139‑1 (MQ=255)                                                  gTGAAGCCTTCAAGTGCCGAACTAAAATTGAAGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                            <  1:574671/139‑1 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                           >  2:1037843/1‑139 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                           >  1:1419882/1‑139 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                           <  1:285311/139‑1 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             <  1:1469268/82‑1 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             >  2:1469268/1‑82 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                            <  2:37911/99‑1 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTa                                                                                                            >  1:37911/1‑99 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCAt                                                                    >  2:1204366/1‑139 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCTACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             >  1:1322705/1‑80 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCGACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                             <  2:1322705/80‑1 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTGTCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATaa                                                                  >  2:13229/1‑139 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                   <  2:888197/90‑1 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                                   >  1:888197/1‑90 (MQ=255)                                                                  ccGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGa                                                            >  1:1491453/1‑139 (MQ=255)                                                                    gAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                            >  2:295415/1‑121 (MQ=255)                                                                    gAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                            <  1:295415/121‑1 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                      >  1:907787/1‑78 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                      <  2:907787/78‑1 (MQ=255)                                                                       cTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAATTATGAAAaat                                                       <  2:282065/139‑1 (MQ=255)                                                                       cTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAaat                                                       <  1:288326/139‑1 (MQ=255)                                                                           aaTTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAtt                                                   <  1:1037843/139‑1 (MQ=255)                                                                                aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                            >  2:933126/1‑109 (MQ=255)                                                                                aTGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                            <  1:933126/109‑1 (MQ=255)                                                                                      ttGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTGTCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAcc                                        <  1:13229/139‑1 (MQ=255)                                                                                                          aTTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                          >  2:1477731/1‑53 (MQ=255)                                                                                                          aTTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                          <  1:1477731/53‑1 (MQ=255)                                                                                                              agtaagCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                          >  1:519315/1‑49 (MQ=255)                                                                                                              agtaagCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                          <  2:519315/49‑1 (MQ=255)                                                                                                              agtaagCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGcaca                >  1:21365/1‑139 (MQ=255)                                                                                                                  agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                                     >  1:686806/1‑50 (MQ=255)                                                                                                                  agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                                     <  2:686806/50‑1 (MQ=255)                                                                                                                  agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTcc                                                                                  <  2:179360/69‑1 (MQ=255)                                                                                                                  agCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTcc                                                                                  >  1:179360/1‑69 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCa                                                                                 <  2:507676/69‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCa                                                                                 >  1:507676/1‑69 (MQ=255)                                                                                                                      gTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGt        <  2:21365/139‑1 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgttg                                               >  2:630482/1‑99 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAAttgttg                                               <  1:630482/99‑1 (MQ=255)                                                                                                                         gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATc                                           <  2:1320941/101‑1 (MQ=255)                                                                                                                         gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATc                                           >  1:1320941/1‑101 (MQ=255)                                                                                                                         gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTg     <  2:1495815/139‑1 (MQ=255)                                                                                                                          aTGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGa    >  2:1316867/1‑139 (MQ=255)                                                                                                                            ggAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACg  <  1:1316867/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGTTTGACG  >  NC_000913/3719639‑3719900

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