Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	540,515	C→T	intergenic (+7/‑50)	allB → / → ybbY	allantoinase/putative uracil/xanthine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	540,515	0	C	T	100.0%	130.0 / NA	41	intergenic (+7/‑50)	allB/ybbY	allantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/19);  total (22/19)

TATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACAC  >  NC_000913/540383‑540638                                                                                                                                     |                                                                                                                            tATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGTCTCCTGCa                                                                                                                       >  1:1079981/1‑139 (MQ=255) tATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCa                                                                                                                       >  2:450323/1‑139 (MQ=255)    gTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGGCATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCTGCAGTAATCTGGCTCCTGCAATg                                                                                                                    <  2:821031/139‑1 (MQ=255)       ggCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGccc                                                                                                                 >  2:860829/1‑139 (MQ=255)                  tGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgcggc                                                                                                      >  1:54251/1‑139 (MQ=255)                     cgcgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCggg                                                                                                   <  1:450323/139‑1 (MQ=255)                      gcgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                  <  2:54251/139‑1 (MQ=255)                                  aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                      <  1:860829/139‑1 (MQ=255)                                                 ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGtt                                                                       >  1:359635/1‑139 (MQ=255)                                                     atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAAt                                                                   >  2:207004/1‑139 (MQ=255)                                                             tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGt                                                           >  2:874860/1‑139 (MQ=255)                                                                    ttGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgc                                                    >  1:676449/1‑139 (MQ=255)                                                                     tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCGTTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                        <  2:797485/134‑1 (MQ=255)                                                                     tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                        >  1:797485/1‑134 (MQ=255)                                                                       aaCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGaa                                                 <  1:754486/139‑1 (MQ=255)                                                                        aCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGaaa                                                <  1:207004/139‑1 (MQ=255)                                                                            ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTg                                                                                                                         >  1:170705/1‑62 (MQ=255)                                                                            ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTg                                                                                                                         <  2:170705/62‑1 (MQ=255)                                                                                    tGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                    <  1:874860/139‑1 (MQ=255)                                                                                      ttGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                  <  2:551153/75‑1 (MQ=255)                                                                                      ttGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                  >  1:551153/1‑75 (MQ=255)                                                                                        gcgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTg                                                                                                                         >  2:593563/1‑50 (MQ=255)                                                                                        gcgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTg                                                                                                                         <  1:593563/50‑1 (MQ=255)                                                                                         cgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                      <  1:153831/84‑1 (MQ=255)                                                                                         cgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                      >  2:153831/1‑84 (MQ=255)                                                                                           ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                      >  2:722217/1‑82 (MQ=255)                                                                                           ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                      <  1:722217/82‑1 (MQ=255)                                                                                           ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                             >  2:1223302/1‑139 (MQ=255)                                                                                           ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                             >  1:873945/1‑139 (MQ=255)                                                                                                     aaTTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAg                            >  2:300423/1‑130 (MQ=255)                                                                                                     aaTTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAg                            <  1:300423/130‑1 (MQ=255)                                                                                                                aaaCATCAGCAGTAATCTGTCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg        <  2:1079981/139‑1 (MQ=255)                                                                                                                aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCg                                                                                                                <  1:735120/35‑1 (MQ=255)                                                                                                                aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCg                                                                                                                >  2:735120/1‑35 (MQ=255)                                                                                                                aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGTTTAAGGTGTGTTTCTGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg        <  2:359635/139‑1 (MQ=255)                                                                                                                      cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                    >  2:380174/1‑105 (MQ=255)                                                                                                                      cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                    <  1:380174/105‑1 (MQ=255)                                                                                                                      cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac  >  2:387674/1‑139 (MQ=255)                                                                                                                      cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac  >  2:1076083/1‑139 (MQ=255)                                                                                                                            tAATCTGGCTCCTGCAATGCCCGTCCTTGcggcgg                                                                                                    >  1:191142/1‑35 (MQ=255)                                                                                                                            tAATCTGGCTCCTGCAATGCCCGTCCTTGcggcgg                                                                                                    <  2:191142/35‑1 (MQ=255)                                                                                                                                     |                                                                                                                            TATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACAC  >  NC_000913/540383‑540638

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑