Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	197,496	T→C	14.6%	I1001T (ATC→ACC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	197,496	0	T	C	14.6%	134.0 / 11.3	60	I1001T (ATC→ACC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (18/33);  new base C (5/4);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCA  >  AE015451/197263‑197657                                                                                                                                                                                                                                          |                                                                                                                                                                  aGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATc                                                                                                                                                                  >  1:1008362/1‑235 (MQ=35)                 ggCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACACCAACGGTTCGCCCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGCACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCCCCGCGGCCCTGACGGCGAGCCCGTCGGTCACCGAAGGTGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCa                                                                                                                                                 <  1:181041/236‑1 (MQ=255)                 ggCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCCACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCa                                                                                                                                                 <  1:259072/236‑1 (MQ=18)                       gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  1:734316/236‑1 (MQ=14)                         cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCgg                                                                                                                                         >  1:864441/1‑236 (MQ=34)                          gATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:864438/233‑1 (MQ=34)                              tCCAGACCCCGGCCAATGACG‑CTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAAc                                                                                                                                     >  2:700719/1‑234 (MQ=37)                               ccAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGACCCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGccc                                                                                                                                                  <  1:100483/221‑1 (MQ=21)                                                      acaacaACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTg                                                                                                                                                               >  1:68101/1‑185 (MQ=33)                                                               gTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAGCACCTACACCGTGATcc                                                                                                                                                                 >  1:117376/1‑174 (MQ=18)                                                                     ccGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCag                                                                                                          >  1:670031/1‑222 (MQ=37)                                                                      cGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  1:700722/236‑1 (MQ=37)                                                                           gCGTCACCACCGAGAACGCCACAGGCGGCAACTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:28519/184‑1 (MQ=14)                                                                            cGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAAc                                                                                                             >  2:989695/1‑213 (MQ=37)                                                                            cGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAAc                                                                                                             <  1:989698/213‑1 (MQ=37)                                                                                                  ggcggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCa                                                                                                                                                 <  2:563757/155‑1 (MQ=32)                                                                                                       cAAGTTCGACCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGCCCCTGACGGCGAGCCCGTCGGTCACCGCAGGCGGCGTGATCACCTACACCGTGATCCTGa                                                                                                                                                              <  2:517272/137‑1 (MQ=11)                                                                                                           ttCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCg                                                         >  2:909644/1‑234 (MQ=35)                                                                                                            tCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCCCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCg                                                                                                                                   <  1:257226/159‑1 (MQ=18)                                                                                                                  aacTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:1059980/143‑1 (MQ=17)                                                                                                                          ccGCACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  1:136468/184‑1 (MQ=25)                                                                                                                               cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTc                                                                                                                         <  2:314799/150‑1 (MQ=37)                                                                                                                               cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTc                                                                                                                         >  1:314802/1‑150 (MQ=17)                                                                                                                                 cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGcc                                                                                                                                                   <  1:249743/122‑1 (MQ=32)                                                                                                                                     cccGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  2:1077393/172‑1 (MQ=25)                                                                                                                                             aGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGc                                                                                              >  1:141105/1‑163 (MQ=35)                                                                                                                                               aCCCCGATCAACGCCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                    <  2:355523/187‑1 (MQ=14)                                                                                                                                               aCCACGATCAACGACTCGGTCAATACCACCACCGCGAACCCGACGGCGAGCCCCTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGcgg                                                                                            <  1:952204/163‑3 (MQ=11)                                                                                                                                                  aCGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGACGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            >  1:588854/1‑160 (MQ=9)                                                                                                                                                                ggTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGcc   <  1:909647/235‑1 (MQ=37)                                                                                                                                                                   cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 >  1:102897/1‑138 (MQ=37)                                                                                                                                                                   cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGc                                                                                                                  >  2:658833/1‑121 (MQ=35)                                                                                                                                                                   cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGc                                                                                                                  <  1:658836/121‑1 (MQ=35)                                                                                                                                                                   cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 <  2:102895/138‑1 (MQ=37)                                                                                                                                                                   cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGcgg                                                                                            >  1:529216/1‑141 (MQ=38)                                                                                                                                                                   cGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGcgg                                                                                            <  2:529213/143‑3 (MQ=16)                                                                                                                                                                     acgCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  1:529446/92‑1 (MQ=18)                                                                                                                                                                     acaccactaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           >  2:529443/4‑94 (MQ=9)                                                                                                                                                                        ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTac                                          <  1:141711/188‑1 (MQ=35)                                                                                                                                                                               cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCca  >  1:302187/1‑221 (MQ=17)                                                                                                                                                                               cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                   <  2:943474/172‑1 (MQ=14)                                                                                                                                                                                gcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  1:357482/83‑1 (MQ=33)                                                                                                                                                                                    cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:621841/79‑1 (MQ=34)                                                                                                                                                                                    cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                 <  2:396659/89‑1 (MQ=38)                                                                                                                                                                                    cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGCTGactct                                                                                             >  2:142786/1‑121 (MQ=255)                                                                                                                                                                                      cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCt                                                                                                                                                                <  1:747210/56‑1 (MQ=38)                                                                                                                                                                                      cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:512419/77‑1 (MQ=34)                                                                                                                                                                                      cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGacc                                                   >  2:141709/1‑163 (MQ=14)                                                                                                                                                                                         aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:977832/74‑1 (MQ=35)                                                                                                                                                                                         aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           >  1:977835/1‑74 (MQ=35)                                                                                                                                                                                                   cGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 >  2:232097/1‑106 (MQ=17)                                                                                                                                                                                                   cGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 <  1:232100/106‑1 (MQ=17)                                                                                                                                                                                                             cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 <  2:105615/96‑1 (MQ=37)                                                                                                                                                                                                              gAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  1:648563/100‑1 (MQ=40)                                                                                                                                                                                                              gAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCg                                                                                            >  2:648560/1‑100 (MQ=11)                                                                                                                                                                                                                      cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACg                                                                                                                                             <  2:169676/43‑1 (MQ=38)                                                                                                                                                                                                                          aTCACTTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  1:405265/88‑1 (MQ=255)                                                                                                                                                                                                                             aCCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:380702/38‑1 (MQ=40)                                                                                                                                                                                                                                        gATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAcc                                                                                                                              >  2:879729/1‑40 (MQ=40)                                                                                                                                                                                                                                        gATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgag         >  2:649319/1‑157 (MQ=25)                                                                                                                                                                                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAcc                                                                                                                              >  1:380705/1‑38 (MQ=33)                                                                                                                                                                                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTTCCGTTGCAGCCGGCCAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGACTGa                                                 >  1:1010934/1‑115 (MQ=11)                                                                                                                                                                                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTAcc                                                                                                      <  1:822223/62‑1 (MQ=39)                                                                                                                                                                                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCg                              >  1:865401/1‑134 (MQ=37)                                                                                                                                                                                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCg                              <  2:865398/134‑1 (MQ=35)                                                                                                                                                                                                                                          tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGACAAAACTATTACCGTTGACGCGGGCCATACCCAggg                                                                               >  1:1052492/1‑85 (MQ=11)                                                                                                                                                                                                                                          |                                                                                                                                                                  AGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCA  >  AE015451/197263‑197657

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