Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	201,664	C→T	8.9%	G2390G (GGC→GGT)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	201,664	0	C	T	8.9%	130.6 / 4.9	45	G2390G (GGC→GGT)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (12/29);  new base T (2/2);  total (14/31)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACC  >  AE015451/201431‑201893                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      cAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGgcggcg                                                                                                                                                                                                                                      <  1:259943/235‑1 (MQ=255)  aGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCgag                                                                                                                                                                                                                                    <  1:85331/236‑1 (MQ=255)  aGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCAACCCCTTCGGGGACCGAAGgcggcg                                                                                                                                                                                                                                      >  1:240188/1‑234 (MQ=255)                 ctggcgAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                                     <  1:544509/231‑1 (MQ=255)                                 gcCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCGTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGCGACCGAAGGCGGCGAGATCACCTACACCATCACGCt                                                                                                                                                                                                                >  1:346279/1‑225 (MQ=255)                                          ccACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCt                                                                                                                                                                                              <  2:408873/234‑1 (MQ=255)                                                cgcgCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGCTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCGACACCATCACGCTGa                                                                                                                                                                                                              >  1:206540/1‑212 (MQ=255)                                                  tgcAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTc                                                                                                                                                                                               <  2:598244/224‑1 (MQ=255)                                                     aaGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACAcc                                                                                                                                                                                                                        >  2:622639/1‑197 (MQ=255)                                                       ggCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                                     <  2:514340/198‑1 (MQ=255)                                                        gCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCAc                                                                                                                                                                              <  1:246099/236‑1 (MQ=255)                                                         ccaTGCCGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTg                                                                                                                                                                                             <  2:531403/218‑1 (MQ=17)                                                         cGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACt                                                                                                                                                                             <  1:238296/236‑1 (MQ=255)                                                           aTGACGTCGATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCg                                                                                                                                                                           <  1:40186/236‑1 (MQ=255)                                                                    aTCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATc                                                                                                                                                                                                                     <  2:705871/185‑1 (MQ=255)                                                                                    gATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGAcc                                                                                                                                                      <  2:680940/232‑1 (MQ=255)                                                                                                              tGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTAGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGCCCAACAAAGACGGGCTGCCGATCAATAAcc                                                                                                                                                                                <  2:540275/180‑1 (MQ=255)                                                                                                                taCGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTg                                                                                                                                                                                             <  2:2026/164‑1 (MQ=255)                                                                                                                 aCGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTac                                                                                                                                                                                                                           >  1:44360/1‑134 (MQ=255)                                                                                                                              caccttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCg                                                                                                                                                                           <  2:565837/165‑1 (MQ=255)                                                                                                                                cGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGgcggcg                                                                                                                                                                                                                                      <  1:711121/108‑1 (MQ=255)                                                                                                                                  ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                                                 >  2:103282/1‑127 (MQ=255)                                                                                                                                  ttCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGc                                                                                                                                                                                                                 <  1:103282/127‑1 (MQ=255)                                                                                                                                       gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTg                                                                                                  <  2:648541/233‑1 (MQ=255)                                                                                                                                             gCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCCTCGTGGCAGCCAACAGCACCACCGGTTCgg                                                                                                        >  1:40415/1‑221 (MQ=255)                                                                                                                                                     gcggcAACGCTTCGGTTCAAGTGACCGACACCCGCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAAt                                                                                                                                                                                    <  2:16091/137‑1 (MQ=255)                                                                                                                                                       ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCaa                                                                                                                                                             >  1:105781/1‑158 (MQ=255)                                                                                                                                                       ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCaa                                                                                                                                                             <  2:105781/158‑1 (MQ=255)                                                                                                                                                                     tCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCa                                                                 <  1:55473/236‑1 (MQ=255)                                                                                                                                                                      cAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTg                                                                                                                                 <  1:88683/171‑1 (MQ=255)                                                                                                                                                                                 cacCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCggt                                                       >  2:638796/1‑234 (MQ=11)                                                                                                                                                                                 cacCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                        <  1:436779/201‑1 (MQ=255)                                                                                                                                                                                 cacCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                        >  2:436778/1‑201 (MQ=255)                                                                                                                                                                                  acCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCggtgg                                                     <  1:597446/235‑1 (MQ=255)                                                                                                                                                                                             aGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGc                                              >  2:472886/1‑231 (MQ=255)                                                                                                                                                                                                 gTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                        >  1:460625/1‑185 (MQ=255)                                                                                                                                                                                                 gTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGAc                                       <  2:309364/234‑1 (MQ=255)                                                                                                                                                                                                 gTGGCGAAGCTGACCGCCACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGcccc                                                                                        <  2:460624/185‑1 (MQ=37)                                                                                                                                                                                                   ggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCaa                                                                                                                         >  2:474796/1‑150 (MQ=255)                                                                                                                                                                                                   ggCGAAGCTGACCGCGACCCCTTCAGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCaa                                                                                                                         <  1:474797/150‑1 (MQ=255)                                                                                                                                                                                                        aGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTg                                                                                                  <  2:549198/168‑1 (MQ=39)                                                                                                                                                                                                           tGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggc                           <  1:311806/236‑1 (MQ=40)                                                                                                                                                                                                              ccGCGACCCCTTCAGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAgcggcg                          <  2:722727/234‑1 (MQ=255)                                                                                                                                                                                                                     cccTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACaa                                                                                  >  2:315129/1‑171 (MQ=34)                                                                                                                                                                                                                          cGGTGACCGAAGGCGGCGAGATCACCTACACCTTCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGAc                                       <  2:397683/209‑1 (MQ=25)                                                                                                                                                                                                                          cGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCATGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCg                                                                                       <  1:552707/161‑1 (MQ=14)                                                                                                                                                                                                                           ggTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAAcc                                                                                                                                                                                <  1:653856/71‑1 (MQ=34)                                                                                                                                                                                                                                    aGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacc  >  1:147031/1‑236 (MQ=18)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      CAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAAGTCGTGGCGAAGCTGACCGCGACCCCTTCGGTGACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACC  >  AE015451/201431‑201893

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