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%	141.2 / NA	40	A71A (GCA→GCG)	glnE	fused deadenylyltransferase/adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/18);  total (22/18)

ATATGCTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTA  >  NC_000913/3199296‑3199558                                                                                                                                      |                                                                                                                                  atatGCTCTCTTCAGTAACCAGTGCCAGCTTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCTTTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCAcgccgc                                                                                                                              <  1:968164/139‑1 (MQ=255)                   ccAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATtcgtc                                                                                                            >  1:255007/1‑139 (MQ=255)                       tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCcgcg                                                                                                                             <  2:663212/118‑1 (MQ=255)                       tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCcgcg                                                                                                                             >  1:663212/1‑118 (MQ=255)                       tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGcc                                                                                                        >  1:512453/1‑139 (MQ=255)                         ccAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTg                                                                                                      <  2:770856/139‑1 (MQ=255)                                 tGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGtt                                                                                              >  2:565397/1‑139 (MQ=255)                                                          tGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGc                                                                          >  2:1222505/1‑134 (MQ=255)                                                          tGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGc                                                                          <  1:1222505/134‑1 (MQ=255)                                                               cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTcc                                                                                <  2:874572/123‑1 (MQ=255)                                                               cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTcc                                                                                >  1:874572/1‑123 (MQ=255)                                                                      aTAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAAt                                                         >  1:564947/1‑139 (MQ=255)                                                                          ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     <  2:564947/139‑1 (MQ=255)                                                                               gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGC‑TTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTcc                                               >  2:163065/1‑139 (MQ=255)                                                                                               ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                                >  1:1063648/1‑139 (MQ=255)                                                                                               ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                                <  1:565397/139‑1 (MQ=255)                                                                                                 cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTaaa                              >  2:98463/1‑139 (MQ=255)                                                                                                   tCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGt                                                                              >  1:754186/1‑89 (MQ=255)                                                                                                   tCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGt                                                                              <  2:754186/89‑1 (MQ=255)                                                                                                     gTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     <  1:457132/112‑1 (MQ=255)                                                                                                     gTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     >  2:457132/1‑112 (MQ=255)                                                                                                     gTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAaagta                          >  1:420/1‑139 (MQ=255)                                                                                                                   agagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTcc                                               <  1:863718/104‑1 (MQ=255)                                                                                                                   agagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTcc                                               >  2:863718/1‑104 (MQ=255)                                                                                                                    gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc           <  2:1063648/139‑1 (MQ=255)                                                                                                                    gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGcc           <  2:512453/139‑1 (MQ=255)                                                                                                                       cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgt        >  1:1355390/1‑139 (MQ=255)                                                                                                                       cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgt        >  2:215796/1‑139 (MQ=255)                                                                                                                        gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     >  1:1069513/1‑93 (MQ=255)                                                                                                                        gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     <  2:1374129/93‑1 (MQ=255)                                                                                                                        gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     <  2:1069513/93‑1 (MQ=255)                                                                                                                        gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                     >  1:1374129/1‑93 (MQ=255)                                                                                                                        gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCt          >  2:1160286/1‑136 (MQ=255)                                                                                                                        gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCt          <  1:1160286/136‑1 (MQ=255)                                                                                                                          ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGTGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTc                                                <  1:976937/96‑1 (MQ=255)                                                                                                                          ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGTGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTc                                                >  2:976937/1‑96 (MQ=255)                                                                                                                          ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCa                                                                            <  2:1319566/68‑1 (MQ=255)                                                                                                                          ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCa                                                                            >  1:1319566/1‑68 (MQ=255)                                                                                                                          ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgtg       >  2:504975/1‑137 (MQ=255)                                                                                                                          ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgtg       <  1:504975/137‑1 (MQ=255)                                                                                                                            cctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGC‑TTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTa  <  1:163065/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                  ATATGCTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTA  >  NC_000913/3199296‑3199558

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