Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	365,306	G→A	R334C (CGC→TGC)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	365,306	0	G	A	84.9%	187.3 / 23.0	73	R334C (CGC→TGC)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base A (31/31);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TGGCTTCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTC  >  NC_000913/365063‑365549                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     tGGCTTCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCTCAGAAAcc                                                                                                                                                                                                                                                <  2:324109/250‑1 (MQ=255) tGGCTTCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCg                                                                                                                                                                                                                                               >  1:190525/1‑251 (MQ=255)   gCTTCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGAc                                                                                                                                                                                                                                             >  2:137289/1‑251 (MQ=255)      tCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACa                                                                                                                                                                                                                                            <  1:41749/249‑1 (MQ=255)      tCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACAt                                                                                                                                                                                                                                           <  2:215342/250‑1 (MQ=255)      gcATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCTGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATc                                                                                                                                                                                                                                          <  1:437708/250‑1 (MQ=255)       cATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCCGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATc                                                                                                                                                                                                                                          <  1:226265/250‑1 (MQ=255)          ccaccaCATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCa                                                                                                                                                                                                                                       <  2:230928/250‑1 (MQ=255)           caccacATACAGGCCGTAGCGGGCGCACAGCGTGTACCACAGCGGATGATTCGGATAATGCGAACAGCGCCCGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATAGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCGTTTTCAATCCGCCCCTTGCATAAACCGACATCGcgtg                                                                                                                                                                                                                                     >  2:93781/1‑248 (MQ=255)                   aCAGGCCGTAGCGGTCGCACCGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTTCACCATCGTCTGCTCATCCATGACCTGACCATGCAGCGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACATAACTCGCAGGCTTCTGCt                                                                                                                                                                                                                             >  2:216962/1‑251 (MQ=255)                        ccGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAAt                                                                                                                                                                                                                        <  1:203183/251‑1 (MQ=255)                            aGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACTGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCGTCTGCTTCAATCAg                                                                                                                                                                                                                     <  2:340330/250‑1 (MQ=255)                                    cacaGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTc                                                                                                                                                                                                            >  1:232462/1‑251 (MQ=255)                                         cGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTcggcgg                                                                                                                                                                                                       <  1:447602/251‑1 (MQ=255)                                          gtgtACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGgt                                                                                                                                                                                                      <  2:81566/251‑1 (MQ=255)                                                  caGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGtt                                                                                                                                                                                              <  2:456923/251‑1 (MQ=255)                                                        aTGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAacc                                                                                                                                                                                         >  2:376717/1‑250 (MQ=255)                                                          ggTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAacc‑ac                                                                                                                                                                                      <  1:15849/250‑1 (MQ=255)                                                              cGGATAATGCGAACAGCGCACGGCGTTAAAGGTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCGCATCCATGACCCGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCATCGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAGCCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGTTGTGCAGTTCAACC‑ACCGCa                                                                                                                                                                                  <  1:343314/250‑1 (MQ=255)                                                                 aTAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGATAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGAt                                                                                                                                                                              >  1:253575/1‑251 (MQ=255)                                                                            aGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCggg                                                                                                                                                                   <  2:399014/251‑1 (MQ=255)                                                                               gcACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGAtt                                                                                                                                                                <  1:287124/251‑1 (MQ=255)                                                                               gcACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATCCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACCTCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCCCGATAGAGATTCGGGAtt                                                                                                                                                                >  1:287556/1‑251 (MQ=255)                                                                                  cGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCg                                                                                                                                                             <  1:222547/251‑1 (MQ=255)                                                                                      gTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGgcgc                                                                                                                                                         <  1:376717/251‑1 (MQ=255)                                                                                        tAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCGCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGGGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTc                                                                                                                                                       >  2:127495/1‑251 (MQ=255)                                                                                             ttgttCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGTATGATGCTCGTGACTGTTAACTCCTCTAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAACCCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCCGCGGTCTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCcaca                                                                                                                                                   <  2:92276/250‑1 (MQ=255)                                                                                                 tCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTc                                                                                                                                              >  1:44019/1‑251 (MQ=255)                                                                                                    gCTTAATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCggg                                                                                                                                           >  1:85836/1‑251 (MQ=255)                                                                                                        catcaGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCTCAGTTTCGGGtttt                                                                                                                                       <  1:141575/251‑1 (MQ=255)                                                                                                        catcaGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGttt                                                                                                                                        <  1:367842/250‑1 (MQ=255)                                                                                                         atcaGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCTACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTc                                                                                                                                      <  1:278962/251‑1 (MQ=255)                                                                                                         atcaGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTc                                                                                                                                      <  1:448844/251‑1 (MQ=255)                                                                                                           cagcagGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCg                                                                                                                                     <  2:232462/250‑1 (MQ=255)                                                                                                             gcagGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACg                                                                                                                                  >  1:239656/1‑251 (MQ=255)                                                                                                                  atatCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTc                                                                                                                               >  1:418836/1‑249 (MQ=255)                                                                                                                         gCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAgt                                                                                                                      <  1:415215/251‑1 (MQ=255)                                                                                                                          cACNATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAgtg                                                                                                                     <  2:442974/251‑1 (MQ=255)                                                                                                                           aCCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAgggt                                                                                                                    >  1:269186/1‑248 (MQ=255)                                                                                                                                      cTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCgg                                                                                                         >  1:263701/1‑251 (MQ=255)                                                                                                                                        cATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCgg                                                                                                         <  1:407781/249‑1 (MQ=255)                                                                                                                                         aTCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCgg                                                                                                                                                            >  2:351905/1‑197 (MQ=255)                                                                                                                                         aTCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCgg                                                                                                                                                            <  1:351905/197‑1 (MQ=255)                                                                                                                                            cATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAac                                                                                                   >  1:102849/1‑251 (MQ=255)                                                                                                                                            cATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAac                                                                                                   >  1:484373/1‑251 (MQ=255)                                                                                                                                            cATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCGCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGAGTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAac                                                                                                   <  2:190525/251‑1 (MQ=255)                                                                                                                                                  ctgaccATGCAGAGGATGATGCTCGTGACGGCTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACcaccac                                                                                             >  2:335716/1‑251 (MQ=255)                                                                                                                                                      ccATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTc                                                                                         >  1:56161/1‑251 (MQ=255)                                                                                                                                                          gCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCg                                                                                     <  2:418836/251‑1 (MQ=255)                                                                                                                                                            agagGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCGTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGAt                                                                                   <  2:61152/251‑1 (MQ=255)                                                                                                                                                              agGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGCCTTCCGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGCGCTTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAAACACCACGCTCATCGATa                                                                                  >  1:292803/1‑250 (MQ=255)                                                                                                                                                              agGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATaa                                                                                 <  2:241954/251‑1 (MQ=255)                                                                                                                                                                gatgatGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCAATCGGCATAACCACCACGCTCATCGATAAtt                                                                               >  1:12429/1‑251 (MQ=255)                                                                                                                                                                         cGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCAccgcc                                                                       <  1:434444/250‑1 (MQ=255)                                                                                                                                                                         cGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCAccgccg                                                                      >  1:263676/1‑251 (MQ=255)                                                                                                                                                                                gTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑CCCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGgcg                                                               <  2:102849/251‑1 (MQ=255)                                                                                                                                                                                gTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGGGAAGCGATCGGCATAAACAACACGCCCATCGATAATTTCACCGCCGAAAGgcg                                                               >  2:8783/1‑251 (MQ=255)                                                                                                                                                                                      gCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGc                                                         >  2:129018/1‑251 (MQ=255)                                                                                                                                                                                                 gCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCGCGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTTCCGCAGGCGAcc                                              >  2:144870/1‑251 (MQ=255)                                                                                                                                                                                                      ggCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCAACGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGt                                         >  2:99383/1‑251 (MQ=255)                                                                                                                                                                                                          tGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTc                                      >  1:118335/1‑250 (MQ=255)                                                                                                                                                                                                                 cagcagcagcagACCATTTTCAATCCGCCCCCCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGGGCCGTCGGCGGTGTGCAGTTCAACC‑CCCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGgc                                                                <  2:423710/217‑1 (MQ=255)                                                                                                                                                                                                                 cagcagcagcagACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGgc                                                                >  1:423710/1‑217 (MQ=255)                                                                                                                                                                                                                   gcagcagcagGCCAGTTTCAACCCGCACCTCGAAGAACCCGACATCGCAGGTTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑CCCGCACGATAGAGATTCGGGATTTCGGCGCTCCACATTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCCCCACGCTCATCGATAATTTCCCCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCAt                            <  1:341683/251‑1 (MQ=255)                                                                                                                                                                                                                   ccagcagcagACCATTTTCATCCCGCCCCTCGCAAAACCCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCACCC‑CCCGCACGATAGAGATTCGGGTTTTCGGCGCTCCACAGTTTCGGTTTTTCGCCGTTCAGACTTATTGTGCCGCGATCGGCATACCCCCCCCGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCAt                            <  1:8783/250‑1 (MQ=255)                                                                                                                                                                                                                        agcagACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGCTTTCGGCGCTCCACAGTTTCGGGGTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATaa                          >  2:488525/1‑248 (MQ=255)                                                                                                                                                                                                                         gcagACCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGa                       >  1:243629/1‑250 (MQ=255)                                                                                                                                                                                                                             aCCATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTg                  >  1:347592/1‑251 (MQ=255)                                                                                                                                                                                                                               cATTTTCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGCTTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGACAGGCGCGGTGCCGATGGAGACCTGCGTTTCACCCTGCCATAAAGAAACTGtt                >  2:470167/1‑251 (MQ=255)                                                                                                                                                                                                                                   ttCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGt                                                                                                                                                 <  2:146573/118‑1 (MQ=255)                                                                                                                                                                                                                                   ttCAATCCGCACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGt                                                                                                                                                 >  1:146573/1‑118 (MQ=255)                                                                                                                                                                                                                                        tcctcACCTCGCAGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGCTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCgtagg       >  1:163315/5‑251 (MQ=255)                                                                                                                                                                                                                                             acctggaggAAACGGACATCACAGGGTTCTGGTTAAAGAAGGGTGCGGGCGGCGGGGTG‑AGGTCAACCCACCGAAAGATAGAGATTCGGGATTGGGGAGCCCCACAGTTGCGGGTTTTCGAGGTGAAGACGTGGTGGGAGGCGATGGGCAAAACCACCCCGCTCATCAAAAATTTCACCGCCGAAAGGCGCGGTGCCGCGGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTc  <  2:149880/244‑1 (MQ=255)                                                                                                                                                                                                                                              ccTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCgtaggtattc  >  1:375133/1‑247 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TGGCTTCATCCACCACATACAGGCCGTAGCGGTCGCACAGCGTGTACCACAGCGGATGGTTCGGATAATGCGAACAGCGCACGGCGTTAAAGTTGTTCTGCTTCATCAGCAGGATATCCTGCACCATCGTCTGCTCATCCATGACCTGACCATGCAGAGGATGATGCTCGTGACGGTTAACGCCTCGAATCAGCAACGGCTTGCCGTTCAGCAGCAGCAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACC‑ACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCCGTAGGTAGTC  >  NC_000913/365063‑365549

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