Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	196,027	T→C	11.1%	T511T (ACT→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,027	0	T	C	11.1%	131.5 / 3.7	55	T511T (ACT→ACC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (24/21);  new base C (3/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

CAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCG  >  AE015451/195820‑196238                                                                                                                                                                                                                |                                                                                                                                                                                                                    caacGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                  >  1:661839/1‑227 (MQ=38)       cTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACAGTGACCCTGTCCAACGGCCaa                                                                                                                                                                                                >  1:479859/1‑223 (MQ=25)                         cATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCaa                                                                                                                                                                                                >  1:863939/1‑205 (MQ=38)                             gagaGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                  <  1:893943/215‑1 (MQ=255)                              agacagcCACAGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                             <  1:809370/214‑1 (MQ=255)                                     cACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCaaa                                                                                                                                                                                               >  1:505454/1‑194 (MQ=37)                                                         aacTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                      >  1:907695/3‑167 (MQ=2)                                                         aacTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                      <  2:907692/165‑1 (MQ=9)                                                              aCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATg                                                                                                                                   >  1:891622/1‑229 (MQ=40)                                                               ccccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                              <  1:950961/169‑1 (MQ=38)                                                                 ccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                <  2:581383/213‑1 (MQ=37)                                                                          accccGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                    <  1:437687/165‑1 (MQ=25)                                                                            gCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                      <  2:891619/212‑1 (MQ=39)                                                                                         aCGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                        >  2:903653/1‑229 (MQ=255)                                                                                                          cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACgg                                                                                >  1:485451/1‑236 (MQ=21)                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAAc                                                              >  1:961951/1‑235 (MQ=255)                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                <  2:609954/153‑1 (MQ=38)                                                                                                                                  ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACtgt                                                                                                                                                                                                                   <  2:144631/81‑1 (MQ=38)                                                                                                                                      tgCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACg                                                                                      <  2:485448/201‑1 (MQ=38)                                                                                                                                             cGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGaa                                                    >  1:203924/1‑229 (MQ=255)                                                                                                                                               gtgaCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                    <  1:51765/96‑1 (MQ=255)                                                                                                                                               gtgaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcccc                                                                                                                                           <  2:718575/137‑1 (MQ=25)                                                                                                                                                tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                >  1:34266/1‑134 (MQ=39)                                                                                                                                                tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                <  2:34265/134‑1 (MQ=39)                                                                                                                                                           ggcgTTATCACCTACACCGTGACCCTTAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                                      <  1:761620/69‑1 (MQ=255)                                                                                                                                                             cgTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAAcccca                                               >  2:352744/1‑216 (MQ=255)                                                                                                                                                                 aTCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                             <  2:829740/88‑1 (MQ=255)                                                                                                                                                                 aTCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                             >  1:829743/1‑88 (MQ=18)                                                                                                                                                                        acgccGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGc                                        <  1:522526/211‑1 (MQ=255)                                                                                                                                                                           ccGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGt                <  1:86682/235‑1 (MQ=255)                                                                                                                                                                            cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACg                                             >  1:703335/1‑205 (MQ=255)                                                                                                                                                                              tGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACgg                                                                                                                                                                                                    >  2:846155/1‑52 (MQ=35)                                                                                                                                                                               gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGAc                                                                                                                                                                                                                >  1:36218/1‑39 (MQ=40)                                                                                                                                                                               gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                   <  2:307971/68‑1 (MQ=255)                                                                                                                                                                               gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGcc                                       >  1:248758/1‑207 (MQ=255)                                                                                                                                                                                 cccTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTc                   >  1:60655/1‑226 (MQ=255)                                                                                                                                                                                   cTGAGCAATCCTGCCCAGACGCCGGTAACt                                                                                                                                                                                                                     >  2:1057119/1‑30 (MQ=255)                                                                                                                                                                                   cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                                      >  2:1043294/1‑45 (MQ=38)                                                                                                                                                                                   cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                                      <  1:1043297/45‑1 (MQ=38)                                                                                                                                                                                     gAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGCCACCATCGGGAGCGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCg                                          <  2:1025087/199‑1 (MQ=255)                                                                                                                                                                                       gCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGAATACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGCCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCg                                                                             >  1:893045/1‑162 (MQ=255)                                                                                                                                                                                         aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAAc                                                                                                                                                                                                      >  2:257223/1‑39 (MQ=40)                                                                                                                                                                                         aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                             <  1:48076/64‑1 (MQ=255)                                                                                                                                                                                         aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCCCAGACCACGATCAACGACTCGGTCGACACCCCTACcg  >  1:328757/1‑235 (MQ=255)                                                                                                                                                                                         aaTCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGTAAAACCCAGGCCAGCGACAATTTCCACACCCCGGCGAATACCGTCCACAACAACAGCTCGACGGCCAGCGTCAGCATCGAAAACGCCACGGGCGATAActt                                                                      >  1:411513/1‑167 (MQ=14)                                                                                                                                                                                          acccTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCg                                          <  2:961948/192‑1 (MQ=255)                                                                                                                                                                                           ccccgCCCAGACGCCGGTAACCGTGACCCTGGCCAACGGCCAAACGATTTCCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAAttt                                                                      <  2:300062/161‑1 (MQ=255)                                                                                                                                                                                              tGCCCAGACGCCGGTAACTGTGACCCTGTCCAACgg                                                                                                                                                                                                    >  1:1079244/1‑36 (MQ=40)                                                                                                                                                                                              tGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCTGGGAAAACCCAGGGCATGGTCGTTATACAGACCCCGGCTACTGACAAAGACAACAACCGTTCGACGGCCAGCGTCAGCATCGAGAAAGCCACAGGCGGTa                                                                          >  1:888250/1‑158 (MQ=14)                                                                                                                                                                                                   aGACGCCGGTAACTGTGACCCTGTCCAACTGCCAAACGAATGCCGTTGACGCGGGTAAAACCCAGGACAGCGTCGGTTTCCGGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAAttt                                                                      >  2:990750/1‑157 (MQ=14)                                                                                                                                                                                                        ccGGTGACAGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                    <  2:463829/90‑1 (MQ=255)                                                                                                                                                                                                        ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGt                                                                                                                                                                                     >  1:61571/1‑41 (MQ=255)                                                                                                                                                                                                        ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCaag                              >  1:5235/1‑190 (MQ=255)                                                                                                                                                                                                          ggTGACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTaaa                                                                                                                                                                        <  2:676091/52‑1 (MQ=255)                                                                                                                                                                                                           gTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                  <  1:859457/41‑1 (MQ=255)                                                                                                                                                                                                             aaCTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                           <  2:616088/62‑1 (MQ=40)                                                                                                                                                                                                             aaCTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                                                                                           >  1:616091/1‑62 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                    CAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCG  >  AE015451/195820‑196238

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