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

CTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCT  >  NC_000913/3199301‑3199557                                                                                                                                 |                                                                                                                                 ctctctTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAAt                                                                                                                        >  1:408788/1‑139 (MQ=255)              ccAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATtcgtc                                                                                                           >  1:703293/1‑139 (MQ=255)              ccAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATtcgtc                                                                                                           >  1:133613/1‑139 (MQ=255)                        cGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTgcggc                                                                                                 >  1:301274/1‑139 (MQ=255)                             gcgcCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGcc                                                                                                                               >  2:783560/1‑104 (MQ=255)                             gcgcCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGcc                                                                                                                               <  1:783560/104‑1 (MQ=255)                              cgcCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                           <  2:1264000/139‑1 (MQ=255)                                  cAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCttt                                                                                       >  1:535860/1‑139 (MQ=255)                                  cAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCttt                                                                                       >  1:547327/1‑139 (MQ=255)                                    ggCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTcc                                                                                     >  1:1039009/1‑139 (MQ=255)                                     gCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCa                                                                                    >  1:1026445/1‑139 (MQ=255)                                        aTGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGACTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGtt                                                                                 <  1:239171/139‑1 (MQ=255)                                             cACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTc                                                                            <  2:1039009/139‑1 (MQ=255)                                              aCCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCa                                                                           <  2:904526/139‑1 (MQ=255)                                               ccATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAg                                                                          >  1:512452/1‑139 (MQ=255)                                                cATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGc                                                                         <  2:673996/139‑1 (MQ=255)                                                cATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGc                                                                         <  2:301274/139‑1 (MQ=255)                                                    aTGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCAct                                                                     <  2:1026445/139‑1 (MQ=255)                                                       cgccgcCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTg                                                                  >  2:69201/1‑139 (MQ=255)                                                          cgccgGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGAt                                                               >  1:56043/1‑139 (MQ=255)                                                                aaTAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCaa                                                         <  2:512452/139‑1 (MQ=255)                                                                aaTAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCaa                                                         <  2:535860/139‑1 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:133613/139‑1 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:56043/139‑1 (MQ=255)                                                                     ccGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:360014/139‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                           >  1:581134/1‑95 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTgg                                                                                           <  2:581134/95‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCAc                                                                      <  2:790145/116‑1 (MQ=255)                                                                          gCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCAc                                                                      >  1:790145/1‑116 (MQ=255)                                                                             cgcgCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTg                                            <  2:408788/139‑1 (MQ=255)                                                                               cggaTTAACCCTGCTTCGTCACTCACACTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  2:875964/4‑129 (MQ=255)                                                                               cgcATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCTTTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  1:938979/1‑129 (MQ=255)                                                                               cgcATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:938979/129‑1 (MQ=255)                                                                               cgcATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGca                                          <  2:547327/139‑1 (MQ=255)                                                                               cgcATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATGCGTCGGCTTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCATTCACGCTGTCCTGca                                          <  2:703293/139‑1 (MQ=255)                                                                               cgcATTAACCCGGCTTCGTCACTCACACTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  1:875964/129‑1 (MQ=255)                                                                                           gCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTc                                               >  2:53825/1‑122 (MQ=255)                                                                                           gCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTc                                               <  1:53825/122‑1 (MQ=255)                                                                                                                 gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCa                                                                                    <  2:1306926/63‑1 (MQ=255)                                                                                                                 gcgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCa                                                                                    >  1:1306926/1‑63 (MQ=255)                                                                                                                  cgcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCtgt       >  2:490023/1‑139 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCa                                                                                    <  1:1013267/61‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCa                                                                                    >  2:1013267/1‑61 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  1:951572/93‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  2:951572/1‑93 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAaagt                          <  2:1014745/119‑1 (MQ=255)                                                                                                                   gcCTCCTGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAaagt                          >  1:1014745/1‑119 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCg                                                                              <  2:636843/65‑1 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCg                                                                              >  1:636843/1‑65 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    <  2:777681/91‑1 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACg                                                    >  1:777681/1‑91 (MQ=255)                                                                                                                     ctcctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTgcg    >  1:334083/1‑139 (MQ=255)                                                                                                                       cctGCAACCACGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCt  <  1:490023/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                 CTCTCTTCAGTAACCAGTGCCAGCGTTTGCGCCCAGGCGATGCGCACCATAATGCGCCGCCGGAATAGCCGCAGCTCGCGCATTAACCCGGCTTCGTCACTCACATTACAGAGCGCCTCCTGCAACCATGCCGCGTAATGCTGCCATTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGCACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCT  >  NC_000913/3199301‑3199557

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