Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	203,470	T→C	9.0%	V2992V (GTT→GTC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	203,470	0	T	C	9.0%	151.9 / 3.5	54	V2992V (GTT→GTC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base C (3/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACC  >  AE015451/203243‑203693                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                gTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATTACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGATCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGTGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAacg                                                                                                                                                                                                                          <  1:878314/235‑1 (MQ=255)               gggCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCa                                                                                                                                                                                                           <  1:454944/236‑1 (MQ=40)                   aCCCCAGGCAACGAAGGCGACTTGGCCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTgg                                                                                                                                                                                                                   <  1:173015/224‑1 (MQ=39)                                        ttGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGt                                                                                                                                                                                     <  2:908929/233‑1 (MQ=35)                                              aaGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCa                                                                                                                                                                             >  2:636222/1‑235 (MQ=37)                                                   cACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTgg                                                                                                                                                                                                                   >  1:301468/1‑192 (MQ=37)                                                          acaGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCtt                                                                                                                                                                <  1:818391/236‑1 (MQ=35)                                                            aGCCGACCAGACCGCGGTGGCTGAGAACGTCAACCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTg                                                                                                                                                                                    <  1:156217/214‑1 (MQ=18)                                                              ccGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGc                                                                                                                                                                                                                      <  1:109559/178‑1 (MQ=14)                                                                    aGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGCGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCTCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGtt                                                                                                                                                                                    >  1:233738/1‑205 (MQ=21)                                                                      aCCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGc                                                                                                                                                                               >  2:728739/1‑209 (MQ=35)                                                                                   agaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCTCAGGGTGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                                                                <  1:720257/163‑1 (MQ=17)                                                                                   agaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCTCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCa                                                                                                                                                                                                           >  2:537678/1‑168 (MQ=18)                                                                                   agaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTcggcg                                                                                                                                                                                       >  2:868705/1‑188 (MQ=34)                                                                                   agaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTcggcg                                                                                                                                                                                       <  1:868708/188‑1 (MQ=33)                                                                                   agaACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                                                                >  2:720254/1‑163 (MQ=14)                                                                                     aaCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGGGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGc                                                                                                                                                                               >  2:246534/1‑194 (MQ=18)                                                                                      aCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGAGCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCcgcg                                                                                                                                    >  1:575307/1‑236 (MQ=25)                                                                                             aCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCTCAGGGTGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                                                                <  2:762611/153‑1 (MQ=11)                                                                                             aCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCa                                                                                                                                                                                                                >  1:762614/1‑153 (MQ=11)                                                                                                                    aCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcgacga                                                                                                       <  1:320514/235‑1 (MQ=37)                                                                                                                     ccAGCCGCGGGCCCACGACCTGGTCGTGACCCTGGGCAACAATGACCAGGTCACAATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCAACCTTCgaga                                                                                                                                                           <  1:28290/182‑1 (MQ=9)                                                                                                                         ccgccgGCCCACGGCATGGCCGTGACCCTGAGCAACAATGCCCAGGTCACAATCAAGGCTGGCGAAACCACCGCGCCGC‑‑AGACCCGCACTACGCAGGGTGATGACGATAACAACGGCGCAGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTg                                                                                                                                                                          <  1:536298/163‑1 (MQ=11)                                                                                                                                              gTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCAccc                                                                             >  1:833107/1‑235 (MQ=37)                                                                                                                                                    cTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCtt                                                                                                                                                                >  1:198113/1‑146 (MQ=32)                                                                                                                                                    cTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGt                                                                      >  1:375328/1‑236 (MQ=35)                                                                                                                                                          aacaaTGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGt                                                                                                                               >  2:785327/1‑173 (MQ=38)                                                                                                                                                                gCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggc                                                            >  1:671282/1‑236 (MQ=21)                                                                                                                                                                 cccAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                         <  1:333079/236‑1 (MQ=34)                                                                                                                                                                  ccAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcgg                                                          >  2:333076/1‑234 (MQ=34)                                                                                                                                                                        caccacCCAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAgg                                                             <  1:922891/225‑1 (MQ=21)                                                                                                                                                                              cAAGGCCGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGgcggcg                                                        >  2:889726/1‑224 (MQ=11)                                                                                                                                                                              cAAGGCCGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGgcggcg                                                        <  1:889729/224‑1 (MQ=14)                                                                                                                                                                               aaGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                 <  1:134898/230‑1 (MQ=34)                                                                                                                                                                               aaGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTacac                                           <  1:299199/236‑1 (MQ=35)                                                                                                                                                                                 ggCTGGCGAACCCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACACCGACGCTGGCCAGATCAGCCTGGGCATCGCCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCCCCGACGAAGTCGTGGCCAAGTTGACCGCCAcc                                                                              <  1:152629/199‑1 (MQ=255)                                                                                                                                                                                 ggCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCcgc                                                                                                                                     >  2:856656/1‑144 (MQ=32)                                                                                                                                                                                  gCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                         <  1:220792/219‑1 (MQ=34)                                                                                                                                                                                               gcgcgcCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCcgcg                                                                                                                                    >  2:220790/1‑131 (MQ=33)                                                                                                                                                                                                                aCGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                 <  2:571443/199‑1 (MQ=33)                                                                                                                                                                                                                aCGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCg            <  1:535642/236‑1 (MQ=38)                                                                                                                                                                                                                aCGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCg            <  1:871556/236‑1 (MQ=39)                                                                                                                                                                                                                 cGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTg                                 <  1:554573/214‑1 (MQ=38)                                                                                                                                                                                                                 cGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTg                                 >  2:554570/1‑214 (MQ=37)                                                                                                                                                                                                                  gCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt          <  1:636225/236‑1 (MQ=40)                                                                                                                                                                                                                    tGCGCAGGGTGATGACGTTTACAACGAGGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATca        >  1:493022/1‑236 (MQ=255)                                                                                                                                                                                                                    tGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATc         <  1:730036/235‑1 (MQ=40)                                                                                                                                                                                                                          gggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCaacaac   <  1:615137/235‑1 (MQ=39)                                                                                                                                                                                                                           ggTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAAcc  <  1:57396/235‑1 (MQ=38)                                                                                                                                                                                                                              gatgaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGaca                                                                                                                 <  1:269172/121‑1 (MQ=35)                                                                                                                                                                                                                                 gaCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCCTCGCGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGt                                                                                                    >  2:50459/1‑131 (MQ=21)                                                                                                                                                                                                                                   cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcgacga                                                                                                       <  2:7240/126‑1 (MQ=35)                                                                                                                                                                                                                                   cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcgacga                                                                                                       >  1:7240/1‑126 (MQ=35)                                                                                                                                                                                                                                   cGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGt                                                                                         >  2:843354/1‑140 (MQ=34)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                GTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGT‑‑ACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACC  >  AE015451/203243‑203693

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