Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,266,974	G→A	intergenic (+113/‑140)	alr → / → tyrB	alanine racemase 1/tyrosine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,266,974	0	G	A	98.2%	181.4 / NA	57	intergenic (+113/‑140)	alr/tyrB	alanine racemase 1/tyrosine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/25);  minor base C (1/0);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACAAAGTGAATTTAAGTATCGGTCTG  >  NC_000913/4266729‑4267218                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATTCTCTTATAAAGGTCATGCct                                                                                                                                                                                                                                                 <  1:336001/251‑1 (MQ=255) cAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCct                                                                                                                                                                                                                                                 >  1:129089/1‑251 (MQ=255)           gATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATg                                                                                                                                                                                                                                       >  2:12511/1‑251 (MQ=255)                 gTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGtt                                                                                                                                                                                                                                 >  1:376427/1‑251 (MQ=255)                   cATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAAGGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGGCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGGTGATGTTCTTTCATGCTGTGATAAAGGTCATGCCGGTGGCGGATGTACGttt                                                                                                                                                                                                                                <  2:398221/250‑1 (MQ=255)                        tATGGGTCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTATGGCGGATGTACGTTTGTCATg                                                                                                                                                                                                                          >  1:315643/1‑251 (MQ=255)                        tATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCGAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCAt                                                                                                                                                                                                                           <  1:393530/250‑1 (MQ=255)                         aTGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAACGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTACTCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATg                                                                                                                                                                                                                          >  1:218517/1‑250 (MQ=255)                                 aaGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCAct                                                                                                                                                                                                                 <  1:269025/251‑1 (MQ=255)                                                aaCGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGcc                                                                                                                                                                                                  <  1:346439/251‑1 (MQ=255)                                                   gTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGtt                                                                                                                                                                                               >  1:223215/1‑251 (MQ=255)                                                    tATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAACATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTc                                                                                                                                                                                              >  1:44346/1‑251 (MQ=255)                                                                gaagAAAGTCAGCGCTTACGAAATTATTACGCGCCTGACTTCAAGGATCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCTGGCCAATAACATCCTAAAAATTCAATAAGTTGATATTCTTTCATGCGCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGACTCCCTCTGTTGCTTATTGCCGTTCGCTCCTCAAcct                                                                                                                                                                                  >  1:69508/4‑249 (MQ=255)                                                                  cGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATcc                                                                                                                                                                                <  2:355907/251‑1 (MQ=255)                                                                   gAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCa                                                                                                                                                                               <  1:70159/251‑1 (MQ=255)                                                                       gTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAATCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCg                                                                                                                                                                           >  1:217845/1‑251 (MQ=255)                                                                       gTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTCTGCCGGACGCGGCGTGAATGCCTTATCCGGACAATAAAATCCTAAAAATTCCATAAGTTGATGTTCTTTAATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCg                                                                                                                                                                           >  1:274960/1‑251 (MQ=255)                                                                         aaGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCCAATTGCCGTTCGCTCCTGAACATTCACTCGAt                                                                                                                                                                         >  2:106334/1‑251 (MQ=255)                                                                              cTTACGAACTTATTACGCGTCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATAGTTGATTGTCTGTGACGGATGCGGCGTGAATGCCTTATTCGGCCAATAAAACCCTAAAAATTCAATAAGTTGATGTTCTTGCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTAc                                                                                                                                                                                                                                    <  1:384435/187‑1 (MQ=255)                                                                                 aCGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCct                                                                                                                                                                 >  2:306260/1‑251 (MQ=255)                                                                                  cGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCctt                                                                                                                                                                <  1:243576/251‑1 (MQ=255)                                                                                    aaCTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCcttct                                                                                                                                                              <  2:393900/251‑1 (MQ=255)                                                                                    aaCTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGACTGTGCCGGATGCGGAGTGAATGCATTATCCGGCCAATAAAATACTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGGCATGCCTCTGGCGGCTGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCcttct                                                                                                                                                              >  1:408498/1‑251 (MQ=255)                                                                                                  cTGACTTCACGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGCTGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACAt                                                                                                                                                                                  >  1:331778/1‑217 (MQ=255)                                                                                                  cTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTGCGCTCCTGAACAt                                                                                                                                                                                  <  2:331778/217‑1 (MQ=255)                                                                                                    gACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGtt                                                                                                                                              <  1:20716/251‑1 (MQ=255)                                                                                                    gACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATGTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATtgtg                                                                                                                                                <  2:154763/249‑1 (MQ=255)                                                                                                          aaGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAacc                                                                                                                                        >  1:270247/1‑251 (MQ=255)                                                                                                                cgcgATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGc                                                                                                                                  <  1:399899/251‑1 (MQ=255)                                                                                                                 gcgATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTgac                                                                                                                                 >  2:362192/1‑249 (MQ=255)                                                                                                                   gATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTCTCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTCTAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTTCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCg                                                                                                                               <  2:301292/251‑1 (MQ=255)                                                                                                                          tACGTGGATTAATCGGTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGGCATATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGCTGTTCTTTCATGCTCTTATAAAGGACATGCCTCGGGCGGATGTACGTTTGTCATGAGTCTCACTATGTTGCCAATTGCCGTTCGCTCCTGAACATACAATCGATCTTCGCCTTCTTCCGGTTTTTTGTGGTTTAAACAACTGCCCGGAAACCt                                                                                                                        >  2:130895/1‑251 (MQ=255)                                                                                                                               ggATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTTTGTTTTAACCCCCTGCCCGTAAACCTgg                                                                                                                      >  2:334943/1‑248 (MQ=255)                                                                                                                                gATTAATCGATCCGTAATCTTTGATTGTCTGTGCCGGTTGCGGCGTGAATGCCTTATTCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTTTTTCATGCTCTTCTCAAGGTCATGCCTCTGGCGGATGTACGATTGTCATTAATCTCACTATGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCACCTTCTTCCGGTTTATTGTGTTTTAACAACCTGCCCTTAAACCTGgaaa                                                                                                                   >  2:31081/1‑248 (MQ=255)                                                                                                                                    aaTCGTTCTGTAATATTTGATTGTCTGTGCCGGTTGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCAt                                                                                                              >  1:140708/1‑251 (MQ=255)                                                                                                                                          tCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCgtg                                                                                                        >  1:31667/1‑251 (MQ=255)                                                                                                                                           cTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCgtgt                                                                                                       >  1:312297/1‑251 (MQ=255)                                                                                                                                             gTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCTGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGttt                                                                                                     <  2:270247/251‑1 (MQ=255)                                                                                                                                             gTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCgtgt                                                                                                       <  2:376427/249‑1 (MQ=255)                                                                                                                                                          gTCTTTGCCGGATGCGGCGTGAATGCCTTTTCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACGCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACg                                                                                        <  1:391602/251‑1 (MQ=255)                                                                                                                                                                   ggATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAacc                                                                                                                                        <  2:61077/194‑1 (MQ=255)                                                                                                                                                                   ggATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAacc                                                                                                                                        >  1:61077/1‑194 (MQ=255)                                                                                                                                                                         ggcgTGCATGCCTTATCCGGCCAATAAAATCCTAAAAATTAACTAAGTTGATGTTCTTTCATGCTCTTATAAAAGTAATGCCTCTGACGGATGTACGTTTGTCATGAGTATCACTCTGTTACTAATTGCCGTTCGCTCCTGACCATCCACTCGATCTTCGCCTTCTACCGGTTTATTGCGTTTTAACCACCTGCCCGTAAACATGAAGAACCATCGCGTGTTTACAAACGTTGCCGCCTACGCTGGCGAcc                                                                         >  2:197092/1‑251 (MQ=255)                                                                                                                                                                                 tGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTa                                                                                                                                                     >  1:405101/1‑167 (MQ=255)                                                                                                                                                                                 tGCCTTATCCGGCCAATAAAATCCTAAAAAGTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCGTCTTCCGGTTTa                                                                                                                                                     <  2:405101/167‑1 (MQ=255)                                                                                                                                                                                  gCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGAGCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTa                                                                <  2:218517/251‑1 (MQ=255)                                                                                                                                                                                                  aaaaTCCTCAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCTTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGttt                                                >  1:359793/1‑251 (MQ=255)                                                                                                                                                                                                   acaTCCTAAAAATTCAATAAGTTGATGTTCTTTCCTGCTCTTATAAAGGACATGCCTCTGGCGGCTGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTTACGCCTACGCTGGCGCCCCGATTCTTACGCTTATGGAGCGtttt                                               >  2:353805/3‑251 (MQ=255)                                                                                                                                                                                                               ttCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGACGACCCt                                   >  2:173308/1‑251 (MQ=255)                                                                                                                                                                                                               ttCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCt                                   >  1:206455/1‑251 (MQ=255)                                                                                                                                                                                                                  aaTAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGACGACCCtagc                                >  2:332878/1‑248 (MQ=255)                                                                                                                                                                                                                          gATGTTCTTTCATGCTCTTCTTAAGGTACTGCCTCTGGCGGATGCACGTTTGTCATGAGTCTCACTCTGTTGATAATTGAGGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTac                                                                                                                                          >  2:203706/1‑136 (MQ=255)                                                                                                                                                                                                                          gATGTTCTTTCATGCTCTTATAAAGGTCATGCCGCTGGCGGATGTACGTTTGTCATGAGTCGCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTaa                                                                                                                                          <  1:203706/137‑1 (MQ=255)                                                                                                                                                                                                                          gATGTTCTTTCATGCTCTTAGAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACa                         <  2:129089/250‑1 (MQ=255)                                                                                                                                                                                                                                   tCATGCTCTTATAAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACAAAGTGAAtt                <  2:134963/250‑1 (MQ=255)                                                                                                                                                                                                                                             ataAAGGTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACAAAGTGAATTTAAGTATCgg      <  1:12511/250‑1 (MQ=255)                                                                                                                                                                                                                                                atggTCATGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGATCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACATGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGAGACACGCTTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCCGCGACAAAGTGAATTTAAGTATCGGTCTg  >  2:349613/3‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CAAAGCCGGGGATCCGGTCATTTTATGGGGCGAAGGTTTGCCCGTAGAACGTATCGCTGAAATGACGAAAGTAAGCGCTTACGAACTTATTACGCGCCTGACTTCAAGGGTCGCGATGAAATACGTGGATTAATCGTTCTGTAATATTTGATTGTCTGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCAATAAAATCCTAAAAATTCAATAAGTTGATGTTCTTTCATGCTCTTATAAAGGTCGTGCCTCTGGCGGATGTACGTTTGTCATGAGTCTCACTCTGTTGCTAATTGCCGTTCGCTCCTGAACATCCACTCGATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACAAAGTGAATTTAAGTATCGGTCTG  >  NC_000913/4266729‑4267218

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