Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,468,768	C→T	W362* (TGG→TGA)	chiA ←	endochitinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,468,768	0	C	T	97.2%	244.5 / NA	71	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (43/26);  minor base G (0/2);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTT  >  NC_000913/3468533‑3469014                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          cgtGATCGACAGAGGTTTTTTAATGTGCATCATTATAATCTACGATGCTGTCCGCGCTCGACGCTTTGCTGGCCTGGAACCAATCCGCAGCAATGCCGCCGCCGTTATCAATTTCACAGGGAGTCGGGTTACCGTACTGACTCAATGCCGCAGCGGTGGCTTGACGCTCGAGACGCCATGGGTTATTGGAGTTTTCATCGCTGGTACCCTGTGTGCAATGTTTAGGCCCAACCCATCAGGCGTTCTTGTaa                                                                                                                                                                                                                                            <  2:43477/249‑1 (MQ=255)  ttGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTaa                                                                                                                                                                                                                                            >  1:484920/1‑250 (MQ=255)   ggttCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCGGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGGAAAc                                                                                                                                                                                                                                          <  2:7888/248‑1 (MQ=255)      gcGACAGAAGTTTTATAATGTGCAGCGTTAAAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGGCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTGTAGACCCAACCCATCAGGCGTTCTTGTAAACAtt                                                                                                                                                                                                                                       <  2:84285/250‑1 (MQ=255)        gACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTAc                                                                                                                                                                                                                                     <  1:403672/251‑1 (MQ=255)         aCAGAAGTTTTTCAATGTGCATCGTTATAATCGACGATGCTGTCCGCGCGGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATGTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTTCCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTAcc                                                                                                                                                                                                                                    <  2:449533/251‑1 (MQ=255)           aGAAGTTTTATAAGGTGCATCGTTATAATCTTCGATGCGGGCCGCGCAGTACGCTGTGCGGGCCTGGAAACCATCCGCAACAATGACGCCGCCGGGATCAATGTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCTCAGCGGGAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCATCGCTGGTGCCTCGTTGGCAATCTTGAGACCCAACCCATCGGGGGTTCTTGTAAACATTACCAt                                                                                                                                                                                                                                  <  2:82863/251‑1 (MQ=255)              aGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGCCCCAACCCATCACGCGTTCTTGTAAACTTTACCATTg                                                                                                                                                                                                                                >  2:13290/1‑250 (MQ=255)                      aaTGTGCATCGTTATAATCTACGATACTGTCCGCGCTGTAAGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGAACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTAGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGCTCTTGTAAACATTACCATTGAAAATacc                                                                                                                                                                                                                        >  2:385906/1‑248 (MQ=255)                       aTGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACAt                                                                                                                                                                                                                      >  1:176060/1‑251 (MQ=255)                        tgtgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCCTGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATg                                                                                                                                                                                                                     >  2:463924/1‑251 (MQ=255)                          tgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTGACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑t                                                                                                                                                                                                                  >  2:485227/1‑251 (MQ=255)                               cGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAACCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAAt                                                                                                                                                                                                              <  2:181768/250‑1 (MQ=255)                               cGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCATCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGGACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAAtt                                                                                                                                                                                                             <  1:651260/251‑1 (MQ=255)                               cGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTGGCTGGCCTGGAAACCACCCGCAACAATGACGCCTCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAAt                                                                                                                                                                                                              <  2:43924/250‑1 (MQ=255)                                gTTATCATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCCTCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGACACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAAttt                                                                                                                                                                                                            >  1:108529/1‑251 (MQ=255)                                gTTATACTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGATACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCCGCGCTGGTACCCCGTGGGCCATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTCCCATTGAAAATAACATGG‑TAAAttt                                                                                                                                                                                                            >  2:525907/1‑251 (MQ=255)                                   ataATCTACGATGCTTTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCg                                                                                                                                                                                                          <  1:118321/250‑1 (MQ=255)                                        caaCGATTCTGTCCGCGCGGTGCGCTTTGCTGGCCTGAAAACCATCCGCAACAATGATGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAAGTCCGCAGCGGTAGCTGTACGCGCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGAGCAATCTTTAGACCCAACCCAGCAGGCGTGCGTGTAAACATTACCATTGAAAATAACATGG‑GAAATTGCGCTAcc                                                                                                                                                                                                     <  2:329850/248‑1 (MQ=255)                                        cTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAACATAACATGG‑TAAATTTCGCTACCt                                                                                                                                                                                                    >  1:376557/1‑251 (MQ=255)                                          aCGATGCGGTCCGCTCTTTACGCTTTGCTGGCCGGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAAGTCCTCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGTCAAGCTGTAGACCCAACCCATCAGTCGTTCTTGTAAACATTACCATTGTAAATAACATGG‑TAAAGTTCGCTACCtt                                                                                                                                                                                                   <  2:650971/250‑1 (MQ=255)                                                       cgcTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCAc                                                                                                                                                                                     >  1:157242/1‑251 (MQ=255)                                                          tGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCt                                                                                                                                                                                  >  1:613125/1‑251 (MQ=255)                                                           gTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGGACCCAGTGGGCAATCTTTAGACCCAACCCGTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTAACTTCCTGGCCGCTCCACGCtt                                                                                                                                                                                 >  2:555736/1‑251 (MQ=255)                                                                            ggAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGGCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCTCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGtt                                                                                                                                                                 <  1:287680/250‑1 (MQ=255)                                                                            ggAAACCATCCGCAACAATGACGACGCCGTTATCAATTTCCCAGGTAGTCGGGTTCCCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGATCCACGCTTGCAATGAGTAATCGtt                                                                                                                                                                 >  2:425921/1‑250 (MQ=255)                                                                              aaaCCAGCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCTCTCCACGCTTGCAATGAGTAATCGTTAAc                                                                                                                                                              <  1:529100/251‑1 (MQ=255)                                                                                        caacaaTGACGCCACCGTTATAGATTTAACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGACTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACCTTCCCATTGAAAATAACATGG‑TAACTTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAACGAGTAATCGTTAACGCTA‑TCCGCt                                                                                                                                                   >  2:263697/1‑251 (MQ=255)                                                                                         aacaaTGACGCCGCCGTTCTCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTg                                                                                                                                                  >  1:536484/1‑251 (MQ=255)                                                                                                    cgccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGc                                                                                                                                      <  2:500123/251‑1 (MQ=255)                                                                                                         ttATCACTTTCACTGGCAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGAAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCCGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCt                                                                                                                                 >  2:309259/1‑251 (MQ=255)                                                                                                              aaTTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTAt                                                                                                                            <  1:463924/251‑1 (MQ=255)                                                                                                              aaTTTCACAGGTAGTCGGGTTAACGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTAt                                                                                                                            >  1:182481/1‑251 (MQ=255)                                                                                                                tttCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGAAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGAGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGCCGTTCTTGTAAACATTACCATTGAAATTAAACTGG‑TACATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGATA‑TCCGCTGGTGTTG‑ATGGCCCGCTGTTATcc                                                                                                                          >  2:346489/1‑251 (MQ=255)                                                                                                                     caGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAACTAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCggcgt                                                                                                                     >  2:392574/1‑248 (MQ=255)                                                                                                                          aGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTggg                                                                                                                 <  1:364078/250‑1 (MQ=255)                                                                                                                                      gTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGttt                                                                                                    >  1:29055/1‑251 (MQ=255)                                                                                                                                      gTACTGACTAAATACAGAAGCGGTAGCTGTAAGCTAGAGACGCCATGGGTTCTTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCCACCCATCAGGCGTTATTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGttt                                                                                                    >  2:650311/1‑251 (MQ=255)                                                                                                                                       tACTGACTCAATCCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACACATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTCTCCGGAGTTGGGGTCACCGGTGTTTc                                                                                                   >  2:644712/1‑251 (MQ=255)                                                                                                                                        aCTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTcc                                                                                                  >  1:343860/1‑251 (MQ=255)                                                                                                                                           gACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCgg                                                                                                >  1:380147/1‑250 (MQ=255)                                                                                                                                             cTCAATTCCGCAGCGGTAGCTGTACGCTAGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGt                                                                                             >  1:544672/1‑251 (MQ=255)                                                                                                                                                 atgccGCAGCGGTAGCTGTACGCTGGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTTCCCCGTGGGCAATCGTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGTAAAGAACATGG‑GAAATTTCGCGACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTGAACTCTA‑TCCGCTGGTTTGG‑GTGGCTCGCTTTTATCCGGAGTTGGGGTCACCGGTTTTTCCGGAGTTGGa                                                                                         <  1:455920/248‑1 (MQ=255)                                                                                                                                                 aTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGa                                                                                         <  2:314912/251‑1 (MQ=255)                                                                                                                                                   gccGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATc                                                                                       <  2:468678/250‑1 (MQ=255)                                                                                                                                                      gCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGCGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATccac                                                                                    >  2:288021/1‑248 (MQ=255)                                                                                                                                                      gCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAAc                                                                                    <  2:484920/251‑1 (MQ=255)                                                                                                                                                       cAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTGATCCTGAGTTGGGGTCACCGGTGTTTCCGGAGTTGTATCAAcc                                                                                   <  1:111289/251‑1 (MQ=255)                                                                                                                                                         gCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTTGACCCAGCCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATATCGCTCCCTTCCTGGCCGTTCCACGCTTGCAATGAGTAATCGTTAACTTTT‑TCAGCTGGTGTTGTATG‑CTATCTGTTATCTGGAGTTGAGGTCACCAGTGTTTCCAGCGATGGATCAAcagg                                                                                 >  2:368684/1‑248 (MQ=255)                                                                                                                                                           ggTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCggcg                                                                               >  2:473891/1‑249 (MQ=255)                                                                                                                                                                  gTACGCTCGAGACGCCATGGGATATTGGAGTTTTCATCGCTTGTACCCCGTGGGCAATCTTTAGACCCAGCCCATCAGGCTTTCGTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATAGCTCGCTGTTATCCGGAGTTGGGGTCAACGGCGTTTCCGGAGTTCGATCAATCGGAGTCACAg                                                                         <  1:381276/250‑1 (MQ=255)                                                                                                                                                                    aCGCTCGAGACGCCAAAGGATATTGGAGTTTTCAGCGCTGGCAACCCGTGGGCCATCTTTAGACCCAACCCATAAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTCAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGCGTTTCCGGAGTTGGATCAACCGGAGTCACAGCAc                                                                      >  2:411487/1‑251 (MQ=255)                                                                                                                                                                                     ggTTATTGGAGTTTTCCGCGCTGGTACACCGCGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAACTAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCACTGGGTAATCGCTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGCTCAACCGGAGTCACAGCACCGCCATTATTATCCGGc                                                     >  2:242933/1‑251 (MQ=255)                                                                                                                                                                                         aatGGAGTTTTCAGCGCTGGTACCCCGGGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTATAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCa                                                                       <  2:232993/227‑1 (MQ=255)                                                                                                                                                                                         aTTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTATAAACATTACCATTGAAAATCACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCa                                                                       >  1:232993/1‑229 (MQ=255)                                                                                                                                                                                          ttGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAAc                                                >  2:452843/1‑251 (MQ=255)                                                                                                                                                                                                 tttcgcccttGGTACCCCGTGGACAATTTTTAGCCCCAACCCCTCGGGCTTTCTTATAAACATTCCCATT‑AAAAAAAAATGGGTAATTTTCCCCCCTTTCCGGGCGCCCCCCCTCTTGAAAGGAGTAATCTTTAAC‑CTAATCCGCGGGTGTGG‑AGGGTTCGCTGTTATCCGGAGTTGGGGTCACCGTTGTTCCCGGAGTGGGATCACCCGGATTCACGGAACCGCCATTATTATCCGGCTTAACGGAACa                                          <  1:571532/241‑1 (MQ=255)                                                                                                                                                                                                        gcTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGcc                                    >  2:180215/1‑249 (MQ=255)                                                                                                                                                                                                               ccccGGGGTCAATCTTAAGCCCCACCCCATCGGGCGTTCTTGTAAACATTACCATAAAAAAAAACAGGT‑TAAATTTCGCTCCCTCCCTGGCCCCTCCACGCTTGAAATGAGTAATCGTTAACGCTA‑TCCGCGGGTGTTG‑AGGGCCCGCTGTTACCCGGAGTGGGGTTCCCGGGTGTTTCGGTAGTTGGATCACCCGGAGCCCCAGCCCCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAcc                           <  1:234177/251‑1 (MQ=255)                                                                                                                                                                                                                   gTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATAc                       >  2:529903/1‑251 (MQ=255)                                                                                                                                                                                                                         aaTCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTACTCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTg                 >  2:5272/1‑251 (MQ=255)                                                                                                                                                                                                                          aTCTTTAGACCCAAACCATCAGACGTTCTTGTAAACATTACCATTGAAAATCACATGG‑TACATCTCGCTACATTCCTGGCAGCTCCACGCTTGCAATGAGTAATAGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTg                 >  1:90594/1‑250 (MQ=255)                                                                                                                                                                                                                              ttAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCCTTGAAAATAACCTGG‑TAAATTAAGCTAACTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGATAACGCTA‑TCCGCTGGTGTTG‑CTGGATCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACAGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTc            >  1:643377/1‑251 (MQ=255)                                                                                                                                                                                                                                aGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATACCCTGG‑TAACTTTAGCTACCTTACTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCAGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGTTCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGAttttcg          >  2:67697/1‑248 (MQ=255)                                                                                                                                                                                                                                aGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTg          >  1:641627/1‑251 (MQ=255)                                                                                                                                                                                                                                   cccaacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttg       >  1:6008/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccaacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgc      >  2:2793/1‑251 (MQ=255)                                                                                                                                                                                                                                     caacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgcc     >  1:585190/1‑251 (MQ=255)                                                                                                                                                                                                                                      aacccaTCAGGCTTTCTTTTAAACATTCCCATTGAAAATAACTAGG‑TAATTTGCGCTACCTTCCTGGCCTCTCCCCGCTTGCAATGCTTAATCTTTAACCCTT‑CCCGCGGGTGTTG‑GTGGCTCGCTGTTATCCGGAGTTGGGTTCCCCGGTTTTCCCGGATTTGGATCACCCGGTGTAACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccg    <  1:507860/251‑1 (MQ=255)                                                                                                                                                                                                                                      aacccaTCAGGCTTTCGTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGATAACACTA‑TCCGCTGGTGTTG‑ATGGCTCGCGGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGGGTTGGATCAACCGGAGTCACAGCACCGCCATTATTCTCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGCTTTCTgttgccg    >  2:415141/1‑251 (MQ=255)                                                                                                                                                                                                                                        cccaGCAGGCGTTCTTGTAAACATTACCATTAAAAATAACATGG‑AAAATTTCGTTCCTTTCCGGGCCGCTCCACGCTTGAAATGAGTAACCTTTAACGCTA‑TCCGCTGGGGTTG‑ATGGCTCTCTGTTATCCGGAGTTGGGTTCACCGGTGTTTCCGGAGTTGGATCACCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccgtt  <  1:497671/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          CTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGG‑TAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTA‑TCCGCTGGTGTTG‑ATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTT  >  NC_000913/3468533‑3469014

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