Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	113,695	G→A	W84* (TGG→TGA)	guaC →	GMP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	113,695	0	G	A	97.7%	121.2 / NA	46	W84* (TGG→TGA)	guaC	GMP reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (23/20);  minor base C (1/0);  total (25/21)
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

ATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTG  >  NC_000913/113450‑113939                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aTTGAAGAAGATCTGAAGTTAAGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTGCCGATGTTGAACTGGAACGTCAATTCACCGTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAg                                                                                                                                                                                                                                                  <  2:190805/250‑1 (MQ=255)                            aaaGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGAtg                                                                                                                                                                                                                       <  1:243704/250‑1 (MQ=255)                                          ccGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACAtgtg                                                                                                                                                                                                         <  1:520082/250‑1 (MQ=255)                                                   aCGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCt                                                                                                                                                                                                >  1:141008/1‑250 (MQ=255)                                                    cggtcccctctGAAAGGCCGTGCCGATGGTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGTGCTGGTCCGGCGCTCCGATTATCGCCGCCAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGTGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGGGATGGTTTCTAc                                                                                                                                                                                              <  2:53129/244‑1 (MQ=255)                                                     ccTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCTCCCTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTAcc                                                                                                                                                                                             <  2:151296/250‑1 (MQ=255)                                                        ccACCCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCgg                                                                                                                                                                                           <  1:428373/250‑1 (MQ=255)                                                         cACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTATTCCGCTGATGTGCTGAAAAATGTGCTGGTTTTTACAGGtc                                                                                                                                                                                         >  2:331328/1‑250 (MQ=255)                                                          aCTCTTAAAAGCCGTTCCGATGTTGAACTGAAACGTCACTTCACCTTCAAACGTTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCTCAAATATGAACACCGGAGGCACAGTTTCTATGGCCGCTGCGCTGGCTTCTTTTGATGTTTAGAATGCTGTGCATAAGCACTATTAGGTCGAAGAGAGCCAATCGTTTATCAAAAATTCTTCTGCTGTAGTGCAGAAACATGTGATTGTTTCTACCgatcc                                                                                                                                                                                        >  2:222458/1‑247 (MQ=255)                                                               tAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCt                                                                                                                                                                                    <  1:87713/250‑1 (MQ=255)                                                                    gCCGTTCCGATGGAGAACTGGAACGGCAATTCCCCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTTATGCg                                                                                                                                                                              <  2:268021/251‑1 (MQ=255)                                                                     ccGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCgg                                                                                                                                                                             <  1:485276/251‑1 (MQ=255)                                                                                    aCTGGAACGTCAATTCACCTTCAACCATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGaaacacat                                                                                                                                                               >  2:284118/1‑245 (MQ=255)                                                                                       ggAACGGCAATTCACCTTCAAACATTCAGGTCAAATCTGGGCCGGCGTGCCGATTATCGCCGCAATTATGGACACCGTAGGCACATTTGCTATGGCCTCTGCGCTTGCTTCTTATGATATGTTGACTGCTGGGCATAAACACTAGTCTGTCGAAGAGTGACAAGTGTTTAGCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCTGTACGTCTGAGGCGGATTGCGaaaaacccaaa                                                                                                                                                            <  2:192010/250‑7 (MQ=255)                                                                                           cgtaaaTTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGCTTATCGCCGCAAATATGGACACAGTAGGCACATTTTCAATGGCCTCTGCGCCTGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGCGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTTCTGAAACATGTGATGGTTTCTCCCGGTACGTCTGATGCGGAGTTCGAAAAAACTacaacgat                                                                                                                                                       >  2:306674/5‑244 (MQ=255)                                                                                                             cATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGc                                                                                                                                     >  2:465061/1‑251 (MQ=255)                                                                                                             cATTCAGATCAGAGCGGGTCCGGAGTGCCGATTATCGCCGCAAATTTGGACACCATAGGCACCTTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGTCAAGCGTTTTTCAACAATTCTTCAGCTGATGTGCTGCCACATGTTATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGaccgtgtc                                                                                                                                     >  2:72300/1‑244 (MQ=255)                                                                                                                        gagCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCg                                                                                                                          >  1:284792/1‑251 (MQ=255)                                                                                                                             ggTCCGGCGTTCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCGCGACCTGAACCCGGCATTAAACTTCGttt                                                                                                                       <  2:141315/249‑1 (MQ=255)                                                                                                                                       gCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATTGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTTACGTgg                                                                                                           <  2:473476/251‑1 (MQ=255)                                                                                                                                           aTTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCg                                                                                                                                                                                                                             >  1:262150/1‑133 (MQ=255)                                                                                                                                           aTTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTAGCAACAATTCTTCCg                                                                                                                                                                                                                             <  2:262150/133‑1 (MQ=255)                                                                                                                                                    gcAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTc                                                                                              >  1:259456/1‑251 (MQ=255)                                                                                                                                                              ccacCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATGCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTc                                                                                    <  1:273668/250‑1 (MQ=255)                                                                                                                                                              acacCGTAGGCACATTTTCTCTGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTCACCAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGACTGGTTATTCCGAACACTTc                                                                                    >  2:44451/1‑251 (MQ=255)                                                                                                                                                               aacCGTAGGCACATTTTCTATGGCCTATGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGCAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTAGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCg                                                                                   >  2:504373/2‑251 (MQ=255)                                                                                                                                                                      ggCACATTTTCTATTGCCTCTGCGCTGGCTTCTTTTGATATGTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTTCGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACGTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACTCTTCGTGCAGt                                                                             <  2:401338/250‑1 (MQ=255)                                                                                                                                                                      ggCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGtt                                                                            <  1:465061/251‑1 (MQ=255)                                                                                                                                                                       gCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTc                                                                           <  1:331328/251‑1 (MQ=255)                                                                                                                                                                            ttttCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGCCCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGc                                                                      >  1:562051/1‑251 (MQ=255)                                                                                                                                                                              ttCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCg                                                                     >  1:26288/1‑250 (MQ=255)                                                                                                                                                                                   tGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAgcg                                                               >  2:538/1‑251 (MQ=255)                                                                                                                                                                                                ggCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTg                                                    <  1:284118/249‑1 (MQ=255)                                                                                                                                                                                                ggCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGc                                                  <  1:370387/251‑1 (MQ=255)                                                                                                                                                                                                   ttcttTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCg                                                >  1:99159/1‑250 (MQ=255)                                                                                                                                                                                                    tcttTTGATATTTTGACTGCCGTGCATAAACACTATTCTGTCGACGATTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTATCGACCTGAACCCGGCATTAACCTTCGTTTGTATTGACGTGGCGGATGGTTATTCCGCACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCTTGGCCgcc                                              >  2:340847/1‑249 (MQ=255)                                                                                                                                                                                                      ttttGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTAGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGAccc                                            >  2:183951/1‑250 (MQ=255)                                                                                                                                                                                                        ttGATATTTTGAATGCTGTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAAATAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCAGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCaa                                           >  1:329155/1‑250 (MQ=255)                                                                                                                                                                                                                     tgctgTGCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGAGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTg                             >  2:412524/1‑251 (MQ=255)                                                                                                                                                                                                                      gctgggcttAAATACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAACCCATTTGTGCt                              <  2:407676/241‑1 (MQ=255)                                                                                                                                                                                                                         gtgCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTa                          >  1:141537/1‑250 (MQ=255)                                                                                                                                                                                                                           gCATAAACACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACg                       >  1:350669/1‑251 (MQ=255)                                                                                                                                                                                                                                aaCACTATTCTGTTGAAGAGTGACAAGCGTTTATCAACAACTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATGCTCGACCTGAACCCGGCATTAAACTTCGGTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTg                  <  2:141008/251‑1 (MQ=255)                                                                                                                                                                                                                                aaCACTATTCTGTCGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTg                  >  1:61456/1‑251 (MQ=255)                                                                                                                                                                                                                                             cGAAGAGTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGTTTTAGACAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTAGTATTGACGTGGAGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGGGAAGCGTGGCAGACCAAAACCATTTTTGCTGGTAACGTAGGGAGTGGGGAAAtgt     >  2:304686/1‑251 (MQ=255)                                                                                                                                                                                                                                                  agTGACAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAAtgtgtg  >  1:459841/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ATTGAAGAAGATCTGAAGTTAGGTTTTAAAGACGTTCTCATCCGCCCTAAACGCTCCACTCTTAAAAGCCGTTCCGATGTTGAACTGGAACGTCAATTCACCTTCAAACATTCAGGTCAGAGCTGGTCCGGCGTGCCGATTATCGCCGCAAATATGGACACCGTAGGCACATTTTCTATGGCCTCTGCGCTGGCTTCTTTTGATATTTTGACTGCTGTGCATAAACACTATTCTGTCGAAGAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTGATGCGGATTTCGAAAAAACTAAACAGATTCTCGACCTGAACCCGGCATTAAACTTCGTTTGTATTGACGTGGCGAATGGTTATTCCGAACACTTCGTGCAGTTCGTTGCGAAAGCGCGTGAAGCGTGGCCGACCAAAACCATTTGTGCTGGTAACGTAGTGACTGGTGAAATGTGTG  >  NC_000913/113450‑113939

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