Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	336,125	+T	coding (210/237 nt)	PP_0278 →	conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	336,124	1	.	T	100.0%	125.9 / NA	37	D70V (GAT→GTT)	PP_0278	conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (14/23);  total (14/23)

TTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCAT  >  AE015451/335826‑336370                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       ttCACTCGGTTTCTCACCCTGGGTAAGGCTGCCGCCCCAAATCCGGTAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑tt                                                                                                                                                                                                                                                      <  2:263715/301‑2 (MQ=255)    ccTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTTAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑ttacc                                                                                                                                                                                                                                                   <  2:659054/300‑5 (MQ=255)                 cccTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGc                                                                                                                                                                                                                                       >  2:628437/1‑301 (MQ=255)                  ccTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGAGTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGcc                                                                                                                                                                                                                                      <  2:688487/301‑1 (MQ=255)                       ggaCGGCTGCCGCCCAAAAACCGTGAGAAAGAATTTTGGCAGTAGCCTGCCGTTAATAAGGAATTGCCATGAAAAAAGATGGTTCTGAGCCAATCACAACCGCAGGCGTGAAAACCTTCCTCGCAGCCGTTAAGCCATCCCTCAATCTCCTACGCTTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGGGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGgc                                                                                                                                                                                                                                 <  2:819010/299‑1 (MQ=255)                            gCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCTCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                                                                                                                            >  1:829407/1‑301 (MQ=255)                               cccgccCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtgatga                                                                                                                                                                                                                         <  1:253705/300‑1 (MQ=255)                                 cgccCCAAATCCGTTGGCAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAAc                                                                                                                                                                                                                       <  2:65689/301‑1 (MQ=255)                                    cccAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGaa                                                                                                                                                                                                                    >  2:480341/1‑300 (MQ=255)                                      cAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGatat                                                                                                                                                                                                                  <  1:245816/301‑1 (MQ=255)                                                         tttGGAAGTCTCCTTCCGAGAATAAGGATTTGCCATGAAAAAAGAGGTTTCAGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAg                                                                                                                                                                                               <  2:829117/301‑1 (MQ=255)                                                          ttGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAgg                                                                                                                                                                                              <  1:108651/301‑1 (MQ=255)                                                            ggCAGTCTCCTTCCGTGAATAAGGATTTGCAAGGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                                                                                                                                            <  2:9336/301‑1 (MQ=255)                                                                        ccGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTGCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTg                                                                                                                                                                                <  2:712684/301‑1 (MQ=255)                                                                         cGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGa                                                                                                                                                                               >  1:845293/1‑301 (MQ=255)                                                                                   gATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAtt                                                                                                                                                                     <  1:382634/301‑1 (MQ=255)                                                                                    aTTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                                                                                                                                    >  1:633760/1‑301 (MQ=255)                                                                                                  aaGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                                                                                                                                      <  2:844992/301‑1 (MQ=255)                                                                                                   aGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTg                                                                                                                                                     >  1:106162/1‑301 (MQ=255)                                                                                                     aTGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTTTGGGAGCGGCCTTGcc                                                                                                                                                   >  1:799536/1‑300 (MQ=255)                                                                                                       ttttCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTc                                                                                                                                                 <  2:47821/300‑1 (MQ=255)                                                                                                               gCCAATCCCAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATgg                                                                                                                                         <  1:440190/301‑1 (MQ=255)                                                                                                                   aacacaACCGCCGGGTTGGAAACCGTCCGCGCCGCCGTTAACCCATCCCCCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACAGCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                                                                                                                                     <  2:32649/299‑1 (MQ=255)                                                                                                                        cgatctaGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAag                                                                                                                                <  1:730803/295‑1 (MQ=255)                                                                                                                              cggCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGccc                                                                                                                          <  1:297112/300‑1 (MQ=255)                                                                                                                              aGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGTAGCGGCCTTGCGGCGCGACTGGGCTGCAAagagatcg                                                                                                                          >  2:730550/1‑295 (MQ=255)                                                                                                                                      aaaCCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAtt                                                                                                                  <  1:480499/301‑1 (MQ=255)                                                                                                                                       aaCCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt                                                                                                                 >  1:483389/1‑301 (MQ=255)                                                                                                                                                          aaCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATTGGCTGCAAAGCAGGCCCGGAGATTTcgacgattgggatgaaaaa                                                                                              >  2:490757/1‑283 (MQ=255)                                                                                                                                                                     ccgtgccccccgcgTGGAGCCCGGCGTCCCCGGCGAGGCGGCCTACGGGCATGGCGGGTTGCAGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGTTATGGAGGACGACCATAAAGTCCGGGGGGCAGCCGATTACCGGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGGATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCtt                                                                                   <  2:262756/291‑1 (MQ=255)                                                                                                                                                                     ccaTCTCCTACGCGTCGATCCCGGCGTCCCCGTCGGCGCCGCCTACGGGCATGTCTCGATCCTGCTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCtt                                                                                   <  2:136871/299‑1 (MQ=255)                                                                                                                                                                              cggcgTAGAGCCCGGCGTCCCCGTCGACGCGGCGCACGAGCAGCTCTAGATGCTGATGGGCTACATCAAGCGCCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCGGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCAGAGGAACGATATCGAGATCGCCAAGGGTAAGGCGCATGAGGCATGGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCGGAAAACGGGGGCCGCTTTGCGGCCCt                                                                          <  2:345174/299‑1 (MQ=255)                                                                                                                                                                                   tAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCGGATGGGCTACATCAAGCACCTGGTGTGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGGCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTcgc                                                                     <  1:467272/301‑1 (MQ=255)                                                                                                                                                                                                                                    gATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAAAGGGGGACGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACCAAAGCCGCTTGACCTAACCACCAGAAATg                    >  1:493829/1‑301 (MQ=255)                                                                                                                                                                                                                                       gggCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCGAGGCCGCTCCCACAAAAGCCGCTTGCGCTAAACACAAGAAATGttt                 >  1:387204/1‑301 (MQ=255)                                                                                                                                                                                                                                         cctacaaaaagaagcgggTGAGCGAGGGGGATAGGGAGGAGGAGCATAAATTCCTGGGAGCAGACGATTACATGTGCGACATGGCCAAGGCGCTGAAGAACGATATCGGGATCGCGAAGGGTAAGGCGCATTAGGCAGTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCCGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAACCGGGGACCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTAt               <  1:502767/285‑1 (MQ=255)                                                                                                                                                                                                                                          aaacataaatgacccggTGGGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATGGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGGCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATg              <  2:106131/286‑1 (MQ=255)                                                                                                                                                                                                                                                 aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGAGGACAGGCCCCGCCAACAAAAGCCCCTCGCGCTAACCCAAAGAAATGTTTATgaccacgt       >  2:320945/1‑294 (MQ=255)                                                                                                                                                                                                                                                      ccTGGTGCGCGAGGGTGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGCCGCATTAGGCATTGACGCGGTCATTTTTGGGAGCTGCCTTGCGTCGCGAATGGGCTTGAAATCAGTCCTGGCGACTTAAACGATTGCGCTGAAAACGGGGGCCGCTTTGCAGCCCTTTTCCGACGCGAGGCAGCACACAAAAAAGCCGCTTGCGCGAACCAACAGAAATGTTTAAGAACAATTGGCaa  >  2:832676/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       TTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACAACCGCAGGCGTGGAAACCTTCCTCGCAGCCGTTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCAT  >  AE015451/335826‑336370

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