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	95.8%	205.2 / ‑3.4	71	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (41/27);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

CCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCC  >  NC_000913/3468524‑3469011                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cccATGCGTCTTGATCGACAGACGTTTTATAATACGCATAGTTATAATCTACGCTGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAACTGTACGCCCGCGACGCCATGGGTTATTGGAGTTTTCAGTCCTGTTAACCCGTGGGCACTTTTTAGAACCAAAccacacggcg                                                                                                                                                                                                                                                 >  2:487743/1‑245 (MQ=255) cccATGCGTCTTGATCGACAGAATTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACTCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCg                                                                                                                                                                                                                                                 <  2:671593/251‑1 (MQ=255) cccATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCg                                                                                                                                                                                                                                                 <  1:557113/251‑1 (MQ=255)   caggcttCTTGATCGCCAGAAGTTTTATAATGTGCATCGTTATATTCTACGATGTTGTCCGCGCTGTACGCTTGGCTTGCTGGGAAACCATCCGCAACAATGTCGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATGCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGtt                                                                                                                                                                                                                                               <  2:736004/245‑1 (MQ=255)      gCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCtt                                                                                                                                                                                                                                            <  1:417573/251‑1 (MQ=255)             gATCGACAGAAGTTTTATAATGTGCATCGTTATAACCTACGATGCTGTCCGCGCTTTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAAc                                                                                                                                                                                                                                      <  1:817388/250‑1 (MQ=255)                   cAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCa                                                                                                                                                                                                                               <  1:568958/251‑1 (MQ=255)                    aGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTTCTGGCCTGGAAACCGTCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCAt                                                                                                                                                                                                                              <  1:757742/251‑1 (MQ=255)                    aGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCgagatcggaagagcgtctgg                                                                                                                                                                                                                              >  2:33379/1‑232 (MQ=255)                    aGAAGTTTTATAAGGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCg                                                                                                                                                                                                                                                 <  1:33379/232‑1 (MQ=255)                         ttttATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGaa                                                                                                                                                                                                                          >  1:709750/1‑250 (MQ=255)                           ttATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAAt                                                                                                                                                                                                                       >  1:41792/1‑251 (MQ=255)                               aaTGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAAc                                                                                                                                                                                                                    >  1:634840/1‑250 (MQ=255)                                 tgtgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATg                                                                                                                                                                                                                 >  1:791233/1‑251 (MQ=255)                                    gCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGATCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAAAATGGt                                                                                                                                                                                                               >  2:653893/1‑250 (MQ=255)                                      ttctttATCATCTACGATGCCGTCCGTGCCTTTCGCTGCGCTGGCCTTAAAACCCTCCTCAACATTATCTCCGCCCTTATCAATTTCACATGTAGTCGGGTGTCCGTACTGGCTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGTTTATTGTTGTTTTCAGCGCTGGGCCCCCGTGGGCAATCTTTATACCCAACCCATCAAGCGTTCTTGTAAACATTACCATCGAAAAGAACATGGTa‑a                                                                                                                                                                                                            <  2:484372/246‑1 (MQ=255)                                         gTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AAttt                                                                                                                                                                                                        <  1:612871/251‑1 (MQ=255)                                          tgataATCTACGATGCTGAACGCGTTGTACGCTTTGCTGGCCTTGAAACCATCCGAAACACTGCCGCCGCCGTTATCCATTTATCAGGTAGTCGGGTTACCGTACTGACTCAAATCCGCAGCGGTAGCAGTACGCTCGAGACGACATGGGTTATTGGAGTTTTCAGCGCTGGTACACCGTGGGCCATCTTTAACCCCAACCCATCAGGCGATCTTGTAAACATTACCATTGAACATACCATGGTA‑ATTTTc                                                                                                                                                                                                       >  2:526026/3‑251 (MQ=255)                                            ataATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTCTTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCACACCCATCAGGCGTTCTTGTACACATTACCATTGAAAATAACCTGGTA‑AATTTcac                                                                                                                                                                                                     >  2:777982/1‑249 (MQ=255)                                              aaTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGATCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCt                                                                                                                                                                                                    <  1:653893/250‑1 (MQ=255)                                               aTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGACACCATCCGCAACAATGACGCCGCCGTTATCAATTTCAAAGGTAGTCGGGTTACCGTACTGACTCCATTCCGCAGCGGTAGCTGTACGCTCGCGACGCCATGGGCTATTGGAGTTTTCAGCGCTGGCACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGACACTAACATGGTA‑ATTTTCGCttc                                                                                                                                                                                                  >  2:378078/1‑249 (MQ=255)                                                 cTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCATTCTTTAGACCAAACCGACCAGGCGTTCTTGTAAAATTTACCATTGAAAATAAAATGGGA‑AATTGCGCTACCt                                                                                                                                                                                                >  2:467495/1‑251 (MQ=255)                                                 cTACGATGCTGTCAGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTTAAAATCACATG‑TATAATTTCGCTACCt                                                                                                                                                                                                >  2:821061/1‑251 (MQ=255)                                                           gTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGc                                                                                                                                                                                      >  1:148024/1‑251 (MQ=255)                                                           gTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCACCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGcc                                                                                                                                                                                        <  1:729211/249‑1 (MQ=255)                                                              cgcgCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTcc                                                                                                                                                                                   >  1:287553/1‑251 (MQ=255)                                                                      aCGCTTTGCTGGCATGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGc                                                                                                                                                                           <  1:102004/251‑1 (MQ=255)                                                                            tGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTCGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGTGTTTTCAGCGCTGGTAGCCCGTTGGCAATATTTAGACCCAACCCATCAGGCGTTCTTGTAAACCTTACCATTGAAAATAACATGGTA‑AATTTCGCTTCCTTCCTGGCCGCTCCACGCTTGCaagggg                                                                                                                                                                     >  2:350776/1‑247 (MQ=255)                                                                                ggCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAAt                                                                                                                                                                 <  1:518465/251‑1 (MQ=255)                                                                                gcccTGCAAACCACCCTCACCAATGACGCCGCCGTTAGCAAGTTCACAGGGAGTCGGGTGACCGTACTGACTCTATTCCGCAGCGGTAGCTGTACTCTCGAGACGCCATGGGTTATTTGAGTTGTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACTCATCAGGTGTTCTTGTAAACATTACCATTGATAATATCATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTTCAATTAGTaa                                                                                                                                                                  <  2:538204/248‑1 (MQ=255)                                                                                 gCCTAGAAACCACCCGAAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATc                                                                                                                                                                <  2:42536/251‑1 (MQ=255)                                                                                    tGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGtt                                                                                                                                                             >  1:669884/1‑251 (MQ=255)                                                                                      gAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTATTCGTTac                                                                                                                                                           >  1:862809/1‑250 (MQ=255)                                                                                                  aacaaTGGCCCCGCCGTTAACAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTTGCCGTTCGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTGTAGACCCAATCCATCAGGCGTTCTTGAAAACATTACCATTGAAAATAACATGGTA‑AATGTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATAGTTAACGCTATCCGCTg                                                                                                                                               <  2:243/251‑1 (MQ=255)                                                                                                  aacaaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGAGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCAGTGGGCAATCTTTAGCCCCAACCCATCAGGCGTTCTTGTAACCATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATAGTTAACGCTATCGGCTg                                                                                                                                               >  2:742228/1‑251 (MQ=255)                                                                                                        gccgccgcCATTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCATCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTg                                                                                                                                         <  2:41792/249‑1 (MQ=255)                                                                                                           gccgccGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGCGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTATTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGTTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGa                                                                                                                                        >  2:106284/1‑249 (MQ=255)                                                                                                              gccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCt                                                                                                                                   <  1:365124/251‑1 (MQ=255)                                                                                                               ccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTc                                                                                                                                  >  1:444368/1‑251 (MQ=255)                                                                                                                   tATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGCGACGCCATGGGTTATTGGAGTTTTCAGCGCTTGTACCCCGTGGGCAATCTTTAGCCCCAACCCATCAGGCGTTCTTGTAACCATTACCATTGCAACTAACCTTGTA‑AATTTCGCTCCCTTCCGTGCCGCTCCACGCTTGCAATGTGTACTCGTTAACGCTATCCGCTGGTGTTGATGGCTctctt                                                                                                                              >  2:347888/1‑247 (MQ=255)                                                                                                                      caagttCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGGACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTGCCTTCCTGGCCGCTCCACGCGTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTa                                                                                                                           <  1:84332/247‑1 (MQ=255)                                                                                                                          tgcacaGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGGACGCTCGAGACGTCAGGGGTGATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTTGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATcc                                                                                                                        <  2:634840/248‑1 (MQ=255)                                                                                                                                gttAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGTTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTg                                                                                                                 <  1:600748/249‑1 (MQ=255)                                                                                                                                     tCGGGTTACCTTACTGGCTCATTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTg                                                                                                                 <  2:719442/246‑1 (MQ=255)                                                                                                                                     tCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTg                                                                                                                 >  1:719442/1‑246 (MQ=255)                                                                                                                                      cGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTc                                                                                                            >  2:113357/1‑250 (MQ=255)                                                                                                                                        ggTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAcc                                                                                                         <  1:481253/251‑1 (MQ=255)                                                                                                                                            acagtaccgagtCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGgtg                                                                                                     <  1:79802/240‑1 (MQ=255)                                                                                                                                                     aCTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTGTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAg                                                                                            <  1:487249/251‑1 (MQ=255)                                                                                                                                                        cAATTCCGCAGCGGTAGCTGTACGCTCGAGACGACATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTCCCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGTCTCGCTGTTATCCGGCGTTGGGGTAACCGGTGTTTCGGTAGTTg                                                                                         >  2:419559/1‑251 (MQ=255)                                                                                                                                                             cctcAGCGGTAGGTGGACTCTCGAGATGCCATGGGTGATTGGAGTTTTCCGCGATGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATTGTA‑AATTTCGCTACCTTCCTTGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTCATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGAta                                                                                     <  2:73907/247‑2 (MQ=255)                                                                                                                                                              cGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGACCGCTCCACGCTTGCCATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGCGTTGGGGTCACCGGTGTTTCCGGAGGTGGCTCa                                                                                    >  2:536921/1‑250 (MQ=255)                                                                                                                                                                    ggTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCggcg                                                                             >  2:50944/1‑249 (MQ=255)                                                                                                                                                                     gTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGt                                                                            >  1:336368/1‑251 (MQ=255)                                                                                                                                                                      tAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCATGGCCGCTCCACGCGTGCAATGAGCAATCGTTAACGCAATCCTCTGATGTTAGTGGCTCGCTGTTATCCGGAGTAGGGGTCACCGGTGTTTCAGGAGTTGGATGAACCGGAGTc                                                                           >  2:779081/1‑251 (MQ=255)                                                                                                                                                                       aGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCAGCTGGTGTTGATGGCTCGCTGTTATCCGGTGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTcc                                                                          >  1:67529/1‑250 (MQ=255)                                                                                                                                                                                 tCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGTCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGcc                                                               >  1:751741/1‑251 (MQ=255)                                                                                                                                                                                   gagaCGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCat                                                             >  1:791113/1‑251 (MQ=255)                                                                                                                                                                                           aTGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATcc                                                     >  1:760136/1‑251 (MQ=255)                                                                                                                                                                                                                     cccccccGTGGGCAATCTATAGACCCACCCCACCAGTCTTTTTTGAAAAATTTACTATAAAAAATAACAGTGTA‑AATTCCGCTCCCTTCCTGCCCCCTCCCCGCTTGCAAGGAGTAATCTTTACCGCTTCCCGCTGTTGTTGATGGCTCGCTGTTATCCGGATTGGGGGCCCCCGGTTTTCCCGGAGTGGGATCAACGGGAGCC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGtt                           <  1:706643/248‑1 (MQ=255)                                                                                                                                                                                                                       aCCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAAc                                                                                  <  1:706114/195‑1 (MQ=255)                                                                                                                                                                                                                       aCCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGGGTTTCCGGAGTTGGATCAAc                                                                                  >  2:706114/1‑195 (MQ=255)                                                                                                                                                                                                                         cccGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAcc                        >  1:66958/1‑250 (MQ=255)                                                                                                                                                                                                                         cccGTGGGCAATCTTTAGACACAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGAAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGG‑GTCCACAGCACCGCCATTATTTTCCGGCTTACAGGAAAAGGAGCCAGGGta                           >  2:706643/1‑246 (MQ=255)                                                                                                                                                                                                                            gTGGGCAAACTTTAGAACCAACCCATCAGGCGTTCTTGTAAACATTCCCATTGAAAATAACATGGTA‑AATTTCGCTCCCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGGGTTGATGGCTCGCTGTTATCCGGAGTTGGTGCCACCGGAGTTTCCGTAGTTGGCTAAACCGGAGTC‑ACAGAACCGCCATTATTTTCCGGCTTAACGGAACAGGACCCAGGGATACCATAc                    >  2:576572/1‑251 (MQ=255)                                                                                                                                                                                                                                 cAATCATTAGCCCCAACCCATCAGGCGGTCTTGTCAACATTACCATTGGAATTAACATGATA‑AATTTCGCTACCTTCCTGGCCGCACCACGCTTGCAATGAGTAATCGTTAACGCTATCAGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCAGCGTTGGATCAACCGGTGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGATTACCATACTGGCt               >  2:800599/1‑251 (MQ=255)                                                                                                                                                                                                                                         aGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTAGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACTGAACAGGAGCCAGGGTGACCATACTGGCTGATTTCTg       >  1:206075/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccccaacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgtt     >  1:644770/2‑251 (MQ=255)                                                                                                                                                                                                                                           acccaacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgtt     >  1:644753/1‑251 (MQ=255)                                                                                                                                                                                                                                           acccaacacaTCAGGCGATCTTGTAAACATTACCATTGCAACTCACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGCGTAATCGTTAACGCTATCCGCTGGTGTTGTTGTCCCGCTGTTATCCGGAGTTGGGGTCACCGGTTTTTCCGGAGTTGGATCAACCGGAGTG‑ACAGCCCCGCCATTCTTATCCGCCTTACAGGAACAGGAGCCAGGGTTAACATACTGGCTTATTTCTgtt     >  2:639864/1‑251 (MQ=255)                                                                                                                                                                                                                                            cccaacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttg    >  2:147227/1‑251 (MQ=255)                                                                                                                                                                                                                                              caacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgcc  >  1:767006/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTA‑AATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTC‑ACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCC  >  NC_000913/3468524‑3469011

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