Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	215,395	A→G	10.4%	P6967P (CCA→CCG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	215,395	0	A	G	10.4%	138.9 / 4.4	56	P6967P (CCA→CCG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base G (3/3);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACA  >  AE015451/215168‑215606                                                                                                                                                                                                                                    |                                                                                                                                                                                                                    aCGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACAccc                                                                                                                                                                                                              <  1:320025/235‑1 (MQ=33)               gCCGGTCAGCACGTCGGTGACCGACGAGACAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                <  1:635644/235‑1 (MQ=25)                          cGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCa                                                                                                                                                                                                             >  1:83272/1‑211 (MQ=32)                              ggTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                <  1:580052/236‑1 (MQ=39)                                         aGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTg                                                                                                                                                                                                        >  2:502008/1‑201 (MQ=33)                                                               aGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAAGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                                                                                                                                <  1:502009/235‑1 (MQ=39)                                                                                 cACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACgg                                                                                                                               <  2:557315/234‑1 (MQ=38)                                                                                        cccgccgACCAGACCTCGGGGGCCGAGAACGTCAAGCCGACCTTCACCGGGCCCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGc                                                                                                                                                               <  2:192001/194‑1 (MQ=255)                                                                                          cgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCa                                                                                                                                                                 <  2:556255/191‑1 (MQ=9)                                                                                           gccgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGa                                                                                                                   >  1:597683/1‑236 (MQ=2)                                                                                            ccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAACGCCGGCGAAACCAGCGCGACATACACCCACGCTgcg                                                                                                                                                                                                      >  2:215754/1‑152 (MQ=18)                                                                                            ccgACCAGACCTCGGTGGCCGAGAACGTCAAGACGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGTCGAAACCAGCGCGCCATaca                                                                                                                                                                                                                 >  2:461636/1‑141 (MQ=21)                                                                                             cgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                >  2:706513/1‑157 (MQ=25)                                                                                             cgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                                                                                                                                                <  1:706514/157‑1 (MQ=38)                                                                                               aCCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAAc                                                                                                               >  1:719689/1‑236 (MQ=18)                                                                                                         ggTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                      <  1:340110/235‑1 (MQ=37)                                                                                                            ggCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                     <  1:99460/153‑1 (MQ=35)                                                                                                            ggCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                                                                                                                                                     >  2:99460/1‑153 (MQ=35)                                                                                                                 agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCGGCTGTGGACGTAg                                                                                                                                   >  1:394385/1‑198 (MQ=21)                                                                                                                   aaCGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGc                                                                                              <  2:665148/233‑1 (MQ=34)                                                                                                                     cgccAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTCGCCCACGACCTGGACGTGATCCTGAGCAACCACGCCCAGGTTACCATCAGGGCCGGCGAAACGAGCCCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                <  1:493302/146‑1 (MQ=11)                                                                                                                     cGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGATGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTa                                                                                         <  1:55236/236‑1 (MQ=18)                                                                                                                        cAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTa                                                                                         <  2:328526/233‑1 (MQ=34)                                                                                                                         aaGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGt                                                                                                                                                  <  2:547815/175‑1 (MQ=35)                                                                                                                                 cctcgccGTGCACGCCAACCCGCCGCTGGCCCACGACCTGGTCGTGCCCCTGAGCAACAACGCCCAGGTCGCCAGCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAacg                                                                                                                                                                                <  2:252749/132‑1 (MQ=11)                                                                                                                                 cTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGt                                                                               <  2:452964/234‑1 (MQ=33)                                                                                                                                 cTTCACCGGGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                               <  2:558078/170‑1 (MQ=21)                                                                                                                                      ccGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCt                                                                                                                                                                           <  1:405252/137‑1 (MQ=37)                                                                                                                                            aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTg                                                                                                                                                                                                        <  2:14207/102‑1 (MQ=17)                                                                                                                                            aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTg                                                                                                                                                                                                        >  1:14207/1‑102 (MQ=35)                                                                                                                                            aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTCGACGGCCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGCGAc                                                                            >  1:124980/1‑226 (MQ=17)                                                                                                                                            aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTAGa                                                                                                                                        >  2:508186/1‑166 (MQ=25)                                                                                                                                             cGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAGCAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACg                                                                                                                                                   <  2:99314/154‑1 (MQ=18)                                                                                                                                              gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTAt                                                                                                                                                       >  2:569767/1‑149 (MQ=25)                                                                                                                                              gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACcga                                                                >  1:334852/1‑236 (MQ=33)                                                                                                                                                 aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTcc                                                                                                                                                >  2:467949/1‑153 (MQ=11)                                                                                                                                                 aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATaccacc                                                                  >  2:54249/1‑231 (MQ=32)                                                                                                                                                         gcaggCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                                                                                                                                               <  2:534694/143‑1 (MQ=14)                                                                                                                                                                cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATa                                                                                                                                                                                                                   <  1:613766/71‑1 (MQ=34)                                                                                                                                                                           cGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATacc                                                                     >  1:580500/1‑202 (MQ=32)                                                                                                                                                                              gACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg                                                                                                      <  1:661480/166‑1 (MQ=255)                                                                                                                                                                                  cTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGt                            >  1:254725/1‑236 (MQ=32)                                                                                                                                                                                      gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGc                                                                                           >  2:89491/1‑169 (MQ=33)                                                                                                                                                                                      gCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTCGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGCAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGt                                                                               >  2:96671/1‑181 (MQ=18)                                                                                                                                                                                       caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTctg                                                                                                                                               >  2:623814/1‑114 (MQ=11)                                                                                                                                                                                       caacaaCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGGGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGa                                            >  2:73630/1‑215 (MQ=17)                                                                                                                                                                                             cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAg                    >  2:560669/1‑233 (MQ=32)                                                                                                                                                                                             cGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCCGCTGTGGAGGTAGACGGGCGCACCTTCGAAAAcaggca                                                                                                          >  2:231695/1‑142 (MQ=255)                                                                                                                                                                                                ccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGc               <  1:724201/235‑1 (MQ=14)                                                                                                                                                                                                ccAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCg              <  1:54249/236‑1 (MQ=31)                                                                                                                                                                                                           tCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCGGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAccc                                   >  2:352473/1‑204 (MQ=11)                                                                                                                                                                                                           tCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCGGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGAccc                                   <  1:352474/204‑1 (MQ=11)                                                                                                                                                                                                             aaGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGtt                                                                                    <  2:52408/153‑1 (MQ=31)                                                                                                                                                                                                                gCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTaca  <  1:680810/232‑1 (MQ=31)                                                                                                                                                                                                                    gCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGAc                                           <  1:319917/187‑1 (MQ=31)                                                                                                                                                                                                                     cGAACCCAGCGCGCCATACACCCGCGCTGCGCAAGGCGATGACGTCTACAGCGACGCTGGCCAGATCAGCCTGGGCCTCACGTCCGCTGTGGATGTAGACGGTCGCGCCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAg                                                <  1:231695/181‑1 (MQ=2)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                    ACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTGGGCGGTGCGGCTACCGTTCAAGTGACCGATACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTCACCGAAGGTGGCGAGATCACTTACA  >  AE015451/215168‑215606

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