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	96.4%	71.7 / NA	28	D228Y (GAT→TAT)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (13/14);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACGCCACAATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATCGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGAA  >  NC_000913/187836‑188104                                                                                                                                      |                                                                                                                                        aTACGCCACAATATATGCAGACATTCGTTTAATTCAGCCCGCTCCGCTGGGGTTAAGAAGCCAAACCCAACCATATCATCCAGCGAGGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTCGATATAGCGCAAgcc                                                                                                                                >  1:533069/1‑143 (MQ=255)       cacaATATATGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAgg                                                                                                                          >  2:103295/1‑143 (MQ=255)                tGCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCtttt                                                                                                                 >  1:18796/1‑143 (MQ=255)                 gCAGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGGGTGGATCTAGCGCAAGCCGCCAGGGCTGCCTTTg                                                                                                                >  1:180403/1‑143 (MQ=255)                   aGACATTCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTAGAt                                                                                                              >  1:414520/1‑143 (MQ=255)                        ttCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAAGGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGCTGTCTg                                                                                                         <  2:447695/143‑1 (MQ=255)                         tCGTTTAATTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTgg                                                                                                        >  2:640878/1‑143 (MQ=255)                               aaTTCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAg                                                                                                  <  2:558012/143‑1 (MQ=255)                                 ttCCGCCCGCTCCGCTGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGt                                                                                                >  2:321077/1‑143 (MQ=255)                                       ccgctccgctGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGc                                                                                          <  2:518615/143‑1 (MQ=255)                                         gctccgctGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACGCATTGCAGAGTGTGGATATAGCGCAAGGCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTg                                                                                        >  1:625177/1‑143 (MQ=255)                                         gctccgctGAGGTTAAGAAGCCAAACCCGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTg                                                                                        <  2:349858/143‑1 (MQ=255)                                                       gagaagCCAAACCCGACCACTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACg                                                                          <  1:578717/142‑1 (MQ=255)                                                             ccAAACCCGGCCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATg                                                                    <  1:570026/143‑1 (MQ=255)                                                                    cGACCATTTCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGTATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTgg                                                             >  2:689701/1‑143 (MQ=255)                                                                          tttCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgttctgt                                                       >  2:245802/1‑143 (MQ=255)                                                                          tttCATCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGgttctgt                                                       >  2:297292/1‑143 (MQ=255)                                                                              aTCCAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCGGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGttctt                                                   <  2:57197/143‑1 (MQ=255)                                                                                ccAGCGATGTTGCGCCAAAATGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCa                                                 <  1:103295/143‑1 (MQ=255)                                                                                     tgcgTTGCGGCAAAGTCGCGGCGGGCCCCCCATTGCACAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTTGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACttt                                            <  1:628337/140‑1 (MQ=255)                                                                                                 acaTGACGGCGGGCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGaa                                <  1:656648/141‑1 (MQ=255)                                                                                                     gACGGGGGGCCACCCATTGCGGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAACATCAACAATGATTTTTCATTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTg                            <  1:519708/143‑1 (MQ=255)                                                                                                           gggCCACCCATTGCAGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCCGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCTGaa                      <  1:640878/143‑1 (MQ=255)                                                                                                                   ctttGCTGAGTGTGGATATAGCGCAAGCCGCCAGGGCTGCTTCTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGaa              <  2:58494/141‑1 (MQ=255)                                                                                                                         agagTGTGGATATAGCGCAAGCCGCCAGGGCTGCTGTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGATCTTGTCGGAAGGCCAGAATCCTTc        <  1:249672/143‑1 (MQ=255)                                                                                                                          tattgcgcatagCGCGCAGGCCGCCCGGGCTGCTTTTGATGTCTGGTTCCAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCg       <  2:583589/135‑1 (MQ=255)                                                                                                                               tgGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGTTCTGTTCTTCAACTTTCGCCGCGTAGAACTTGTCGGCAGGCCAGAATCCTTCGCTGac  >  2:554947/1‑142 (MQ=255)                                                                                                                               tgGATATAGCGCAAGCCGCCAGGGCTGCTTTTGATGTCTGGTTCAAGGTTGTAGCTGGTGCCATGGTAACGCTGATGGCGCTGGATACGTACGACTGGGTTCGCCGCGTAGAACTTGTCGGAAGGCCAGAATCCTTCGCTGaa  >  2: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 < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑