Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	3,951,159	T→C	intergenic (‑123/+72)	PP_3486 ← / ← PP_3487	putative cytochrome c/conserved exported protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	3,951,159	0	T	C	100.0%	153.2 / NA	46	intergenic (‑123/+72)	PP_3486/PP_3487	putative cytochrome c/conserved exported protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/23);  total (23/23)

AACTGGCGCAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCC  >  AE015451/3950874‑3951445                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                 aaCTGGCGCAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAgg                                                                                                                                                                                                                                                                                   <  4:271882/301‑1 (MQ=255)                         gggggACAGGGATGTCGCGGGCCCCGACCCGGGCTGGCAATTTCGCGTCGCTACTGGAGACGCCCGCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCca                                                                                                                                                                                                                                                          <  2:463469/300‑1 (MQ=255)                             gACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCCCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACa‑ga                                                                                                                                                                                                                                                     >  4:514960/1‑301 (MQ=255)                                 cctttccctccccgaccctcttccgctctcAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGcc                                                                                                                                                                                                                                                 <  1:633005/272‑1 (MQ=255)                                        cGCTGGCGCCGACTCGGTCTTGCAAGTGAGCGTCGCAACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCGGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCAGGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTccc                                                                                                                                                                                                                                          <  1:165567/301‑1 (MQ=255)                                          cTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCCCCGTTCGCGACGCAAGCGCGC‑CCCACCAAAAGCGGATGCCACA‑GCGGCCaagtacccc                                                                                                                                                                                                                                        >  1:590743/1‑294 (MQ=255)                                                   aCCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCCCCCA‑CAAAATCCGA‑GCCCCACGCGGCCAATTTCCCCGGCGGcccc                                                                                                                                                                                                                               >  2:450002/1‑300 (MQ=255)                                                          cTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACC‑GAGGCCAATTCCCCAGGCCGCCCATGATAcc                                                                                                                                                                                                                        >  2:477143/1‑300 (MQ=255)                                                              cAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCagatcggaagagcacacgtctgaactcca                                                                                                                                                                                                                    >  3:633019/1‑273 (MQ=255)                                                                    aGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACCCACCGACCgg                                                                                                                                                                                                              >  1:498074/1‑301 (MQ=255)                                                                    aGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGCCGGCCCATGATACTCACCGACCgg                                                                                                                                                                                                              >  4:245091/1‑301 (MQ=255)                                                                      cGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGcc                                                                                                                                                                                                            <  3:498089/301‑1 (MQ=255)                                                                       gTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCCAAAGGGCTGCGCCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGAAAGGCCGCTCCCA‑CAAAAGCGGCTGCCACC‑AAGGCCAATTCGCCAGGCCGCCCATGATACTCACCGACCGGccc                                                                                                                                                                                                           >  1:161822/1‑300 (MQ=255)                                                                                                tAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCCGGGTCAACAGCTGCATgcg                                                                                                                                                                                  >  3:544687/1‑301 (MQ=255)                                                                                                       gTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCCTGATACTCACCGACCGGCCAGCGGCTGGGTCACCAGCTGCagtcgcggggca                                                                                                                                                                           >  4:622153/1‑290 (MQ=255)                                                                                                                   gTTCGTGGTCGGTCAGGTCCAGGGCAAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGc                                                                                                                                                               <  3:378815/301‑1 (MQ=255)                                                                                                                       gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg                                                                                                                                                           <  3:548596/301‑1 (MQ=255)                                                                                                                                cAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTTGGCAATCGCCGCTGGCGGTGCATCCgggg                                                                                                                                                  >  1:482437/1‑300 (MQ=255)                                                                                                                                     ccagggccagggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATc                                                                                                                                             <  4:389306/301‑1 (MQ=255)                                                                                                                                           ccagggtggagttCGAAGTTGCCATGCCAAGCCGGCAGGTAAGCCATATCGCAAGCCCTGTGGGGGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTc                                                                                                                                       <  1:544673/289‑1 (MQ=255)                                                                                                                                            aagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTcc                                                                                                                                      <  3:482451/300‑1 (MQ=255)                                                                                                                                                    gggTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTa                                                                                                                              <  4:450111/300‑1 (MQ=255)                                                                                                                                                     ggTCGAAGTCGCCATCCCAAGACCGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTAc                                                                                                                             <  2:514834/301‑1 (MQ=255)                                                                                                                                                           aGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCaa                                                                                                                       >  1:47363/1‑301 (MQ=255)                                                                                                                                                            gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCCGCCGGGTCAACAGCTGCATGCGCTGGGCAATCGTCGCTGGCTGTGCAGCCGGGCGAATCTTGCTCCAGCATTTACCTGCAAc                                                                                                                      >  2:86312/1‑301 (MQ=255)                                                                                                                                                                  atgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGACAGCGGCGGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACAGCAcc                                                                                                                >  4:177476/1‑301 (MQ=255)                                                                                                                                                                   cgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCg                                                                                                               <  3:47364/300‑1 (MQ=255)                                                                                                                                                                          ccggagggtaggcagtgcgcagaGCCCTGTGGGGGCGGCATGTCGGCGGGAAAGGGCTGCGCAGCAGCCCCGGAGGGGGCAGCGCAGGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTAGCGAGGCAAGGCCGCTCCCA‑AAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCa                                                                                                        <  2:284164/279‑1 (MQ=255)                                                                                                                                                                          cACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCCTCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATcc                                                                                                        >  1:410124/1‑300 (MQ=255)                                                                                                                                                                                          atCGAAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCGGTAGCAGCTCAGGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAAtgctg                                                                                        <  4:86331/301‑1 (MQ=255)                                                                                                                                                                                               aaGCCCTGTGGGAGCGGCCTGGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCgg                                                                                   <  3:161828/301‑1 (MQ=255)                                                                                                                                                                                                                 cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGcc                                                                 >  4:248002/1‑301 (MQ=255)                                                                                                                                                                                                                   ggcggcgcgaaggggcggctgggcagcCCCGGCGGTGGCAGCGAAAGCGCAGACAGAGGGGAGCGCTGGGCGCTCCTTTCGCGCCGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCt                                                               <  2:247958/279‑1 (MQ=255)                                                                                                                                                                                                                      gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCtata                                                            >  1:536465/1‑301 (MQ=255)                                                                                                                                                                                                                       gcgggaaagggtggagcggcggcgccggaggTGGCAGCACAAGCGCGGAAAGAGGGGAGCGCTGGGCGCTCCTTTCGGGAGGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATAc                                                           <  2:245047/271‑1 (MQ=255)                                                                                                                                                                                                                       cacgacgctcttccgatctcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAAATGCCCTATAc                                                           <  1:601732/282‑1 (MQ=255)                                                                                                                                                                                                                              aGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAg                                                    <  4:295214/301‑1 (MQ=255)                                                                                                                                                                                                                                  cTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCAcc                                                >  4:469452/1‑301 (MQ=255)                                                                                                                                                                                                                                     cgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGTCAGGAACTGCCCTATACGGCCAAGCACCAgg                                             >  4:98994/1‑301 (MQ=255)                                                                                                                                                                                                                                      gcagcaGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGt                                            >  3:111632/1‑301 (MQ=255)                                                                                                                                                                                                                                          cagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAAATGCCCTATAcagatcggaagagcacacgt                                        >  3:601746/1‑282 (MQ=255)                                                                                                                                                                                                                                               ccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTgg                                   >  3:364950/1‑301 (MQ=255)                                                                                                                                                                                                                                                              cAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCa                    <  4:477256/301‑1 (MQ=255)                                                                                                                                                                                                                                                                aGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAgg                  >  1:589759/1‑301 (MQ=255)                                                                                                                                                                                                                                                                      aaaaGAGGGGCGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGc            <  3:410134/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                ggcgcTGCGCGCTCCTTTCGCGAGGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGGTACTCACCGACCAGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGcc  <  3:536479/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                 AACTGGCGCAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCA‑CAAAAGCGGATGCCACA‑GAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCC  >  AE015451/3950874‑3951445

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑