Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	187,969	C→A	D228Y (GAT→TAT)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	187,969	0	C	A	100.0%	72.8 / NA	27	D228Y (GAT→TAT)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/14);  total (13/14)

ATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAA  >  NC_000913/187836‑188104                                                                                                                                      |                                                                                                                                        aTACGCCACAATATATGCAGACATTCGTTTAATTCAGCCCGCTCCGCTGGGGTTAAGAAGCCAAACCCAACCATATCATCCAGCGAGGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTCGATATAGCGCAAgcc                                                                                                                                >  2:533069/1‑143 (MQ=255)       cacaATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAgg                                                                                                                          >  1:103295/1‑143 (MQ=255)                tGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCtttt                                                                                                                 >  2:18796/1‑143 (MQ=255)                 gCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGGGTGGATCTAGCGCAAGCCGCCAGGGCTGCCTTTg                                                                                                                >  2:180403/1‑143 (MQ=255)                   aGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTAGAt                                                                                                              >  2:414520/1‑143 (MQ=255)                        ttCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAAGGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGCTGTCTg                                                                                                         <  1:447695/143‑1 (MQ=255)                         tCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTgg                                                                                                        >  1:640878/1‑143 (MQ=255)                               aaTTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAg                                                                                                  <  1:558012/143‑1 (MQ=255)                                 ttCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGt                                                                                                >  1:321077/1‑143 (MQ=255)                                       ccgctccgctGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGc                                                                                          <  1:518615/143‑1 (MQ=255)                                         gctccgctGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACGCATTGCAGAGTGTGGATATAGCGCAAGGCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTg                                                                                        >  2:625177/1‑143 (MQ=255)                                         gctccgctGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTg                                                                                        <  1:349858/143‑1 (MQ=255)                                                       gagaagCCAAACCCGACCACTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACg                                                                          <  2:578717/142‑1 (MQ=255)                                                             ccAAACCCGGCCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATg                                                                    <  2:570026/143‑1 (MQ=255)                                                                    cGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGTATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTgg                                                             >  1:689701/1‑143 (MQ=255)                                                                          tttCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgttctgt                                                       >  1:297292/1‑143 (MQ=255)                                                                          tttCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgttctgt                                                       >  1:245802/1‑143 (MQ=255)                                                                              aTCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCGGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGttctt                                                   <  1:57197/143‑1 (MQ=255)                                                                                ccAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCa                                                 <  2:103295/143‑1 (MQ=255)                                                                                     tgcgTTGCGGCAAAGTCGCGGCGGGCCCCCCATTGCACAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTTGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACttt                                            <  2:628337/140‑1 (MQ=255)                                                                                                 acaTGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGaa                                <  2:656648/141‑1 (MQ=255)                                                                                                     gACGGGGGGCCACCCATTGCGGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAACATCAACAATGATTTTTCATTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTg                            <  2:519708/143‑1 (MQ=255)                                                                                                           gggCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCCGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCTGaa                      <  2:640878/143‑1 (MQ=255)                                                                                                                   ctttGCTGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTCTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGaa              <  1:58494/141‑1 (MQ=255)                                                                                                                         agagTGTGGATATAGCGCAAGCCGCCAGGGCTGCTGTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGATCTTGTCGGAAGGCCAGAATCCTTc        <  2:249672/143‑1 (MQ=255)                                                                                                                               tgGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGCAGGCCAGAATCCTTCGCTGac  >  1:554947/1‑142 (MQ=255)                                                                                                                               tgGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGATACGTACGACTGGGTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGaa  >  1:165365/1‑143 (MQ=255)                                                                                                                                      |                                                                                                                                        ATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAA  >  NC_000913/187836‑188104

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