Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,730,560	C→T	W69* (TGG→TAG)	xylA ←	xylose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,730,560	0	C	T	95.3%	128.3 / ‑2.0	43	W69* (TGG→TAG)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (19/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

GCCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTG  >  NC_000913/3730313‑3730798                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               gcCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACg                                                                                                                                                                                                                                             >  1:257108/1‑251 (MQ=255)      aGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACgggcg                                                                                                                                                                                                                                         <  2:34175/250‑4 (MQ=255)                tGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGcc                                                                                                                                                                                                                               <  1:153562/250‑1 (MQ=255)                   aacgCAGTTGGCCGTTCCCCACAGCAGTTTCTCGCCGCTCTCTACTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCGCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAGTGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGccccc                                                                                                                                                                                                                            <  1:225809/247‑1 (MQ=255)                   aaaGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGccccc                                                                                                                                                                                                                            >  2:403348/1‑250 (MQ=255)                              ccGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAAGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACat                                                                                                                                                                                                                 >  1:355910/1‑250 (MQ=255)                                gTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATc                                                                                                                                                                                                              >  2:147810/1‑251 (MQ=255)                                                        ctctctTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGtgtg                                                                                                                                                                                      >  1:429326/1‑251 (MQ=255)                                                                  gcttgcCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCa                                                                                                                                                                             <  1:147810/250‑1 (MQ=255)                                                                    ttgcCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCACCGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACCTATCCGCCCCGTTCCAGCAGAAGGAGTGCCAGTAGCAggc                                                                                                                                                                          >  2:423426/1‑251 (MQ=255)                                                                     tgcctgccAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAgg                                                                                                                                                                           <  1:403348/249‑1 (MQ=255)                                                                      gcctgccAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAggcg                                                                                                                                                                         <  1:422436/250‑1 (MQ=255)                                                                             aGGACATCAACCATTTGCG‑AAAATGTATTGATGTACTCGTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACg                                                                                                                                                                 <  2:144812/251‑1 (MQ=255)                                                                                     aacaaTTTGCGCAAAAT‑TATTGATGTAATCTTTTAACGACGCGCCCTCAGGGGACCCATCCACATCGTGGACGCAATAAAATGGCCCGTGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCCCCAAGCTGCTGCTACGGACGATTACACGCCCCCACAACACACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGct                                                                                                                                                         >  2:63278/5‑251 (MQ=255)                                                                                       ccATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAACCATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGctct                                                                                                                                                       >  1:436500/1‑251 (MQ=255)                                                                                        cATTTGAGAAAAAT‑TATTGATGTACTATTTTAACGACGCGCCCTCAGGGGCAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAACACTCCAATGCAACATCTGCTTTACGCTTCGCCAACACCAGTGCCTCACCTGGCTGCTGCTACGGACGATTAAACGCACCCACACCAAACATATCCGCCCAGTTCCAGCAGAAGGTGTGCCAGTAGCAAGCGGCACACAGCAAGTGATCtt                                                                                                                                                      >  2:393743/1‑251 (MQ=255)                                                                                            tGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCa                                                                                                                                                   <  1:461659/250‑1 (MQ=255)                                                                                               gcAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACg                                                                                                                                               <  1:292925/251‑1 (MQ=255)                                                                                                 aaaaT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCt                                                                                                                                             <  1:289163/251‑1 (MQ=255)                                                                                                  aaaT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCtt                                                                                                                                            <  1:297009/251‑1 (MQ=255)                                                                                                                aCTCGTTTAACGACGCGCCCTCAGGGGAACCATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGt                                                                                                                               <  1:264328/251‑1 (MQ=255)                                                                                                                     tttAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTtc                                                                                                                             <  1:423426/248‑1 (MQ=255)                                                                                                                          cgacgCGCCCTAAGGGGAAACATCCACATCGGGGAAGCAAAAAACTGGCACATGTAACTTGTGGAAAAACTCAAATGCGAAATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGAAGATTAACCGACCCCACACCAAACATATCCCCCCCGTTCTAGCAGAAGGTGTGCAAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCAAATGCTTACCCCACACCAGTTCGTCGGGAt                                                                                                                     >  2:278922/1‑251 (MQ=255)                                                                                                                             cgcgCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCCGAAGGTGTGCCCGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACAACAGTTCGTCGGGATTGt                                                                                                                  >  2:124464/1‑251 (MQ=255)                                                                                                                                   cTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAACCATATCCGCCCCGTTCCAGCGGAAGGTGTGCCAGAAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCAACCCCCAGTTCGTCGGGATTGTAGTGAc                                                                                                            >  2:13911/1‑251 (MQ=255)                                                                                                                                     cAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACg                                                                                                           <  1:13446/250‑1 (MQ=255)                                                                                                                                         gtaaaCATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATg                                                                                                      <  1:428692/249‑1 (MQ=255)                                                                                                                                           aaaCAGCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTTCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTGCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTTCTCTTCCATACGCTTACCCAACACCAGCTCGTCGGGATTGGAGTGACGGAATGCt                                                                                                    <  2:317204/251‑1 (MQ=255)                                                                                                                                             aCATCCACATCGTGGAAGCAATAAAATGGCACATGTACCTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTa                                                                                                   >  1:136620/1‑250 (MQ=255)                                                                                                                                             aCATCCACAGCGGGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTTCTTTACGCTCCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTaa                                                                                                  <  2:429326/251‑1 (MQ=255)                                                                                                                                                 ccACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACggg                                                                                              <  1:102044/251‑1 (MQ=255)                                                                                                                                                     aTCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTg                                                                                          >  1:346141/1‑251 (MQ=255)                                                                                                                                                                aaTAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCa                                                                                                                                                                      >  1:10244/1‑164 (MQ=255)                                                                                                                                                                aaTAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCa                                                                                                                                                                      <  2:10244/164‑1 (MQ=255)                                                                                                                                                                     aaTGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTACCGGGTTTGAGGATTTTGAGCCTTc                                                                          >  1:446523/1‑251 (MQ=255)                                                                                                                                                                      aTGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGGCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTc                                                                          >  2:99758/1‑250 (MQ=255)                                                                                                                                                                                     tgtgGAAAAACTCAAATGCGTCATCTGCTTTACGCTCCGCCAACTCCAGTGCCTCACCAGGCTGCCGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATc                                                           <  2:343876/250‑1 (MQ=255)                                                                                                                                                                                                  aaaTGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAAt                                             <  2:365042/251‑1 (MQ=255)                                                                                                                                                                                                        gACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGCTTGAGGGTTTTGATCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCtt                                       >  1:223319/1‑251 (MQ=255)                                                                                                                                                                                                           attcgCTTTACGCTTCGCCAACGCCAGTCCCTCACCAGGCTGCTGCTACTGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTCTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCa                                    <  2:355910/247‑1 (MQ=255)                                                                                                                                                                                                                              aaCGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAg                 <  1:99758/251‑1 (MQ=255)                                                                                                                                                                                                                                            aCCAGGCTGCTGATACGGACGATTAAAAGACCCAACACCAAACATATCCGCCCCGTTCCAGCAGACGGAGTGCCAGTAGCAGGAGGCAAAACGCAAGTGCTCTTCCATCCGCTTACCCAAAACCAGTTCGTCGGGATTGTAGTGACGGAATGATAACGGGTATGAGGATTTTAAGCCTTCATACAGAACGCGATCGCACTGGTCAAAGTAGGCATGCATATTTAACTCCATAATAAGGTAAGTCCGCGGGt   >  2:170812/1‑251 (MQ=255)                                                                                                                                                                                                                                              cAGGATGCTGCTACGGACGATTAAACGCCACCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTg  >  1:143938/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               GCCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAAT‑TATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGTAATGCCGCGGGTG  >  NC_000913/3730313‑3730798

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