Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,199,429	T→C	A71A (GCA→GCG)	glnE ←	fused deadenylyltransferase/adenylyltransferase for glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,199,429	0	T	C	100.0%	205.8 / NA	57	A71A (GCA→GCG)	glnE	fused deadenylyltransferase/adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/27);  total (30/27)

GCAATATGCTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCG  >  NC_000913/3199293‑3199555                                                                                                                                         |                                                                                                                               gcAATATGCTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGc                                                                                                                              >  1:145249/1‑139 (MQ=255)                        aGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCTCATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCgg                                                                                                       >  1:85331/1‑139 (MQ=255)                          tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGcc                                                                                                     >  2:1567768/1‑139 (MQ=255)                          tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGcc                                                                                                     >  1:979157/1‑139 (MQ=255)                          tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGcc                                                                                                     >  1:1658163/1‑139 (MQ=255)                          tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGCCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGcc                                                                                                     >  1:438553/1‑139 (MQ=255)                                   ttGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGt                                                                                            <  2:979157/139‑1 (MQ=255)                                   ttGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGt                                                                                            <  1:306720/139‑1 (MQ=255)                                        cccAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCt                                                                                       >  2:881763/1‑139 (MQ=255)                                                  gcgcACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTcc                                                                             <  2:199342/139‑1 (MQ=255)                                                          tAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCa                                                                     <  2:1459029/139‑1 (MQ=255)                                                                  cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACCATACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGAt                                                             >  1:629497/1‑139 (MQ=255)                                                                  cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                         >  2:1292663/1‑111 (MQ=255)                                                                  cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                         <  1:1292663/111‑1 (MQ=255)                                                                  cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGAt                                                             >  1:682700/1‑139 (MQ=255)                                                                    ccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCcgcg                                                                                                                          <  2:555455/76‑1 (MQ=255)                                                                    ccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCcgcg                                                                                                                          >  1:555455/1‑76 (MQ=255)                                                                    ccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGc                                                                       >  1:327558/1‑127 (MQ=255)                                                                    ccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGc                                                                       <  2:327558/127‑1 (MQ=255)                                                                        aaTAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCaa                                                       <  2:85331/139‑1 (MQ=255)                                                                             ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  2:682700/139‑1 (MQ=255)                                                                             ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  2:145249/139‑1 (MQ=255)                                                                             ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  2:438553/139‑1 (MQ=255)                                                                                 aGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGt                                              >  1:284065/1‑139 (MQ=255)                                                                                       cgcATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  >  1:227257/1‑129 (MQ=255)                                                                                       cgcATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  2:227257/129‑1 (MQ=255)                                                                                                  ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                             <  2:1466388/139‑1 (MQ=255)                                                                                                  ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                             <  2:284065/139‑1 (MQ=255)                                                                                                    cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCa                                                                     >  2:1479021/1‑97 (MQ=255)                                                                                                    cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCa                                                                     <  1:1479021/97‑1 (MQ=255)                                                                                                    cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTaaa                           >  2:1415125/1‑139 (MQ=255)                                                                                                            ctcacATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                   <  1:1415125/139‑1 (MQ=255)                                                                                                            ctcacATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                   <  1:1506671/139‑1 (MQ=255)                                                                                                            ctcacATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                   <  2:629497/139‑1 (MQ=255)                                                                                                                 aTTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGcaca                                      <  2:484269/115‑1 (MQ=255)                                                                                                                 aTTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGcaca                                      >  1:484269/1‑115 (MQ=255)                                                                                                                       gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc        >  2:1301731/1‑139 (MQ=255)                                                                                                                       gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc        >  2:1519559/1‑139 (MQ=255)                                                                                                                         gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  2:392939/95‑1 (MQ=255)                                                                                                                         gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  2:1172423/95‑1 (MQ=255)                                                                                                                         gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  >  1:1172423/1‑95 (MQ=255)                                                                                                                         gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  >  1:392939/1‑95 (MQ=255)                                                                                                                         gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtg      <  2:758237/139‑1 (MQ=255)                                                                                                                         gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtg      >  1:758237/1‑139 (MQ=255)                                                                                                                          cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgt     >  2:749597/1‑139 (MQ=255)                                                                                                                           gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                         <  1:507780/54‑1 (MQ=255)                                                                                                                           gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                         >  2:507780/1‑54 (MQ=255)                                                                                                                           gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgtg    >  1:1265221/1‑139 (MQ=255)                                                                                                                            cctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc        <  2:1589238/134‑1 (MQ=255)                                                                                                                            cctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc        >  1:1589238/1‑134 (MQ=255)                                                                                                                             ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCa                                                                         <  2:570247/68‑1 (MQ=255)                                                                                                                             ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCa                                                                         >  1:570247/1‑68 (MQ=255)                                                                                                                             ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  >  2:1092195/1‑91 (MQ=255)                                                                                                                             ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                  <  1:1092195/91‑1 (MQ=255)                                                                                                                             ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTgcg  >  1:560851/1‑139 (MQ=255)                                                                                                                                  gCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTgcg  >  1:30856/1‑134 (MQ=255)                                                                                                                                  gCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTgcg  <  2:30856/134‑1 (MQ=255)                                                                                                                                         |                                                                                                                               GCAATATGCTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCG  >  NC_000913/3199293‑3199555

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