Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	186,409	C→T	A748T (GCG→ACG)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	186,409	0	C	T	100.0%	76.6 / NA	29	A748T (GCG→ACG)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/10);  total (19/10)

TGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTG  >  NC_000913/186271‑186547                                                                                                                                           |                                                                                                                                           tGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCCTCGTCATa                                                                                                                                        >  1:76751/1‑143 (MQ=255)         tGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACCAGGAGAAGCCATCGGGGTCCAGCACGATAAAGGAATCCATCGTCATACCGtgcga                                                                                                                                >  2:323167/1‑139 (MQ=255)             gACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAgtg                                                                                                                            <  1:323043/143‑1 (MQ=255)                 gggCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGAGGTGa                                                                                                                        >  1:554739/1‑143 (MQ=255)                     tGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAAt                                                                                                                    <  1:477996/143‑1 (MQ=255)                         gTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTtgt                                                                                                                <  1:236841/143‑1 (MQ=255)                                  aGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCCTCGTCATACAGTCGCGAGCCGTGAACATTTGTGCGTCGTGa                                                                                                       >  2:132489/1‑143 (MQ=255)                                     tACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGCCGCGAGTGGTAAAAATTTGTGCGTCGTGAaca                                                                                                    >  2:350778/1‑143 (MQ=255)                                              tCAGTACTTGCTCCAGACCAAACCGAATCACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTg                                                                                           <  2:598875/143‑1 (MQ=255)                                                              actaaaCCGAATAACCTCATGACGATCTTCGTACAGCGTGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCgg                                                                           <  2:256387/140‑1 (MQ=255)                                                                   aCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGGGCGGTGCAAAGGTGTGCGTCGTGGACACTTGAAGTGCGGGGGGCGAATTTGGccgag                                                                      >  2:476239/1‑139 (MQ=255)                                                                             cTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCaa                                                            >  2:651121/1‑143 (MQ=255)                                                                              tCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGAGAGTGGTGAAAATATGTGCGTCGTGAACACTGAAATGGCGGCGGACGAATTCGGCACAGACGGCGGCaac                                                           >  2:178610/1‑142 (MQ=255)                                                                                  gACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGCGAAAATTCGTGCGTCGTGAACACCTAAATTGCGACGGTCTAATTCGGCACAGACGGAGGCAGACAgc                                                       >  2:461771/1‑142 (MQ=255)                                                                                         tGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAACCTTTGTGCGTCGTGAACACTTAAATTGCGGCGGACTAATTCGGCACAGACGGCGGCACACAGACAAgacc                                                >  2:647416/1‑140 (MQ=255)                                                                                                  cGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCCCAGACGGAGGCCGACAGCTACAGGAGGTCCGGGc                                       >  2:523837/1‑143 (MQ=255)                                                                                                  cGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGc                                       <  1:472824/143‑1 (MQ=255)                                                                                                  cGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGc                                       >  2:389952/1‑143 (MQ=255)                                                                                                   gggTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGCTAAGGGCGGTCCGGGCt                                      >  1:399787/1‑143 (MQ=255)                                                                                                      tgtccATGGGGCTCCCGCACGATAAAGGTATCCATCGTCAGACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTACTTCGGCCCAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCa                                   <  2:572375/140‑1 (MQ=255)                                                                                                         ccATCGGGTTCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAACCAGATAAGGGCGGTCCGGGCTCCAGac                                >  1:319849/1‑142 (MQ=255)                                                                                                             cGGGTTCCAGCACGATACAGGTCCACATCGTCATACCGTCGCGCGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCCGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATaaaa                            <  1:248487/143‑1 (MQ=255)                                                                                                              gggTCCCAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATaaaaa                           <  1:147718/143‑1 (MQ=255)                                                                                                                   cccgCACGATAGCGCTATCCAGCGTCATACCGGCGCGAGGTGTGAAAGTTTGTGCGTCGTGAACCCTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGTCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCg                      <  2:651883/140‑1 (MQ=255)                                                                                                                    cAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAACCCGATAAGGGCGCTCCGGGCTCCAGATAAAAATCGCgg                     >  2:532604/1‑143 (MQ=255)                                                                                                                    cAGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAACCAGATAAGGGCGGTCCGGGCTCCCGAGaaaacgctccg                     >  2:610903/1‑136 (MQ=255)                                                                                                                     aGCACGATAAAGGTATCCATCGTCATACCGTCGCGAGAGGGGAGAATTTGTCCGGCGTGAAGACGTCAATTGCGGCGGCCTAAGTCGGCCCTGACAGCGCCAACCACATAAGCGCACTCCGGGCTCCCGAGACACATCCCggg                    >  2:659903/1‑142 (MQ=255)                                                                                                                                      cATCGTCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATGCAAAGCTCGGTGCCCGCACGCGTAGCCt   >  2:181905/1‑143 (MQ=255)                                                                                                                                       aTCGTCATACCGACGCGAGTGGTGAAAATTTGTGCGTCGTGAAGACTTAGATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCCCGCGTAGCCTg  >  1:642724/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                           TGGCGGGTTGGCGACGGGGCTGCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTG  >  NC_000913/186271‑186547

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑