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

GGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAGACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGT  >  NC_000913/3719639‑3719894                                                                                                                                  |                                                                                                                               ggTTTTGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                       <  2:730269/139‑1 (MQ=255)   ttttGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                     <  2:1393249/139‑1 (MQ=255)   ttttGGGCGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                     <  2:1055553/139‑1 (MQ=255)          cGTATCAGTCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTAtt                                                                                                              <  1:178834/139‑1 (MQ=255)                  tCGCCATGGTAAAGGAATGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTCTCCACTTTATTGCTGTTTa                                                                                                      <  2:980359/139‑1 (MQ=255)                    gCCATGGTAAAGGAATGGTAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACg                                                                                                    <  2:988168/139‑1 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                       >  2:1385924/1‑105 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                       <  1:832449/105‑1 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                       <  1:1385924/105‑1 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATc                                                                                                                       >  2:832449/1‑105 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACtt                                                                                                                  <  1:354340/110‑1 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACtt                                                                                                                  >  2:354340/1‑110 (MQ=255)                                   tGGGAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGa                                                                                     >  2:77547/1‑139 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                             >  1:195762/1‑114 (MQ=255)                                    gggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                             <  2:195762/114‑1 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCtgatga                                                                                          <  1:40719/132‑1 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCtgatga                                                                                          >  2:40719/1‑132 (MQ=255)                                     ggAATCCAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGAcc                                                                                   >  2:552516/1‑139 (MQ=255)                                           cAGACGCGTTAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                            <  1:572382/108‑1 (MQ=255)                                           cAGACGCGTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGc                                                                                                            >  2:572382/1‑108 (MQ=255)                                                 cgTGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCg                                                                       >  2:362186/1‑139 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACttt                                                                                                                           <  2:932013/85‑1 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACttt                                                                                                                           >  1:932013/1‑85 (MQ=255)                                                   tGAAGCCTTCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                     >  2:96144/1‑139 (MQ=255)                                                          ttCAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCAt                                                              >  2:614587/1‑139 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                             <  1:661227/90‑1 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTg                                                                                                             >  2:661227/1‑90 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATaa                                                            >  2:903535/1‑139 (MQ=255)                                                            cAAGTGCCGAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATaa                                                            >  1:1257955/1‑139 (MQ=255)                                                                    gAACTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGaaa                                                    <  1:903535/139‑1 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                >  2:495361/1‑78 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTa                                                                                                                <  1:495361/78‑1 (MQ=255)                                                                     aaCTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGaaaa                                                   >  1:922291/1‑139 (MQ=255)                                                                       cTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                               <  2:175020/93‑1 (MQ=255)                                                                       cTAAAATTGATGCGTTTGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCt                                                                                               >  1:175020/1‑93 (MQ=255)                                                                                     tttGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                     >  2:1400700/1‑105 (MQ=255)                                                                                     tttGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGAt                                                                     <  1:1400700/105‑1 (MQ=255)                                                                                     tttGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAc                                   >  1:434930/1‑139 (MQ=255)                                                                                      ttGATTCAAGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCAcc                                  <  1:614587/139‑1 (MQ=255)                                                                                             aaGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                     <  1:1213579/49‑1 (MQ=255)                                                                                             aaGCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCa                                                                                                                     >  2:1213579/1‑49 (MQ=255)                                                                                               gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAat                                                          <  2:954891/106‑1 (MQ=255)                                                                                               gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAat                                                          >  1:954891/1‑106 (MQ=255)                                                                                               gCCAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTTCATCACCCGCCAATGc                         >  2:976057/1‑139 (MQ=255)                                                                                                ccAACCCGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCg                        <  2:1257955/139‑1 (MQ=255)                                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACtt                                                                                                                  <  1:1378592/45‑1 (MQ=255)                                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACtt                                                                                                                  >  2:1378592/1‑45 (MQ=255)                                                                                                    cccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGc                    >  1:859829/1‑139 (MQ=255)                                                                                                     ccGGCATTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCt                   <  2:922291/139‑1 (MQ=255)                                                                                                          aTTAAGTAAGCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATg              <  1:96144/139‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAg                                                                                       >  1:1435740/1‑57 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAg                                                                                       <  2:1435740/57‑1 (MQ=255)                                                                                                                   gCAGTTGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAAc     >  2:427029/1‑139 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGt  >  1:139393/1‑138 (MQ=255)                                                                                                                       ttGATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGATTAATCATAAATATGAAAAATAATTGTTGCATCACCCGCCAATGCGTGGCTTAATGCACATCAACGGt  <  2:139393/138‑1 (MQ=255)                                                                                                                         gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                      >  1:205968/1‑68 (MQ=255)                                                                                                                         gATGGAATAAACTTTTATCCACTTTATTGCTGTTTACGGTCCTGATGACAGGACCGTTTTCCAACCGa                                                                      <  2:205968/68‑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: ‑