Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	700,038	G→A	H97Y (CAT→TAT)	umpH ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	700,038	0	G	A	97.8%	150.9 / ‑4.1	45	H97Y (CAT→TAT)	umpH	UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/23);  total (21/24)
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.64e-01

GGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCG  >  NC_000913/699791‑700279                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  gggctgggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTGCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCAGCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAAt                                                                                                                                                                                                                                                <  1:426550/251‑1 (MQ=255)  ggctgggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGCGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACAGGCTTTGTAAAGTTCATAAATc                                                                                                                                                                                                                                               >  2:122776/1‑251 (MQ=255)                         gagttGCGCCCGGAGATTTTCGCAACCCCGGCACACTACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTTCGTATTGGTGGCAATAAAACGTGCCCCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCATTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACAGCGGTAATAGTGAAACCGGCTTTGTACTGTTCAGAAATCATTGCGCCTTCGCCCACCACATa                                                                                                                                                                                                                        <  2:419287/248‑1 (MQ=255)                                                         acacaACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAAt                                                                                                                                                                                        <  1:487620/251‑1 (MQ=255)                                                                      caAGCGGGATAAAAACCGCGCCCGTGTGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAACCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCAt                                                                                                                                                                           <  2:19249/251‑1 (MQ=255)                                                                        aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTGTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCg                                                                                                                                                                         <  2:207564/251‑1 (MQ=255)                                                                             gATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAAAGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCag                                                                                                                                                                      >  1:243167/1‑249 (MQ=255)                                                                                cggtctcgcgCCCGTGTGTCTCCGGATTGGGTGCGATCAAACGTGCACCGCTAGCGACGAAATAGGCTGCTTTATGCATCATGGCCCAGCTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGAGAAACCGGCTTTGTACAGTTCAGGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTGGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGta                                                                                                                                                                 <  2:340411/245‑1 (MQ=255)                                                                                           ccGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCt                                                                                                                                                      >  1:110276/1‑251 (MQ=255)                                                                                                ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTGTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAAGCAGGGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGTCAGTCGCCATCGCAGATGTATAAAACACACAGTCAg                                                                                                                                                 <  2:39859/251‑1 (MQ=255)                                                                                                  tgtCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGt                                                                                                                                               >  1:95066/1‑251 (MQ=255)                                                                                                  tgtCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGCTGTAGAAACGCGTTTCGCAAACAATCACGAAATCAGGGTTCCCATCGGTAATAGTGAAACCGGCTTTGTACAGTCCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTACTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGCGGAATAAAACACGCTGTccgg                                                                                                                                                >  1:464037/1‑247 (MQ=255)                                                                                                       ggATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATc                                                                                                                                          <  1:455779/251‑1 (MQ=255)                                                                                                           tggaggCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGACACAGTTGAAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTACTAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATCCGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTAGCCAACGCCGAGGTATAAAACACGCTGTCAGGTAcctagaa                                                                                                                                       >  2:417844/5‑244 (MQ=255)                                                                                                                gCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGc                                                                                                                                 <  1:484035/251‑1 (MQ=255)                                                                                                                  gATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCgg                                                                                                                               >  1:123969/1‑251 (MQ=255)                                                                                                                     aaaaCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTgg                                                                                                                            >  1:174531/1‑251 (MQ=255)                                                                                                                             cACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCgg                                                                                                                    >  1:214874/1‑251 (MQ=255)                                                                                                                              aCCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGt                                                                                                                   >  2:94911/1‑251 (MQ=255)                                                                                                                               gcGTTAGCGACGAAATAGGCTGCTTTATGCAGCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTGCACATCGGTAATAGTGAAACCGGCGTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGTCGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACTCGCGGTCAGGTACCTCGCCTTCTGCGGTGGCAAAGCgg                                                                                                                    <  2:95269/248‑1 (MQ=255)                                                                                                                                  ttAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGTTATAAAACACGCTGTCAGGTACATCGACACCTGACGTGGCAAAGctgatcgc                                                                                                               >  2:12612/1‑244 (MQ=255)                                                                                                                                     gCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAg                                                                                                            <  2:436847/251‑1 (MQ=255)                                                                                                                                          gaAATAGGCTGCGTTAGGCATCATGTCACAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGGTCACATCGGTAATAGTGAAACCGGCTTGGTACAGTTCATAAATCAGTGCGCCTTCGACCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGGACATCGACGCCTGCGGGGGCAAAGCGGTTCGCCAGATCtt                                                                                                       <  2:163921/251‑1 (MQ=255)                                                                                                                                            aaTAGGCTGCTTTATGCAGCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCAGCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACTCCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGc                                                                                                     <  2:51307/251‑1 (MQ=255)                                                                                                                                                    gcattATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCAGGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGGACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTg                                                                                             <  2:158160/248‑1 (MQ=255)                                                                                                                                                         aTGCATCATGTCCCAGTTGTAGGAACGCGGTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTGCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACCCGCTGTCAGGTACATCGACACCTGCGGTGGCAGAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAg                                                                                        <  1:429371/251‑1 (MQ=255)                                                                                                                                                             atcatGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACGCGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATa                                                                                    <  1:315812/251‑1 (MQ=255)                                                                                                                                                              tcatGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGAt                                                                                     >  2:100876/1‑249 (MQ=255)                                                                                                                                                                    cccAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTc                                                                             <  1:168149/251‑1 (MQ=255)                                                                                                                                                                          tGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAac                                                                       >  1:194556/1‑251 (MQ=255)                                                                                                                                                                                  cgcgTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCgg                                                               >  2:111322/1‑251 (MQ=255)                                                                                                                                                                                    cgTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGc                                                              >  1:25001/1‑250 (MQ=255)                                                                                                                                                                                             aacaaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCTTTTAt                                                    <  2:476922/251‑1 (MQ=255)                                                                                                                                                                                                aaTCAAGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCAACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGGCGCCATCGCAGAGGTATAAAACACGCTGTCAGGGACAGCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCa                                                 <  1:101536/251‑1 (MQ=255)                                                                                                                                                                                                   cACGAAATCAGGGTTCACCTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATa                                               >  2:96234/1‑250 (MQ=255)                                                                                                                                                                                                    aCGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGACAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTAAGCAACACCAGCGGCAGGCCTTTATCCATAAt                                             >  1:57119/1‑251 (MQ=255)                                                                                                                                                                                                        aaTCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCCGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCg                                         >  2:172813/1‑251 (MQ=255)                                                                                                                                                                                                         aTCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGt                                        <  1:120249/251‑1 (MQ=255)                                                                                                                                                                                                             ggtttctcaTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGTCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCTACTCCTGCGGTGGCAAAGCGGTTCGCCAGATCTGGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCTTAATCCCGTGCaa                                    <  2:25210/244‑1 (MQ=255)                                                                                                                                                                                                                   acaTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTc                              <  1:23253/251‑1 (MQ=255)                                                                                                                                                                                                                       cGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGc                           <  1:431591/250‑1 (MQ=255)                                                                                                                                                                                                                         gTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGc                        >  2:35265/1‑251 (MQ=255)                                                                                                                                                                                                                                      cGGCTTTGTACAGTTCATAAATACGTGCTCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCCAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGCAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTACGCTGAACCCGGAACGGCg           >  2:406653/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTAt    <  1:175204/251‑1 (MQ=255)                                                                                                                                                                                                                                               aCAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCg  >  2:220320/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCG  >  NC_000913/699791‑700279

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: ‑