Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,121,547	G→A	intergenic (‑60/+14)	dinI ← / ← pyrC	DNA damage‑inducible protein I/dihydroorotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,121,547	0	G	A	100.0%	77.5 / NA	26	intergenic (‑60/+14)	dinI/pyrC	DNA damage‑inducible protein I/dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

GTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTG  >  NC_000913/1121321‑1121777                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       gTCGCGCCAATAACCGATAAATTATTCGCTTCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTa                                                                                                                                                                                                                 <  1:401523/250‑1 (MQ=255)      gcCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGAc                                                                                                                                                                                                           <  1:112568/251‑1 (MQ=255)                 aTAAATTAATCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTAGCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGt                                                                                                                                                                                                <  2:54829/251‑1 (MQ=255)                                                                                                     cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAATTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGAt                                                                                                            <  1:77079/251‑1 (MQ=255)                                                                                                         cagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCa                                                                                                                   <  1:314141/240‑1 (MQ=255)                                                                                                         cagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCa                                                                                                                   >  2:314141/1‑240 (MQ=255)                                                                                                          agGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTa                                                                                                                                                                                                                                 <  1:225909/129‑1 (MQ=255)                                                                                                          agGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTa                                                                                                                                                                                                                                 >  2:225909/1‑129 (MQ=255)                                                                                                          agGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGCGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTGTTTATACAGGTACAATAACCTACTGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCCGCGTGTAATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAAtg                                                                                                       >  2:392306/1‑251 (MQ=255)                                                                                                           gggCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACAAGTTCGATGAAtgt                                                                                                      >  1:453963/1‑251 (MQ=255)                                                                                                                            gCTGGCAATGGAGTAGTTTGCGCTATGGTGACTTCAAGTCGTATAATAGCCCCCTTTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAAGGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCGCGTACCAGTTCGATGAATGTGTCGTTGACCGTCAAcc                                                                                     <  2:477518/251‑1 (MQ=255)                                                                                                                                  aaaggAGAAGTTTTCCTTATGGGGACGTCAATTCGCATAATAGCTGCCGGTTGAATATACTGGTTATATATACAGGTAAAATAACCTAAGGACACCAAGAAGCTACGATTGTTATTGTTTAACGGATCAGCGTACCGGTGCCCCGGCAAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGAAGCTTTCATGATTCTTTTGTTCTTCACGGTCGAGTGCGACGAATGTGTCGTTGTCCTGCAAGCCATAGa                                                                               <  2:278441/248‑1 (MQ=255)                                                                                                                                   aTGGAGACGTTTTCGCTATGGTGACTTCAATTCGCCTAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGa                                                                               >  2:245799/1‑250 (MQ=255)                                                                                                                                      gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTg                                                                           <  1:112053/251‑1 (MQ=255)                                                                                                                                        gaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACACCAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCg                                                                         >  2:430570/1‑251 (MQ=255)                                                                                                                                                           gttCAATTCGCAGAATAGCCTCCTGTTGAATATACTCGTTATTTTTATCGGCAACATCACCTCATGCCAACAAGAAGCTACGGTTTTTATTGTTTAACGGACCAGCCTACCGGTTCCCCGGCGAGGACGGGGACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTTTTGATCATCACGTACCAGTTCGCTGAATGTGTCGTTGAACGGCCACACATAGAACTGCGGGCCGTTTacagaataga                                                       >  2:305119/2‑250 (MQ=255)                                                                                                                                                             tCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTAcagaacagaa                                                      <  1:206002/249‑1 (MQ=255)                                                                                                                                                                                         tataCTGGTTATTTATACAGGTAAACTAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCAtttt                        >  1:150610/1‑250 (MQ=255)                                                                                                                                                                                         tataCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAACGTGCTGCAAAGCATTCATttc                        >  1:391350/1‑251 (MQ=255)                                                                                                                                                                                          ataCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATttc                        <  1:134090/250‑1 (MQ=255)                                                                                                                                                                                            aCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATttcttc                     >  1:325330/1‑251 (MQ=255)                                                                                                                                                                                                     ttatACAGGTAAATTAACCTAATGACAACAAGCAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGCTGAATGTGTCGTTGACCGGCAACCCATAGACCTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCCACGCATTCATTATTTCAAAAACGGt            >  2:231433/1‑251 (MQ=255)                                                                                                                                                                                                     ttatACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGt            >  2:381053/1‑251 (MQ=255)                                                                                                                                                                                                      tataCAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGGGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGt            <  2:173016/250‑1 (MQ=255)                                                                                                                                                                                                       ataCAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCAGTCATTTCTTCAAAAACGGTa           <  2:41597/250‑1 (MQ=255)                                                                                                                                                                                                               aaaaTAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTg  >  2:147517/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       GTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTG  >  NC_000913/1121321‑1121777

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