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	98.1%	178.4 / NA	53	W69* (TGG→TAG)	xylA	xylose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (20/31);  minor base G (0/1);  total (20/33)
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

GCCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGT  >  NC_000913/3730313‑3730785                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  gcCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACg                                                                                                                                                                                                                                <  2:39110/251‑1 (MQ=255)    gTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCCCATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGAc                                                                                                                                                                                                                             >  1:256108/1‑251 (MQ=255)         gAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGAtt                                                                                                                                                                                                                         >  2:18633/1‑250 (MQ=255)         gAGGGTTTGTAAAGCAGTAGGCCGTCCCCCACAGCAGCTTCATGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGAGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCGCAGGGGAAACATCCACATCGTGGAATCAATAAAATGGCACATGTAACTTGTGGAAAACCGCAAATGCGACAGCTGCTTTACGCTGCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGAttc                                                                                                                                                                                                                        <  2:131492/251‑2 (MQ=255)           gggTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTaa                                                                                                                                                                                                                       >  2:364164/1‑250 (MQ=255)                 gTAAAGCATTTGGCCGTTCCCCACCGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGTAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCGGCTGCCACGGACGAGTAAGCGccc                                                                                                                                                                                                                 <  2:375916/250‑1 (MQ=255)                   aaaGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCca                                                                                                                                                                                                              <  1:465773/251‑1 (MQ=255)                     atcAGTTGGCCGTTCCCCACAGCAGCTTCACTCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGTTCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACAGCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGAGTAAACGCCCCcaca                                                                                                                                                                                                            <  1:216881/249‑1 (MQ=255)                     aGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCcaca                                                                                                                                                                                                            <  1:198132/251‑1 (MQ=255)                       cAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACAcc                                                                                                                                                                                                          <  2:180393/251‑1 (MQ=255)                            ggCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAAc                                                                                                                                                                                                      <  1:448432/250‑1 (MQ=255)                                     ccACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGccc                                                                                                                                                                                            <  1:364164/251‑1 (MQ=255)                                            agcgtCACGCCGCTCTCGTCGTGCTGGCCTGCCAGGACATCAACCAATGGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACAACCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAACTGCGACATCTGCGTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGTTGCTACGGACGATTAAACGCCCCCACACCAACCCTATCCGCCCCGTTCca                                                                                                                                                                                     <  2:23878/247‑1 (MQ=255)                                                  aCGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTGGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGTCATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGaa                                                                                                                                                                               <  1:289538/251‑1 (MQ=255)                                                   cgccgcTCTCTTCTTGCTTGCCTTCCAGGACATCAACCAATTGCGCAAAATTATTGATGTACTCTTTGAACGACGCGCCCTCAGGGGAACCAGCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATAGCCGCCCCGGTCCAGCAGAAg                                                                                                                                                                              <  2:481546/251‑1 (MQ=255)                                                            cttcttGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCACTCAGGGGAAACATCCACATCGTGGAAGCAATAACATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTAAGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCACACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGcc                                                                                                                                                                       >  1:550173/1‑249 (MQ=255)                                                                 tgcttgcCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTa                                                                                                                                                                   <  1:260960/248‑1 (MQ=255)                                                                 tgcttgcCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGGACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGc                                                                                                                                                                 <  2:441999/250‑1 (MQ=255)                                                                  gcttgcCTGCCAGGACATCAACCATGTGCGAAAAATTATTGATGTACTCTTTTAACGACGCGCACGCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAg                                                                                                                                                               <  2:176921/251‑1 (MQ=255)                                                                    ttgcCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAACCGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAggc                                                                                                                                                             >  2:155367/1‑251 (MQ=255)                                                                       cctgccAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCGTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAggcggc                                                                                                                                                          <  2:478124/251‑1 (MQ=255)                                                                        ctgccAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAggcggc                                                                                                                                                          >  1:424572/1‑250 (MQ=255)                                                                                             gcgcAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAgg                                                                                                                                                                             <  1:479962/210‑1 (MQ=255)                                                                                             gcgcAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAgg                                                                                                                                                                             >  2:479962/1‑210 (MQ=255)                                                                                                cAAAATTATTGATGTACGCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGc                                                                                                                                 <  2:428182/251‑1 (MQ=255)                                                                                                 aaaaTTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCt                                                                                                                                >  1:130234/1‑251 (MQ=255)                                                                                                    attattGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTAc                                                                                                                             >  2:57044/1‑251 (MQ=255)                                                                                                      tattGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTACAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTAccc                                                                                                                           >  2:175217/1‑251 (MQ=255)                                                                                                         tGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAac                                                                                                                        >  2:89048/1‑251 (MQ=255)                                                                                                                     ttAACGCCGCGCCCTCAGGGGAAACATCCACCTCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCACCGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCACGTGCTCTTCCATACGCTTACCCACCACCAGTtcgtcg                                                                                                            >  2:229841/1‑251 (MQ=255)                                                                                                                       aaCGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCggg                                                                                                          >  2:125868/1‑251 (MQ=255)                                                                                                                                          aaaCATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCt                                                                                       >  2:99890/1‑251 (MQ=255)                                                                                                                                             cATCCACATCGTGGAAGCAATAAAATGGCACATATAACTTGTGGAAAAACTCAAATGATACAGCTGCTTTACGTTTCGCCAACGCCAGTTCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGTGGCAAAACGCAAGTGCTCTTCCATTCGCTTGTCCAACACCAGTTCGTCTGGATTGGAGTGACTGTATGCTAAc                                                                                    <  2:529150/251‑1 (MQ=255)                                                                                                                                                ccACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACgg                                                                                  <  2:256108/250‑1 (MQ=255)                                                                                                                                                               aaTAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCt                                                                                                                                <  1:152258/189‑1 (MQ=255)                                                                                                                                                               aaTAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCt                                                                                                                                >  2:152258/1‑189 (MQ=255)                                                                                                                                                                  acaaTGGCACATGTAACTTTTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCGACGCCAGTGCCTCCCCACTCTGCTTCTCCCGACGATTAACCGCCCCCACCCCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCt                                                               <  2:177037/249‑1 (MQ=255)                                                                                                                                                                  aaaaTGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATGTTGAGcc                                                                <  1:155367/250‑1 (MQ=255)                                                                                                                                                                     aTGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTc                                                             >  1:18836/1‑250 (MQ=255)                                                                                                                                                                     aTGGCACATGTAACTAGTGGAAAAACTCAAATGCGACATCTGCTGTACGCGTCGCCAACGCCAGTGCCTCACCAGGCTTCTGCGACGGACTATTAAACGCCCCCACACCAAACATATCCGCCCCTGTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTTCATGCGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCa                                                            <  2:49899/251‑1 (MQ=255)                                                                                                                                                                       ggCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATa                                                          <  1:462889/251‑1 (MQ=255)                                                                                                                                                                       ggCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATa                                                          <  1:239456/251‑1 (MQ=255)                                                                                                                                                                        gCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCAAGCAGACGGTGTGCCAGTAGCAGGCGGCAAACCGCAAGTGCTCTTCCATCCGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGCGGATTTTGAGCCTTCATa                                                          >  1:552155/1‑250 (MQ=255)                                                                                                                                                                           caTGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATaacga                                                      <  1:125868/251‑1 (MQ=255)                                                                                                                                                                                 aCTTGTGGAAAAACTCAAATTCGACATCTGCTTTACGCTGCTCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCCTACCCAACACCAGTTCGTCTGGATTGTAGTGCCGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGa                                                <  2:550173/251‑1 (MQ=255)                                                                                                                                                                                          aaaaCTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTgg                                       <  1:89048/251‑1 (MQ=255)                                                                                                                                                                                                       gACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCtt                          <  1:513008/251‑1 (MQ=255)                                                                                                                                                                                                              gCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCACAATAGGCTTGCATAtt                   >  2:163107/1‑251 (MQ=255)                                                                                                                                                                                                                     gCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTcc            >  1:402879/1‑251 (MQ=255)                                                                                                                                                                                                                           ccAACGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATc      <  1:18633/251‑1 (MQ=255)                                                                                                                                                                                                                               cGCCAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACGCCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGt  <  2:156699/251‑1 (MQ=255)                                                                                                                                                                                                                                  cAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATa                                                          >  1:194014/1‑192 (MQ=255)                                                                                                                                                                                                                                  cAGTGCCTCACCAGGCTGCTGCTACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATa                                                          <  2:194014/192‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  GCCGTAGCGAGGGTTTGTAAAGCAGTTGGCCGTTCCCCACAGCAGCTTCACGCCGCTCTCTTCTTGCTTGCCTGCCAGGACATCAACCATTTGCGCAAAATTATTGATGTACTCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACGCCCCCACACCAAACATATCCGCCCCGTTCCAGCAGAAGGTGTGCCAGTAGCAGGCGGCAAAACGCAAGTGCTCTTCCATACGCTTACCCAACACCAGTTCGTCGGGATTGTAGTGACGGAATGCTAACGGGTTTGAGGATTTTGAGCCTTCATAACGAACGCGATCGAGCTGGTCAAAATAGGCTTGCATATTGAACTCCATAATCAGGT  >  NC_000913/3730313‑3730785

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