Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	213,619	T→C	8.7%	I6375I (ATT→ATC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	213,619	0	T	C	8.7%	177.9 / 4.3	57	I6375I (ATT→ATC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base C (3/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTG  >  AE015451/213396‑213847                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     cTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATCACCTACACCAt                                                                                                                                                                                                                           >  1:651563/1‑235 (MQ=255)         aCCAGCGCGCCAGACACCCACGCTGCGCAAGGCGATGACGTCTACAGCGATGCTGGCCAGATCAGCCTGGGCTTTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATc                                                                                                                                                                                                                          <  1:104708/228‑1 (MQ=255)                   cATACACCCACGCAGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGAc                                                                                                                                                                                                        <  1:467715/236‑1 (MQ=255)                         cccACGCTGCGCAGGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                                                    <  1:81949/202‑1 (MQ=255)                                      aGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCaacaac                                                                                                                                                                                      <  1:260516/235‑1 (MQ=255)                                                gTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGc                                                                                                                                                                              <  2:403455/233‑1 (MQ=255)                                                         ggtGCTGGCAAGATCAGCCTGGGCATTAACTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                                                    <  1:208294/168‑1 (MQ=255)                                                            gCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCAc                                                                                                                                                                                                                        <  2:52901/179‑1 (MQ=255)                                                             cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGATGGTCTGCCTATcgacaagcatgcggc                                                                                                                                                                              <  2:660930/220‑10 (MQ=11)                                                             cTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGATGGTCTGCCTATcgacaagcatgcggc                                                                                                                                                                              >  1:660931/1‑211 (MQ=11)                                                              tGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGTTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACAcct                                                                                                                                                                                                                            >  1:384703/1‑172 (MQ=255)                                                              tGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCTTCGGTAACCGAAGGCGGCGAGATTAc                                                                                                                                                                                                                                    >  1:308446/1‑165 (MQ=255)                                                                  cAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGa                                                                                                                                                         >  1:320422/1‑236 (MQ=255)                                                                        aGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGAc                                                                                                                                                                                                        <  2:129927/183‑1 (MQ=255)                                                                        aGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACCGTGCATTGACCTTCACGCt                                                                                                                                                               >  2:681507/1‑224 (MQ=255)                                                                               gCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACATACACCATCACGCTg                                                                                                                                                                                                                    >  2:479091/1‑164 (MQ=255)                                                                                       tCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGTCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                <  2:462922/224‑1 (MQ=255)                                                                                       tCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAAc                                                                                                                                                >  1:462923/1‑224 (MQ=255)                                                                                            aGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGAccc                                                                                                                                                                      >  2:135904/1‑196 (MQ=255)                                                                                            aGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACAGCCTGCCGATCAACAACCACAGTGCATTGACCTTCa                                                                                                                                                                   >  2:669106/1‑200 (MQ=255)                                                                                                 aCGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCt                                                                                                                                                                                                   <  2:22415/163‑1 (MQ=255)                                                                                                 aCGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                                                                                                             >  1:178318/1‑217 (MQ=40)                                                                                                 aCGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAAc                                                                                                                          <  1:632367/236‑1 (MQ=255)                                                                                                      aCTGGTGCGATCTTCGAGAATCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATcaacaa                                                                                                                                                                                       <  2:692615/170‑1 (MQ=255)                                                                                                                  ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGCCGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGCAGGCGGCGAGATTACCTACACCATCACGCTGACCAACACAGACGGCCTGCCGACCAACAACcacact                                                                                                                                                                                >  2:11230/1‑163 (MQ=255)                                                                                                                  ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                 >  2:565595/1‑148 (MQ=255)                                                                                                                  tgcGAGAACCTGCAACTGGGCGGCAACGCTTGGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGAGAGATGACCTACACCATCACGCTGACCAACAAAGACGGGCTGCCGATCAACAGCCACAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                     <  2:9628/190‑1 (MQ=255)                                                                                                                     gagaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAcc                                                                                                                                      <  1:686922/204‑1 (MQ=38)                                                                                                                         aCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                    <  1:687930/234‑1 (MQ=255)                                                                                                                                       ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACg                                                                                                                                                                                                       >  1:559977/1‑121 (MQ=255)                                                                                                                                       ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACg                                                                                                                                                                                                       <  2:559976/121‑1 (MQ=255)                                                                                                                                       ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACa                                                                                     >  1:436108/1‑235 (MQ=255)                                                                                                                                        gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                     >  2:720414/1‑138 (MQ=255)                                                                                                                                        gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAAcc                                                                                                                                                                                     <  1:720415/138‑1 (MQ=255)                                                                                                                                            cGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCaa                                                                                                                                                   >  1:346329/1‑168 (MQ=255)                                                                                                                                                          tGACCGGCACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                    <  1:405605/201‑1 (MQ=38)                                                                                                                                                          ggACCGACACCGCCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCGCAGTGCATTGACCTTCACGCTGAGCGACGGc                                                                                                                                                     <  2:182640/151‑1 (MQ=255)                                                                                                                                                             ccGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                    <  1:731757/198‑1 (MQ=37)                                                                                                                                                                   ccaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                    >  2:661395/1‑192 (MQ=37)                                                                                                                                                                   ccaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCg                                                                                                    <  1:661396/192‑1 (MQ=37)                                                                                                                                                                   ccaccGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAAc                                                                                   <  2:54997/209‑1 (MQ=40)                                                                                                                                                                   ccaccGACGAAGTAGTGGCCAAGGTGACCGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACg                                                                <  1:697419/228‑1 (MQ=11)                                                                                                                                                                                          ttGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGa                                 >  1:64844/1‑236 (MQ=40)                                                                                                                                                                                            gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCgg                                                                                        <  1:538125/179‑1 (MQ=255)                                                                                                                                                                                            gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCgg                                                                                        >  2:538124/1‑179 (MQ=40)                                                                                                                                                                                            gACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGACTGCCGATcgacaagc                                                                                                                                                                                     >  2:167238/1‑84 (MQ=255)                                                                                                                                                                                                    ccccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCACCGGCACCGTGGGTACTGCCCCCGTGACTGCTCCGGACac                                                                                    >  2:32975/1‑174 (MQ=255)                                                                                                                                                                                                      ccGTCGGTGACCGAAGGCGGCGAGATTACCTACGCCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCAGTGACCTGCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTc                                                                          <  2:15435/183‑1 (MQ=255)                                                                                                                                                                                                      ccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                 <  2:712854/64‑1 (MQ=35)                                                                                                                                                                                                      ccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                                                                                                                                                                                                 >  1:712855/1‑64 (MQ=35)                                                                                                                                                                                                      ccGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATc                                                                                                                                                                                            <  2:57382/69‑1 (MQ=35)                                                                                                                                                                                                        ggcGGTGACCGAAGGCGGCGAGATTACCTACACCAGCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACg                                                                <  1:3073/189‑1 (MQ=255)                                                                                                                                                                                                        gTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCACCAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGt                                                                       <  2:127673/184‑1 (MQ=39)                                                                                                                                                                                                          cGGTGACCGAAGGCGGCGAGATTACCGACACCATCCCGCTGACCAGCAAAGACGGCCTGCCGGTCACCAACCACAGTGCATTGACTGTCACGCTGAGCGACgg                                                                                                                                                      >  2:346328/1‑103 (MQ=11)                                                                                                                                                                                                                         gcgAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTg  <  1:746550/236‑1 (MQ=40)                                                                                                                                                                                                                            agaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCATGGGTACTGCCACCGTGACTGCCCCGGAcct                                                                                    >  2:340104/1‑149 (MQ=9)                                                                                                                                                                                                                            agaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCg          >  1:191095/1‑225 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     CTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGTTGACCGCGACCCCGTCGGTGACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAGGTGCGGACGTTGGCAAGTTCGAGCAACTG  >  AE015451/213396‑213847

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