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

CTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCT  >  NC_000913/3199301‑3199557                                                                                                                                 |                                                                                                                                 ctctctTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAt                                                                                                                        >  1:1329454/1‑139 (MQ=255)    tctTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAtgct                                                                                                                     >  1:590274/1‑139 (MQ=255)                aGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAtg                                                                                                                       >  2:1407357/1‑125 (MQ=255)                aGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAtg                                                                                                                       <  1:1407357/125‑1 (MQ=255)                  tGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGcc                                                                                                       >  2:987299/1‑139 (MQ=255)                         gTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAt                                                                                                                        >  2:1246389/1‑115 (MQ=255)                         gTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAt                                                                                                                        <  1:1246389/115‑1 (MQ=255)                                                  tAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTaa                                                                                                                         >  1:1189092/1‑89 (MQ=255)                                                  tAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTaa                                                                                                                         <  2:1189092/89‑1 (MQ=255)                                                        gccgccGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTgg                                                                 <  2:1329454/139‑1 (MQ=255)                                                          cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGAt                                                               >  1:550919/1‑139 (MQ=255)                                                          cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGAt                                                               >  1:1292653/1‑139 (MQ=255)                                                          cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGAt                                                               >  1:619851/1‑139 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:619851/139‑1 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  1:580683/139‑1 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:1181894/139‑1 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:276618/139‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGTAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  1:615304/1‑134 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGTAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:615304/134‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCa                                                                       <  2:1090176/115‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCa                                                                       >  1:1090176/1‑115 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTc                                               <  2:1292653/139‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTc                                               >  1:1174734/1‑139 (MQ=255)                                                                                   ttAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACaaa                                      <  2:1174734/139‑1 (MQ=255)                                                                                      aCCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATc                                   >  1:663253/1‑139 (MQ=255)                                                                                          ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGTCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                               >  2:377058/1‑139 (MQ=255)                                                                                          ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                               <  1:987299/139‑1 (MQ=255)                                                                                          ggCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTa                               <  2:1522383/139‑1 (MQ=255)                                                                                            cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTcc                                                                               <  1:477851/89‑1 (MQ=255)                                                                                            cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTcc                                                                               >  2:477851/1‑89 (MQ=255)                                                                                            cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  1:123661/116‑1 (MQ=255)                                                                                            cTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  2:123661/1‑116 (MQ=255)                                                                                               cGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAaagt                          >  1:164494/1‑139 (MQ=255)                                                                                                    ctcacATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                     <  1:851073/139‑1 (MQ=255)                                                                                                    ctcacATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                     <  2:550919/139‑1 (MQ=255)                                                                                                          ttACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACt               >  2:1344084/1‑139 (MQ=255)                                                                                                           tACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACtt              >  1:535209/1‑139 (MQ=255)                                                                                                               gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  2:1107966/1‑97 (MQ=255)                                                                                                               gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  1:1107966/97‑1 (MQ=255)                                                                                                               gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                     <  1:1463761/128‑1 (MQ=255)                                                                                                               gagCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAagtaagta                     >  2:1463761/1‑128 (MQ=255)                                                                                                                 gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtg        <  2:535209/139‑1 (MQ=255)                                                                                                                  cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATtcgtcg                                                                                                          >  1:1205723/1‑40 (MQ=255)                                                                                                                  cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATtcgtcg                                                                                                          <  2:1205723/40‑1 (MQ=255)                                                                                                                  cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTcc                                              >  2:352471/1‑100 (MQ=255)                                                                                                                  cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTcc                                              <  1:352471/100‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTg                                                                  <  2:1284038/79‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTg                                                                  >  1:1284038/1‑79 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGca                                          >  2:752995/1‑103 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGca                                          <  1:752995/103‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACaaa                                      <  2:1037708/107‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgtg      >  1:1012836/1‑139 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgtg      >  1:715685/1‑139 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgtg      >  1:1534252/1‑139 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTAGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACaaa                                      >  1:1037708/1‑107 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAg                                                                                   <  2:800864/60‑1 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAg                                                                                   >  1:800864/1‑60 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGACCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTgcg    >  1:918842/1‑139 (MQ=255)                                                                                                                       cctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCt  <  2:1012836/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                 CTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCT  >  NC_000913/3199301‑3199557

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