Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,455,846	A→C	14.5%	H228P (CAC→CCC)	paaC →	ring 1,2‑phenylacetyl‑CoA epoxidase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,455,846	0	A	C	14.5%	164.1 / 2.7	73	H228P (CAC→CCC)	paaC	ring 1,2‑phenylacetyl‑CoA epoxidase subunit
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (3/10);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGGGCAAAAGATGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTA  >  NC_000913/1455613‑1456081                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             aTCAGGGCAAAAGATGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AAAACCTTg                                                                                                                                                                                                                                         >  1:68416/1‑239 (MQ=255)  tCAGGGCAAAAGATGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTg                                                                                                                                                                                                                                         >  1:743849/1‑238 (MQ=255)  tCAGGGCAAAAGATGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGGGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGGGAAGGCGGATCGCACGGGCGGGAAAAAAGAAATGCAAACC‑C‑AAAACCATgg                                                                                                                                                                                                                                        >  2:1344082/1‑239 (MQ=255)      ggCAAAAGATGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATGGCACTGGCGGGAAAAAAAGCCTACAAACC‑G‑AACAACTTGGAcc                                                                                                                                                                                                                                     >  1:976878/1‑238 (MQ=255)             aaTGCAGCAGGCCTTCAACAAGTTGGGGCGTTTAA‑CCGGCGAACTGTTCGATGGCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTgg                                                                                                                                                                                                                             <  1:428979/238‑1 (MQ=255)               tGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTgg                                                                                                                                                                                                                             <  1:1344082/237‑1 (MQ=255)                tcacgagGCCATTAAAAAGTTGTGGCGTTTTA‑CCGGCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAg                                                                                                                                                                                                                          <  1:535297/234‑1 (MQ=255)                gcagcaGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGCCGTATCGCACTGGCGGTAAAAAAAGAATCCATACC‑G‑AAAACCTCGGAccccgcgtggcgg                                                                                                                                                                                                                          >  1:1090887/1‑229 (MQ=255)                cggtcagGCCAACAAAATGTTGTGG‑GTTTTACCCGCCGAACTGTTCGATGGCGACGGGATTGATATGGCACTGAGTGAAGAGGGTATTGCGGGTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAg                                                                                                                                                                                                                          <  2:1246384/235‑1 (MQ=255)                    agacagaTCAACAAGTTGGGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAAt                                                                                                                                                                                                                      <  1:580758/233‑1 (MQ=255)                         cATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGGCTGCCTACC‑G‑AACACCCTGGACCCCTGCTGGCAGAAAAGCAGt                                                                                                                                                                                                                 >  1:62849/1‑239 (MQ=255)                                   ttGGGGCTTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAg                                                                                                                                                                                                         <  1:715528/237‑1 (MQ=255)                                                 ccgccgAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGAGCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCg                                                                                                                                                                                            <  2:1159606/239‑1 (MQ=255)                                                 cccccgAACTGTTCGGGTGCGAGGAGATTGATATTGCACTGAGTGAAGGGGGTATTGCGGTTGACCCACGCCATTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑t                                                                                                                                                                                         <  2:743849/236‑1 (MQ=255)                                                   ggcggaCTGTTGGGTGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCa                                                                                                                                                                                       <  2:421153/234‑1 (MQ=255)                                                             tCGATGCCGACGAGATGGATATTGCACTGAGTGAAGAGGGTATTGCGGTGGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTaa                                                                                                                                                                             <  1:1085859/239‑1 (MQ=255)                                                               taTGCCGAGGGGATTGTTATTGCACTGGGTGAGGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACa                                                                                                                                                                           <  2:733135/238‑1 (MQ=255)                                                               gttGCCGGCGAGTTTGGTTTTGCCCTGAGTGAAGGGGGTTTTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACa                                                                                                                                                                           <  2:1374887/237‑1 (MQ=255)                                                               gATGCCGGCGGGTTTGTTTTTTCATTGAGTTAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACa                                                                                                                                                                           <  2:529992/239‑1 (MQ=255)                                                               gATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACa                                                                                                                                                                           <  1:626547/239‑1 (MQ=255)                                                                      acgaGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGGCATGGGAACAaaggaag                                                                                                                                                                    >  1:428630/1‑237 (MQ=255)                                                                       cgaGATTTATATTGCACGGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATTCC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGgaga                                                                                                                                                                    <  1:396584/238‑1 (MQ=255)                                                                        gagaTTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATg                                                                                                                                                                  <  1:1186464/239‑1 (MQ=255)                                                                                tatTGCACTGGGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑AGTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGc                                                                                                                                                           <  1:295536/238‑1 (MQ=255)                                                                                tatTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCaa                                                                                                                                                                                   <  1:1360325/214‑1 (MQ=255)                                                                                     cACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGcc                                                                                                                                                                                              >  1:330413/1‑199 (MQ=255)                                                                                     cACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACCGCCATCC‑G‑ACAAACTTTGACCCATGCTGTACGAAATGCAGGATCCTCAACGGGTCTTTCCCGG‑TCCGCACTGGGAAAAAGGGAAATGGGGATGCAAcc                                                                                                                                                       >  2:295536/1‑236 (MQ=255)                                                                                     cACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGT‑AAAAAGGACTGCATACC‑G‑AACACCTTGGACCCCAGCTGGGAGAAATGCAGATTTTCCCGCGGGTCTTGCCCGG‑TCAGCAATGGGAAAAGAGGGAATGGGGATGCAACg                                                                                                                                                       >  2:1180510/1‑236 (MQ=255)                                                                                     cACTGAGGGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGcc                                                                                                                                                                                              <  2:330413/199‑1 (MQ=255)                                                                                         gAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTACCACAGGAGATGGGTATGCAACGTCTgg                                                                                                                                                  >  1:348470/1‑238 (MQ=255)                                                                                          aGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTACCCGG‑TCAGCAATTGTAACACAGGAACTGGGGTTGCAACGTcacgc                                                                                                                                                 >  1:1169904/1‑234 (MQ=255)                                                                                          aGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAAGGATCGACACC‑A‑CACCCCTTGGGACCATGCTGGGCGAAATACAGAATCTCACCCCGGTCTTTCCCGGCTCA‑CAAAGGTAACAAGGGGAATGGGAAAGCAAAAACTGGcc                                                                                                                                                >  2:1091055/1‑239 (MQ=255)                                                                                          aGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGGAAAAAAGGACTGCATACC‑G‑AACCCCTTTGACCCATGCTGGCAGAAATGCAGGATCTCCAGCGGGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGGATGCACCGGCTGGc                                                                                                                                                 >  1:57741/1‑238 (MQ=255)                                                                                             gAAGAGGGTATGGCGGTTTATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGAATCAACGAAGCCAC‑ATTGAACGTACCGAAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATCCC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGccac                                                                                                                                              <  1:1339353/238‑1 (MQ=255)                                                                                                    gTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATTTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACCGAGGGGATGGGGATGAAACGTCTGGCCACCATTGCa                                                                                                                                       >  1:1083939/1‑238 (MQ=255)                                                                                                    gTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGG‑AGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATTGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTTGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGGGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCAc                                                                                                                                      >  1:1106293/1‑238 (MQ=255)                                                                                                      aTTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTTGACCCATGCTGGCAGAAATGCAGTATCTCCAACGGGTCTTTCCCGG‑TCCGCAATGGTAACAAAGGGACTGGGGATGCAACGTCTGCCCCCCATTGcccc                                                                                                                                     >  1:719974/1‑238 (MQ=255)                                                                                                             ttGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCCTACC‑G‑AACAACTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCACTCGTAAAAGAGGAGGTGGGGAATGAAAGACTGGCCACCCCAGAACCGCAACAgg                                                                                                                             >  2:542730/1‑239 (MQ=255)                                                                                                              gtccccccGCCCTGTACGC‑‑GCGGGGGGGGAGGCCGAATTTTTTGCGGGAAACAACGAAGCCAC‑ATGGAACTGACCCAAAGAGCAGGCTTACCGCCCTGGGGGAAAAAAGGGGCTGCATACC‑G‑AACCCCTGGGACCAATGCTGGCAGAAATGCAGTATCCCCGGGGTGTCTTCCCCGG‑TCAGAAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATGGCACCGCCACAgg                                                                                                                             <  2:1188868/231‑1 (MQ=255)                                                                                                                 ccctcGCACTTTACGC‑‑GAAGGGGGGGAAACCGAAGTTTTTGCCGGGATAAACGAAGCCAC‑ATTGAACGTACCGAAAGGGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑ACCACCTGGGACCCATGCTAGCAGAAATGCAGTATCTCCAGCGGGTCTTGCCCGG‑TAAGCAATGATAAAAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCa                                                                                                                         <  1:405524/235‑1 (MQ=255)                                                                                                                 cccACGCACTTTACGC‑‑GCAGCGTGGAAAGCCGAAGTTTTTGCGGGGACAAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGAAAAAAAGGACTGCATCCC‑G‑AACCCCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGAGGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCa                                                                                                                         <  1:312225/238‑1 (MQ=255)                                                                                                                  ccACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGgaga                                                                                                                                                                    <  1:647595/195‑1 (MQ=255)                                                                                                                            tACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCACACC‑G‑AAACCCCTGGGCCCCTGCCGGCGGAAATTCCGATTTCCCACCGCGGCCTTCCCGG‑CCCGCAAAGGGACACGAGGGGATGGGTATGCCACGCTTGGCCCCCCTGTCACCGCCCCAAGTGCagaggaaaaggg                                                                                                              >  2:312225/1‑228 (MQ=255)                                                                                                                                      gCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCa                                                                                                                                                                                       >  1:25318/1‑158 (MQ=255)                                                                                                                                      gCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCa                                                                                                                                                                                       <  2:25318/158‑1 (MQ=255)                                                                                                                                               gCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAAAGGAGATGGGGATGCAACGTCTGGCCCCCATTGCACCGCCACAGGTTCATGAGGAATGGGGGCTACTCAGCCCGAt                                                                                              >  1:232289/1‑238 (MQ=255)                                                                                                                                                 cGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCCTGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATTGGGTTGCAACGTCTGGCCACCATTGCACCGCCACAAGTTCATGGGATATGGGCGCTCCTCAGCCAGATcac                                                                                           >  1:193091/1‑237 (MQ=255)                                                                                                                                                  gAAGTTTTTGCGGTGATCAACGAAGCCAA‑ATTGACCGTACCTCAAGAGCAGCCGTATCGACCTGGCGGTAAAAAAGGCCTGCATCCC‑G‑AACCCCTGGGACCCATGCGGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGTTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATcccg                                                                                          <  1:1141214/239‑1 (MQ=255)                                                                                                                                                         ttcccGGAATCGACGAA‑CCACAATTGAACGTCCCGCAAGAGCAGGCTTATGGCCTGGGCGGT‑AAAAAGGACTGCATACC‑G‑AACCCCTTGGACCCATGCTGGAAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCCCCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATcccggacccgga                                                                                  <  1:1180510/236‑1 (MQ=255)                                                                                                                                                               ggACCAAGGACGCCAA‑ATAGAACGTACCGCAAGAGGAGGCGTAACAGTCTGGGGGAAAAAAGGCACGGCATCCC‑G‑AACCTCTTGGCCCCATGCGGGAAAAAATGCGGTATCCCGAGCGTGTTTTGCCGGG‑TCAGCAAGGGTAACAGGGGAGATGGGTATGCAAGGTCTGCCCACCTTGGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATcc                                                                             <  2:348470/239‑1 (MQ=255)                                                                                                                                                                gaccAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACGGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTGGGACCCATGCTGGCAAAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAAGGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCCCCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATcc                                                                             <  1:1213989/235‑1 (MQ=255)                                                                                                                                                                          gCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAAGACTGCATACC‑G‑AACACCTTGGACCCATGCTTGCAGAAATGCAGTATCTCCAGCGCGGCTTACCCCG‑ACAGCACTGGGAACAGAGGCGATGGGTATGGAACGTCTGGACCCCCTTGAACCGCCACAGGTTCATGAGATATGGGGGCTACAAAGCCAGATCCCGGACCCAGACAACCCGGGGCTgtaa                                                                  >  2:1047940/1‑236 (MQ=255)                                                                                                                                                                              c‑aTTGAACGTACCGAAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATCCC‑G‑AACCCCTTGGACCCATGCTGGCAGAAAGGCAGTATCCCCAGCGGGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCca                                                                                                                                 <  1:154959/172‑1 (MQ=255)                                                                                                                                                                                    aaCGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATgg                                                                                                                                                                                >  1:760590/1‑120 (MQ=255)                                                                                                                                                                                    aaCGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATgg                                                                                                                                                                                <  2:760590/120‑1 (MQ=255)                                                                                                                                                                                          ccGCAAGAGCAGGCGTGTCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGGAACAGAGGAGATGAGTATGCAACGTCTGGCCACCATTGCACCGCCCCAGGTTCATGAGATATGGGCGCTACCCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTCCTGCTTTAgtc                                                   >  1:894279/1‑237 (MQ=255)                                                                                                                                                                                          ccGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAgaggag                                                                                                                                                                     <  2:1172251/125‑1 (MQ=255)                                                                                                                                                                                          ccGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATAAC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGGAACAgaggag                                                                                                                                                                     >  1:1172251/1‑125 (MQ=255)                                                                                                                                                                                          ccGAAAGAGCAGGCTTATCGCATTGGGGGAAAAAAAGGTCTGCATCCC‑G‑AACCCCTTGGACCCATGCT‑GCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCa                                                  <  1:572744/239‑1 (MQ=255)                                                                                                                                                                                           cGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCCGATCCCGGACCCGGAGATCCCGGGGCGGACCCTTACt                                                            >  1:1083316/1‑229 (MQ=255)                                                                                                                                                                                            gaaagcgaaGGCGTATCGCA‑TGGGAGTAAAAAAAGACCAGAGTCC‑GCAACACTTAGGACACATGCTGGAAAAACGGCATTATCACCAG‑GGGTCTTTCCCGGTTCAGAAAGGGTAGGGGAGGAGATGGTTAGCCAACGTGTGCCCACCTTGGAACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCAt                                                 <  1:1047940/231‑1 (MQ=255)                                                                                                                                                                                            gaaaGGGCGGGCGTTTAGCACGGGCGGTAAAAAGGGACTGCATTCC‑G‑AACCCCTTGGACCCAGGCTGGCAGAAAGGCATTAGCTCCAGCTTGTCTTGCCCGG‑TCAGCAAGGGTAACAGAGGAGATGGTTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCAt                                                 <  1:395784/237‑1 (MQ=255)                                                                                                                                                                                                   cAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCAAGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCCTTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGGGATCCCGGGGCTGGCCATTACTGATTTAGGCATGGTGCGt                                          >  1:960998/1‑239 (MQ=255)                                                                                                                                                                                                          agcgaaCTGGCGG‑AAAAAAGGGCCGCAACCCAG‑ACCCCCTTGGACCCATGCTGGCAGAAATGCGTTATCTCAGGCGTGTCTTGCCCGG‑GCAGCAAGGGTAAAAGAGGAGATGGGTATGCAACGTCTGGCCACAATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGa                                   <  1:1315342/234‑1 (MQ=255)                                                                                                                                                                                                          aTCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGGATGCAACGTCTGGCCACCATTGCACCGCCCCAGGTTCATGAGATATGGGCGCTACTCAGCCCGATCCCCGACCCGGGGATCCCGGTGCTGACCCTTACTGATTTAGGCAAGGTGCGGAAtgtg                                    >  1:908210/1‑238 (MQ=255)                                                                                                                                                                                                           tCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGAACCGGAGATCCCGGTGCTGACCCTTACTGATTTAGGCATGGTCGGTAATGTGaa                                  >  1:633111/1‑238 (MQ=255)                                                                                                                                                                                                            gcacatGGCGGTAAAAAAGGACGGCATCCC‑G‑AACACCGTGGACCCATGCTGGCAGAAATGCAGTTTCTCCGGGGTGTCTTGCCGGG‑TCAGAAATGGTAACAGAGGAGATGGTTATGCACCGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGaca                                 <  1:1091055/234‑1 (MQ=255)                                                                                                                                                                                                                tggggtgTAAAAAAGGACTGCATACG‑G‑AACGCCTTGGACCCATGCTGGAAGAAAGGCATGTCCGCAAGCGTGTCCTGCCGGT‑ACAGCAAGGGTAACAGAGGAAAGGGGTATGACACGCCGGCCCACCATTGCCCCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGa                             <  2:233216/233‑1 (MQ=255)                                                                                                                                                                                                                   gCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGGATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGGCCGTAATGGGACACAAATgg                          >  1:1200520/1‑239 (MQ=255)                                                                                                                                                                                                                                   gCATCCC‑G‑AACCCCTTGGACCCATGCTGGAAGAAATGCAGTATCGCCAGCGTGTCTTGCCCGG‑TCAGAATTGGAAACAGAGGAGATGGGTAGGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGAt          <  1:752886/239‑1 (MQ=255)                                                                                                                                                                                                                                    cATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAAt                                                                                                                                                                                  <  2:334279/70‑1 (MQ=255)                                                                                                                                                                                                                                    cATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAAt                                                                                                                                                                                  >  1:334279/1‑70 (MQ=255)                                                                                                                                                                                                                                           g‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGGGATCCCGGTGCTGACCATTACTGATTTAGGAATGGTGCGTAATGGGACCCAGAGGGGAGAAGGATGGGGGATCGGCttt   >  2:250590/1‑239 (MQ=255)                                                                                                                                                                                                                                              acacCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCAAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGATGCGAAATGTGACACAGCTGGGAGAAGGATGGGGGATCGGCTTTa  >  1:789660/1‑238 (MQ=255)                                                                                                                                                                                                                                                acCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCa                                                                                                                                                                                       <  1:642886/55‑1 (MQ=255)                                                                                                                                                                                                                                                acCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCa                                                                                                                                                                                       >  2:642886/1‑55 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             ATCAGGGCAAAAGATGCAGCAGGCCATCAACAAGTTGTGGCGTTTTA‑CCGCCGAACTGTTCGATGCCGACGAGATTGATATTGCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGC‑‑GCAGCGTGGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCAC‑ATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAAGGACTGCATACC‑G‑AACACCTTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGTCTTGCCCGG‑TCAGCAATGGTAACAGAGGAGATGGGTATGCAACGTCTGGCCACCATTGCACCGCCACAGGTTCATGAGATATGGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTA  >  NC_000913/1455613‑1456081

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