Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	211,783	T→C	8.9%	D5763D (GAT→GAC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	211,783	0	T	C	8.9%	187.0 / 6.0	66	D5763D (GAT→GAC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base C (2/4);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGC  >  AE015451/211551‑212018                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            gcaacgCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAAc                                                                                                                                                                                                                                          >  1:333091/6‑236 (MQ=255)             gtCAAGGGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCGTCGGTCACCGAAGGTGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGt                                                                                                                                                                                                                                  <  2:633972/231‑1 (MQ=255)                              accGACGAAGTAGTGGCCAAGTTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTa                                                                                                                                                                                                                                     <  2:333090/212‑1 (MQ=255)                                  acgaAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTg                                                                                                                                                                                                            >  2:985394/1‑233 (MQ=255)                                       gtggtgGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTc                                                                                                                                                                                                       <  2:402616/233‑1 (MQ=255)                                                         gcgACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGGCTACGt                                                                                                                                                                                                                                  >  1:524160/1‑188 (MQ=255)                                                          cgACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCaaag                                                                                                                                                                                                                       >  1:54771/1‑196 (MQ=255)                                                               ccTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGt                                                                                                                                                                                    >  1:1135027/1‑228 (MQ=40)                                                                          cGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTcc                                                                                                                                                                 >  2:796969/1‑235 (MQ=39)                                                                                  gtgAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTg                                                                                                                                                         <  1:1054493/236‑1 (MQ=35)                                                                                  gcgagaTCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACCGTCACTGCCCCGGATAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                                                     <  1:393956/190‑1 (MQ=255)                                                                                      gaTCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGaca                                                                                                                                                      >  1:653801/1‑235 (MQ=38)                                                                                          aCCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTg                                                                                                                                                   >  2:959532/1‑234 (MQ=38)                                                                                            cTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTg                                                                                                                                                                                         >  2:524532/1‑194 (MQ=37)                                                                                             tACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTTATGTca                                                                                                                                                                                                      >  1:327320/1‑179 (MQ=25)                                                                                                       cccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGt                                                                                                                                     >  1:450061/1‑235 (MQ=35)                                                                                                       cccTGACCAACAAAGATGGTCTGCCGACCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGt                                                                                                                                                                                    >  1:1095713/1‑188 (MQ=18)                                                                                                       cccTGACCAACAAAGATGGTCTGCCGAACAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTc                                                                                                                                    >  1:597717/1‑236 (MQ=25)                                                                                                        ccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                                                   >  1:902246/1‑172 (MQ=34)                                                                                                        ccTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                                                   <  2:902245/172‑1 (MQ=34)                                                                                                             ccAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCAc                                                                                                                               >  1:576703/1‑235 (MQ=35)                                                                                                               aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAg                                                                                                                                  >  1:104575/1‑230 (MQ=35)                                                                                                                acaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTc                                                                                                                                    <  2:450060/227‑1 (MQ=35)                                                                                                                         gccTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGACtgtgatact                                                                                                                                                                              >  1:148838/3‑171 (MQ=11)                                                                                                                           cTGCCGATCAACAACCACAGTGCACTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCgg                                                                                                                          <  1:1102371/226‑1 (MQ=21)                                                                                                                            tGCCGACCACCAACCACAGTGCATTGACCTTCACCCGGAGCGACGGCAAAAGCGTCATCCCCGTGCCGGCTAAGGGCACCGTGGGTACTGCCACTGTCGCTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGt                                                                                                                                                                            <  2:384231/175‑1 (MQ=11)                                                                                                                              ccGATCAACATCCACGGTGCATTGACCTTCACGCTGAGCGACGGCAAAGCCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGAt                                                                                                                                                                                                          <  2:148838/143‑1 (MQ=255)                                                                                                                              ccGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGGCGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGt                                                                                                                                                                                                                                  <  2:327319/119‑1 (MQ=17)                                                                                                                              ccGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGc                                                                                                              >  1:644207/1‑235 (MQ=34)                                                                                                                                 ttCAACAACCCCAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCg                                                                                                                                                                                      <  2:871359/159‑1 (MQ=14)                                                                                                                                  tCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAGCCGTCATCACCGTGCCGGCTAATGGCACCGTGGTTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGc                                                                                                                                 <  1:20862/212‑1 (MQ=18)                                                                                                                                         ccACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCAc                                                                                                                               <  2:1097002/207‑1 (MQ=33)                                                                                                                                         ccACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGcc                                                                                                    <  2:525298/234‑1 (MQ=33)                                                                                                                                           acaGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCggg                                                                                                 >  1:344489/1‑235 (MQ=32)                                                                                                                                                 cATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGCCACTgg                                                                                                                                                  >  1:1880/1‑180 (MQ=14)                                                                                                                                                 cATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACg                                                                                            >  1:725571/1‑234 (MQ=32)                                                                                                                                                              ccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTc                                                                              >  1:323775/1‑235 (MQ=32)                                                                                                                                                              ccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCa                                                                             >  1:691954/1‑236 (MQ=33)                                                                                                                                                              ccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGGACGCCGGGCAACGAAGGCGACCTg                                                                                 >  1:130578/1‑232 (MQ=17)                                                                                                                                                              ccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTG‑CCGGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCaa                                                                            >  1:42341/1‑236 (MQ=17)                                                                                                                                                               cTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCaa                                                                            <  1:796970/236‑1 (MQ=33)                                                                                                                                                                         gCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCAcc                                                                     >  2:902644/1‑233 (MQ=33)                                                                                                                                                                          cAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCa                                                                    <  2:653800/233‑1 (MQ=33)                                                                                                                                                                           aaaaCCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCg                                                                                      <  1:416393/214‑1 (MQ=33)                                                                                                                                                                                       aCCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCg                                                                                      >  2:416392/1‑202 (MQ=33)                                                                                                                                                                                       aCCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCa                                                       >  2:1105884/1‑233 (MQ=34)                                                                                                                                                                                       aCCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGa                                                     <  1:902645/235‑1 (MQ=35)                                                                                                                                                                                                   aaTGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTc                                                                                                                                    <  2:994032/144‑1 (MQ=35)                                                                                                                                                                                                   aaTGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTc                                                                                                                                    >  1:994033/1‑144 (MQ=34)                                                                                                                                                                                                      ggCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGaa                                     <  1:959533/236‑1 (MQ=39)                                                                                                                                                                                                        cACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAAc                                    >  1:890732/1‑235 (MQ=37)                                                                                                                                                                                                          ccGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGc                                                             <  2:323774/208‑1 (MQ=34)                                                                                                                                                                                                           cGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGcc                                                            <  2:950288/208‑1 (MQ=35)                                                                                                                                                                                                             tGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGGCTTCGGTAGCCGAGAACGTCaaa                              >  1:1052811/1‑236 (MQ=37)                                                                                                                                                                                                              cggcaCTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGCACTAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTc                                 <  1:977115/228‑1 (MQ=11)                                                                                                                                                                                                                gTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGc                                           <  2:875798/220‑1 (MQ=39)                                                                                                                                                                                                                        aCTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGAccagccag                                                                                                           >  2:925332/1‑148 (MQ=25)                                                                                                                                                                                                                        aCTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACTCCGGGCAACGAAGGCGACCTGGTc                                                                              >  1:1071900/1‑177 (MQ=18)                                                                                                                                                                                                                         ccgTCACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGa                                                          <  2:1099291/194‑1 (MQ=25)                                                                                                                                                                                                                          cgTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTg                 <  1:499962/235‑1 (MQ=255)                                                                                                                                                                                                                               cTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTg                                                                                 <  2:989248/167‑1 (MQ=25)                                                                                                                                                                                                                               cTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTg                                                                                 >  1:989249/1‑167 (MQ=21)                                                                                                                                                                                                                                   cccGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGt                                                                                                         >  1:759836/1‑139 (MQ=37)                                                                                                                                                                                                                                   cccGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGt                                                                                                         <  2:759835/139‑1 (MQ=37)                                                                                                                                                                                                                                   cccGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCCCCATCACGGCCGACCa                                                       >  2:321727/1‑189 (MQ=25)                                                                                                                                                                                                                                   cccGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACAGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCGGGTACGCCGGGCATCGAAGGCGAGCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGTTGGCCGAGAACGTCAAACCGACCTTCACTGTGCAc            >  2:1044009/1‑232 (MQ=11)                                                                                                                                                                                                                                         tAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGc  <  1:75020/236‑1 (MQ=255)                                                                                                                                                                                                                                         tAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGGCACCATCACGGCCGACCAGACTTCGGTat                                            >  1:273884/1‑193 (MQ=37)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            GTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGC  >  AE015451/211551‑212018

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