Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	196,327	A→C	28.5%	T611T (ACA→ACC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	196,327	0	A	C	28.5%	94.1 / 31.2	61	T611T (ACA→ACC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base C (11/7);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACC  >  AE015451/196132‑196534                                                                                                                                                                                                     |                                                                                                                                                                                                                ggTTAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         <  1:546156/236‑1 (MQ=255) ggTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         <  1:938279/236‑1 (MQ=255)                cGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAccc                                                                                                                                                                                                          >  1:363793/1‑188 (MQ=25)                    aaCGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGa                                                                                                                                                                                                             >  2:901132/1‑181 (MQ=255)                         cACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCa                                                                                                                                                                                        <  2:357079/197‑1 (MQ=40)                         cACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCa                                                                                                                                                                                        >  1:357082/1‑197 (MQ=255)                               cGGTAATTTCGACCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         <  2:354266/206‑1 (MQ=38)                                           gCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGc                                                                                                                                                                       >  1:65763/1‑196 (MQ=38)                                           gCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                                  <  2:563327/233‑1 (MQ=21)                                                     ccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                <  2:695321/225‑1 (MQ=14)                                                          cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGa                                                                                                                                                   >  1:435113/1‑201 (MQ=14)                                                                 ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                                  <  1:646368/147‑1 (MQ=25)                                                                    gcccAGACCACCATCAACGACCCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                                  <  2:360745/205‑1 (MQ=9)                                                                        aGACCACGATCACCGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         <  2:124971/165‑1 (MQ=18)                                                                          aCCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                                             <  1:310645/223‑1 (MQ=21)                                                                                 tCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaaga                                                                                                                                                                   >  1:523256/1‑160 (MQ=39)                                                                                        cccGGTCGACACCACACCCGCGACCCTGACGGCGAGCCCGGCGGTCACCGGAGGCGGCGTGATCACTTACCCCCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                                  <  2:167397/123‑1 (MQ=11)                                                                                          cGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                          <  2:324567/130‑1 (MQ=25)                                                                                                acacCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                         >  1:604269/1‑173 (MQ=21)                                                                                                          cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGc                                                                                                                                                                                               <  2:718760/109‑1 (MQ=18)                                                                                                          cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGc                                                                                                                                                                                               >  1:718763/1‑109 (MQ=17)                                                                                                          cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         <  2:348464/131‑1 (MQ=32)                                                                                                          cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAttt                                                                  >  1:606715/1‑234 (MQ=9)                                                                                                          cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                                 >  1:809718/1‑235 (MQ=12)                                                                                                          cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         >  1:348467/1‑131 (MQ=17)                                                                                                              aCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCAca                                                                             >  1:925372/1‑218 (MQ=17)                                                                                                                 cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                                                                                                                                 >  1:599033/1‑100 (MQ=33)                                                                                                                               gTCGGACACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTc                                                                                                                                                     >  2:778978/1‑130 (MQ=255)                                                                                                                                 cGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAGACCATCACCgg                                                                                                                                                                                 >  2:78331/1‑99 (MQ=11)                                                                                                                                           aGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                                      <  2:809715/213‑1 (MQ=2)                                                                                                                                            ggcggcGTGATCTCTTA‑CACCGAGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAg                                                                                                                                                                                                                <  1:844858/58‑1 (MQ=11)                                                                                                                                            ggcggcGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                                                                                                  <  1:852334/72‑1 (MQ=34)                                                                                                                                               ggagTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAcc                                                                                                                                                                                         <  1:614432/75‑1 (MQ=21)                                                                                                                                                 cgTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         <  1:467262/92‑1 (MQ=33)                                                                                                                                                 cgTGATCACTTA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                         >  2:467259/1‑92 (MQ=16)                                                                                                                                                    gATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAGACCATCACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGActgttc                                                                                                                        >  1:819404/1‑133 (MQ=11)                                                                                                                                                        aCTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGAc                                                                                                         >  1:460395/1‑149 (MQ=34)                                                                                                                                                          ttA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAt                                                                    <  1:855368/184‑1 (MQ=2)                                                                                                                                                          gtA‑CACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCg                                      <  1:754575/213‑1 (MQ=255)                                                                                                                                                                 cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGTAAAAc                                                                                                                                                                  >  1:4065/1‑84 (MQ=18)                                                                                                                                                                  gTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcg                                                                            >  2:938975/1‑169 (MQ=37)                                                                                                                                                                      cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                                                                                                                        <  2:682897/89‑1 (MQ=35)                                                                                                                                                                      cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                                                                                                                        >  1:682900/1‑89 (MQ=33)                                                                                                                                                                      cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGt                                                                                                                                                      <  2:170909/91‑1 (MQ=33)                                                                                                                                                                      cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGt                                                                                                                                                      >  1:170910/1‑91 (MQ=33)                                                                                                                                                                      cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGc      >  1:724504/1‑235 (MQ=39)                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGAc                                                                                                                             <  2:56664/116‑1 (MQ=11)                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGAc                                                                                                                             >  1:56664/1‑116 (MQ=11)                                                                                                                                                                        cTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                                                                                                                        <  2:407059/87‑1 (MQ=255)                                                                                                                                                                        cTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                       <  2:674631/88‑1 (MQ=255)                                                                                                                                                                        cTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                       >  1:674634/1‑88 (MQ=255)                                                                                                                                                                        cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCg                                                                                                                                   >  2:552513/1‑108 (MQ=11)                                                                                                                                                                              aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccact  <  1:752166/231‑2 (MQ=12)                                                                                                                                                                              aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccact  >  2:752163/1‑230 (MQ=16)                                                                                                                                                                              aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGGAACTGACCCCgc                                                 >  1:96712/1‑183 (MQ=37)                                                                                                                                                                                  cTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCCGCAAAACCCAGGGCAg                                                                                                                                                         >  1:669217/1‑76 (MQ=21)                                                                                                                                                                                   tGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcca                                                                                                                                  >  1:721847/1‑96 (MQ=255)                                                                                                                                                                                   tGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCaaaa                                                                                                                                                                   >  2:930293/1‑65 (MQ=18)                                                                                                                                                                                       cAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTc               <  1:15724/209‑1 (MQ=255)                                                                                                                                                                                        aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCCTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGc                                                                                 >  1:716361/1‑142 (MQ=255)                                                                                                                                                                                                 tGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAc                                                                                                                                          <  2:669214/76‑1 (MQ=34)                                                                                                                                                                                                     |                                                                                                                                                                                                                GGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA‑CACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACC  >  AE015451/196132‑196534

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