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%	259.6 / NA	72	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (44/26);  minor base G (0/2);  total (44/28)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTT  >  NC_000913/3468526‑3469014                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         cATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCATGCGtt                                                                                                                                                                                                                                                  <  1:170959/251‑1 (MQ=255)  aTGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTc                                                                                                                                                                                                                                                 <  1:110201/251‑1 (MQ=255)    gCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCtt                                                                                                                                                                                                                                               <  1:489163/251‑1 (MQ=255)           gATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAAGTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACa                                                                                                                                                                                                                                        <  1:549225/251‑1 (MQ=255)            aTCGACAAAAGTTTTATATTGGGCATCGTTATACTCTACGACGCTGCCCGCGCTGTCCGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCGTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACAt                                                                                                                                                                                                                                       <  2:57249/251‑1 (MQ=255)               gACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGGTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTAc                                                                                                                                                                                                                                    >  1:414723/1‑251 (MQ=255)                aCAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCATTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTAcc                                                                                                                                                                                                                                   >  2:226779/1‑251 (MQ=255)                  aGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACACTTTGCTGGCCTGGAAACCATCCGCAACAATGACCCCGACGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCCTGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCACCAC‑ATCAGGCGTTCTGGTAACCATTACCAt                                                                                                                                                                                                                                 >  2:519045/1‑251 (MQ=255)                    aaGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGAGCTGTACGCTTTGCTGGCCTGGAAACCATCCGAAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTg                                                                                                                                                                                                                               >  1:258824/1‑251 (MQ=255)                      gTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGa                                                                                                                                                                                                                              >  1:20755/1‑250 (MQ=255)                             aaTGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCACTCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAAc                                                                                                                                                                                                                       >  1:157397/1‑250 (MQ=255)                             aaTGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAAc                                                                                                                                                                                                                       >  1:170086/1‑250 (MQ=255)                              aTGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCCACAATGACGCCGCCGTTATCAATTTCACCGGCAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACAt                                                                                                                                                                                                                     >  1:546144/1‑251 (MQ=255)                                gtgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATg                                                                                                                                                                                                                    <  2:191433/250‑1 (MQ=255)                                   cATCGTTATAACCTACGTTGCTGTCCGCGCTGTCCGCTCTGCTGGCCGGGACACCATCCGAAACAAGGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTGCCATTGAAAATAACATGGTaa                                                                                                                                                                                                                <  2:442358/251‑1 (MQ=255)                                          ataATCTACGATGCTGTCCGCGCTTTACGCTTTGCTGGCCTGTAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCg                                                                                                                                                                                                          <  1:392652/250‑1 (MQ=255)                                          ataATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGc                                                                                                                                                                                                         >  1:500348/1‑251 (MQ=255)                                            aaTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACCATGACGCCGCCGTTATCAATTTCACCGGTCGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTCGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCCGCGCTGGTACCCCGTGGGCAATCTTTAGACCAAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCt                                                                                                                                                                                                        >  2:354556/1‑250 (MQ=255)                                              actccGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGGCATTGGTTATTTGAGTTTTCAGCGCTGAAGCCCCGTTTGCAATCTTTAGGCCCAACCC‑AGCAGGCGTTCTTGGAAACATTACCATTGAAAATAACATGGTAAATTTCGCTAcc                                                                                                                                                                                                     <  1:38648/247‑1 (MQ=255)                                                 ccGATGCGGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTc                                                                                                                                                                                                  <  2:170086/250‑1 (MQ=255)                                                   gATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTAGCTACCTTCCt                                                                                                                                                                                                >  1:197158/1‑251 (MQ=255)                                                      gCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGc                                                                                                                                                                                             >  1:396448/1‑251 (MQ=255)                                                         gTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGc                                                                                                                                                                                          >  1:92561/1‑251 (MQ=255)                                                         gTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGc                                                                                                                                                                                          >  2:458862/1‑251 (MQ=255)                                                                         ttGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAAt                                                                                                                                                                            <  2:157397/249‑1 (MQ=255)                                                                               gCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGCAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGCGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTCAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTactc                                                                                                                                                                    >  2:457235/1‑248 (MQ=255)                                                                                   ggAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGtt                                                                                                                                                                 <  1:458862/250‑1 (MQ=255)                                                                                        ccATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGGACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCGTTAGACCCAACCC‑ATCAGGCGGTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCt                                                                                                                                                           <  1:574483/251‑1 (MQ=255)                                                                                         cATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTa                                                                                                                                                          <  1:124646/251‑1 (MQ=255)                                                                                         cATCCGCAACAAGGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACAGTACCATTGAAAATAACATGTTAAATGTCGCTACCTTCCTGGCCGCTGCACGCTTGCAGTGAGTAATCGTGAACGCTa                                                                                                                                                          <  2:546144/251‑1 (MQ=255)                                                                                                  caaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTggg                                                                                                                                                 >  2:455899/1‑250 (MQ=255)                                                                                                   aaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGgtg                                                                                                                                                >  1:421646/1‑251 (MQ=255)                                                                                                        cgccgcCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGAt                                                                                                                                           >  1:532754/1‑251 (MQ=255)                                                                                                          ccgccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAGTTCGCTACCTTCCTGTCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATgg                                                                                                                                         <  1:354556/251‑1 (MQ=255)                                                                                                           cgccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGc                                                                                                                                        >  1:320897/1‑251 (MQ=255)                                                                                                           cgacgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGc                                                                                                                                        <  1:108811/248‑1 (MQ=255)                                                                                                                  aTCAATTTCCCAGGTAGTCGGGTTACCGTACTGACTCAATTCCGAAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCCCGCTGt                                                                                                                                 >  2:361097/1‑251 (MQ=255)                                                                                                                    cAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTa                                                                                                                               <  2:104956/251‑1 (MQ=255)                                                                                                                         gcaaaGGTAGTCGGGTTGCCGTACTGACTCAATTCCGCAGCGGTAGCTGTACTCGCTAGACGCCATGGGTTCTTGGAGTATTCAGCGCTGGTCCCCACTGGGCAATCTGTAGACCCAGCCC‑ATCAGGCGTGCTTGTTAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCGTCCTGGCCGCTCCACGCTTCCTAGGATTAATCTTTGACGCTACCCGCGGGTGTTGATGGTTCGCTATTAACCgg                                                                                                                          <  2:540648/247‑1 (MQ=255)                                                                                                                           acaGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAg                                                                                                                        <  1:481792/251‑1 (MQ=255)                                                                                                                                  gTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGCGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTggtgt                                                                                                                 >  1:130544/1‑248 (MQ=255)                                                                                                                                     gggTTACCGTCCTGAATCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACACAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAc                                                                                                              >  1:470076/1‑251 (MQ=255)                                                                                                                                      ggTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAcc                                                                                                             >  1:19956/1‑251 (MQ=255)                                                                                                                                               aCTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCACTCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTcc                                                                                                  >  1:59601/1‑251 (MQ=255)                                                                                                                                               aCTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTGGTAAACATTACCATTTAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTGTCCTCTGGGGTTGAGGGCTCGCTGTTATCCGGAGTTTGGGTCACCGG‑G‑GTTTcc                                                                                                  <  2:257817/251‑1 (MQ=255)                                                                                                                                                  gACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCGTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTGCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGATTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTGATCCGGAGTTGGGTTCACCGG‑T‑GGTTCCg                                                                                                 <  2:254003/249‑1 (MQ=255)                                                                                                                                                   aCTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAAACC‑ATCAGGCGTTCTTGTAAATATTACCATTGACAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCgg                                                                                                >  1:307288/1‑249 (MQ=255)                                                                                                                                                       aaTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGAGGGCTCGCAGTTATCTGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTgg                                                                                          >  2:82436/1‑251 (MQ=255)                                                                                                                                                        aTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTgg                                                                                          >  1:385140/1‑250 (MQ=255)                                                                                                                                                           ccGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATc                                                                                       <  1:294614/250‑1 (MQ=255)                                                                                                                                                             gcggCGGTAGCTGTACGCTCGAGCCGCCATGGGTGATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCGTGAGACCCAACCC‑AGCAGGCGTTCTTGTAAACAATTCCATTGAAAATAACATGGTAAATTTCGCTACCTTGCTGGCCGCTCCACGCGTGCAATGAGTAATCGGTAGCGCTATCCGCTGGTGTGTATGGCGCGCCGTTAGCCGGTGTAGGGGTCACGGC‑T‑TTTTCCGGGGTTAGATCAAc                                                                                    <  2:303155/248‑1 (MQ=255)                                                                                                                                                             gCAGCGGAAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAAc                                                                                    >  1:514674/1‑251 (MQ=255)                                                                                                                                                              cagtggTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAAcc                                                                                   <  1:468783/247‑1 (MQ=255)                                                                                                                                                              cAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGGGTTGCTGGCTCGCTGTTATCCGTAGTTGGGGTCACCGG‑T‑GTTTCCGGCGTTGGATCAAcc                                                                                   >  2:532971/1‑251 (MQ=255)                                                                                                                                                                   gTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCgg                                                                                                                          >  1:1164/1‑209 (MQ=255)                                                                                                                                                                   gTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCgg                                                                                                                          <  2:1164/209‑1 (MQ=255)                                                                                                                                                                       cTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCCCCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTcac                                                                           >  2:482003/1‑250 (MQ=255)                                                                                                                                                                                  agaCGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCatt                                                               >  1:222383/1‑251 (MQ=255)                                                                                                                                                                                              ttATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTCGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCtt                                                   >  1:542601/1‑251 (MQ=255)                                                                                                                                                                                                     aGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGa                                             >  1:268707/1‑250 (MQ=255)                                                                                                                                                                                                      gTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAAc                                           >  1:113466/1‑251 (MQ=255)                                                                                                                                                                                                          tCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAgg                                        >  1:164552/1‑250 (MQ=255)                                                                                                                                                                                                                         ccgtggacaatcgtaaacccaaccccATCAGGCTTTCTTGAAAACATTACCATTAAAAATAACATGTAAAATTTCTCCCCCTCCCGGCCCCCCCCCCGCTTACAATGAGTAATGGTTAACGCTGTCGGCTGGGGTGGATGGCCCGCTGTTA‑CCGGGGGTGGGGGCCCCGGGT‑GTTTCCGGAGTTGGACAAACCGGAGCCCGAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCAt                         <  2:268707/235‑1 (MQ=255)                                                                                                                                                                                                                          gTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATAc                       >  2:461659/1‑251 (MQ=255)                                                                                                                                                                                                                                         acccaaccc‑aTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgtt        >  1:329002/1‑251 (MQ=255)                                                                                                                                                                                                                                          cccaacca‑aTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttg       >  1:452996/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccacccc‑aTCAGGGGTTCTTGTAAACATTCCCATTAAAAATAACATGGTAAATTCCGCTACCTCCCGGGCCGCCCCACGCTGGCAAGGATTAAGCTTAAACCCTATCCGCTGGTGTTGATGGTTCGCTTTTATCCGGAGTGGGGGTC‑CCGG‑TTGTTTCCGGAGTTGAATCACCCGGAGCCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgc      <  1:48361/247‑1 (MQ=255)                                                                                                                                                                                                                                            caaccc‑aTCAGGCGTTCTTGTAAACATTACCATTGAAAATCACATGGTAAATTTCGCTACCTTACTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgc      >  2:48361/1‑250 (MQ=255)                                                                                                                                                                                                                                            caaccc‑aTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGCGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgtc     >  2:588624/1‑249 (MQ=255)                                                                                                                                                                                                                                             aaccc‑aTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccg    >  1:224372/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaccc‑aTCAGGCGCTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccg    >  1:135209/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccc‑aTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTTCCATACTGGCTGATTTCTgttgccgtt  >  1:48019/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         CATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGG‑T‑GTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTT  >  NC_000913/3468526‑3469014

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