Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	202,342	C→G	8.9%	T2616T (ACC→ACG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	202,342	0	C	G	8.9%	133.9 / 3.4	44	T2616T (ACC→ACG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base G (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAGC  >  AE015451/202127‑202576                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            ccccACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGAc                                                                                                                                                                                                                              <  1:143471/230‑1 (MQ=37) acCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                                                                          <  1:744455/235‑1 (MQ=39)    cACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCt                                                                                                                                                                                                                         >  2:323370/1‑233 (MQ=40)              aaGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaacaa                                                                                                                                                                                                                                  >  2:356001/1‑214 (MQ=39)                                  gACGCTGGCCAGATCAGCTTGGGCATCACGTCGGCGGTGGATGTAGACGGTCCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                                                                                                                                                                                             <  1:269107/183‑1 (MQ=11)                                     gCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTg                                                                                                                                                                                                                        <  2:436937/201‑1 (MQ=35)                                     gCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTg                                                                                                                                                                                                                        >  1:436938/1‑201 (MQ=38)                                      cTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTCCCGATCAACAACCACTCTGAGCTGTAc                                                                                                                                                                                             >  1:559765/1‑227 (MQ=25)                                       gggCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCgg                                                                                                                                                                                                     >  2:563307/2‑218 (MQ=11)                                         gCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTc                                                                                                                                                                                                                          <  1:377643/195‑1 (MQ=34)                                         gCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGAc                                                                                                                                                                                 >  1:757952/1‑235 (MQ=34)                                         gCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACACggc                                                                                                                                                                                                    >  2:156594/1‑216 (MQ=21)                                         gCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAA‑GCTg                                                                                                                                                                                   <  1:563308/233‑1 (MQ=11)                                             gATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACg                                                                                                                                                                            <  1:432368/236‑1 (MQ=34)                                                      gggCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCa                                                                                                                                                                                                          >  1:111302/1‑198 (MQ=33)                                                      gggCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCaccacc                                                                                                                                                                    >  1:218798/1‑235 (MQ=33)                                                         aaTCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACcgtcg                                                                                                                                                                <  2:374085/235‑1 (MQ=21)                                                            cACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTg                                                                                                                                                                                                                                          <  2:138619/160‑1 (MQ=31)                                                              cggccGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATcaa                                                                                                                                                                                                                <  2:135831/181‑1 (MQ=18)                                                                          aTGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGcac                                                                                                                                                                        >  2:430855/1‑211 (MQ=32)                                                                          aTGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACa                                                                                                                                                 >  2:454537/1‑234 (MQ=255)                                                                                  ggTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACt                                                                                                                                                                                                        <  2:96023/172‑1 (MQ=33)                                                                                       cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACcgtcgtc                                                                                                                                                              >  1:703457/1‑208 (MQ=35)                                                                                        acCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCGTCACCCTGAc                                                                                                                                                                                                                                        >  2:727491/1‑134 (MQ=17)                                                                                              gagaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGCCAGACACCCTCGACGAGGTGGTGGCGAAGCTGGCCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCa                                                                                                                                                        <  1:216869/207‑1 (MQ=21)                                                                                               agaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACgg                                                                                                                                                                           >  2:237654/1‑187 (MQ=33)                                                                                                  aCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAcc                                                                                                                       <  1:323371/236‑1 (MQ=38)                                                                                                     tGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                                                           >  2:177118/1‑150 (MQ=34)                                                                                                     tGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAAcc                                                                                                                                                                                                           <  1:177119/150‑1 (MQ=34)                                                                                                         aCTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTc                                                                                                                                                                                                       >  2:573812/1‑150 (MQ=34)                                                                                                                 gcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                                                                                                                                                                                             >  1:745813/1‑104 (MQ=33)                                                                                                                 gcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                                                                                                                                                                                             <  2:745812/104‑1 (MQ=33)                                                                                                                  cAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCg                                                                                                                      <  2:757951/221‑1 (MQ=38)                                                                                                                          cGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACt                                                                                                                                                                                                        <  1:636021/132‑1 (MQ=25)                                                                                                                                  gTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCg                                                                                                                 <  1:572653/210‑1 (MQ=255)                                                                                                                                          cacCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGAAGGTCTGCCGATCAACAACCCCTCGGAGCTGTACTTCAA‑GCTGACCGa                                                                                                                                                                              >  2:244161/1‑141 (MQ=17)                                                                                                                                             ccTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCa                                                                            >  1:366499/1‑236 (MQ=40)                                                                                                                                               tCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGCCCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGca                                                                                                                                              <  2:374705/169‑1 (MQ=255)                                                                                                                                               tCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCCACCggcc                                                                                                                   >  2:572652/1‑192 (MQ=255)                                                                                                                                                          gtggCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGa                                                                                                               >  1:235139/1‑188 (MQ=40)                                                                                                                                                                     gACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGAc                                                                 >  1:658364/1‑223 (MQ=37)                                                                                                                                                                         gcgACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCa                                                                      <  2:324976/214‑1 (MQ=35)                                                                                                                                                                           gACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCAca                                                                                                                                                                                                        >  1:151032/1‑82 (MQ=25)                                                                                                                                                                                   cGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCg                                                                                                                 <  2:742714/161‑1 (MQ=2)                                                                                                                                                                                          cGAAGGCGGTGAGATCACCTACCCCATCACCCTGACCAACAAAGCCGGTCTGCCGATCAACAACCCCTCGGAGCTGTACTTCAA‑GCTGACCGACGGCCCCACCGTCGTCGTGGCAGCCAACAGCCCCCCCGGTTCGGCCACTGCAACCGcccc                                                                                                                  <  1:212329/153‑1 (MQ=11)                                                                                                                                                                                               gCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTg                          <  1:356002/236‑1 (MQ=34)                                                                                                                                                                                                                  cATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAAACGCCCCGGACAACGACTacgacgcc                                                                                               >  2:229415/1‑141 (MQ=255)                                                                                                                                                                                                                       cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTg                    >  2:39076/1‑218 (MQ=35)                                                                                                                                                                                                                       cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTc     >  2:647023/1‑233 (MQ=37)                                                                                                                                                                                                                       cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAg   >  1:224545/1‑235 (MQ=35)                                                                                                                                                                                                                       cccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAGc  >  1:193890/1‑236 (MQ=35)                                                                                                                                                                                                                        ccTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAAcgaccc                                                                                    >  1:593681/1‑149 (MQ=18)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            ACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAA‑GCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAATCTGGACAAGACCCCGGTCAGC  >  AE015451/202127‑202576

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