Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	212,134	T→C	10.1%	D5880D (GAT→GAC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	212,134	0	T	C	10.1%	153.8 / 5.9	56	D5880D (GAT→GAC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base C (2/4);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACC  >  AE015451/211907‑212348                                                                                                                                                                                                                                     |                                                                                                                                                                                                                       gaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCa                                                                                                                                                                                                                <  1:565383/236‑1 (MQ=40)             ccGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCa                                                                                                                                                                                                           <  2:13986/229‑1 (MQ=34)                                cggTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCAAGATCAGCCTGGGATTCAACTcggcggag                                                                                                                                                                                    >  2:270627/2‑231 (MQ=11)                                                          cAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTg                                                                                                                                                                                                      >  2:180139/1‑189 (MQ=33)                                                          cAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTC‑‑GACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGt                                                                                                                                                                                                   >  2:371046/1‑190 (MQ=18)                                                                           gaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCgg                                                                                                                                      <  1:866118/236‑1 (MQ=34)                                                                                                          aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTg                                                                                                                                                                       <  2:300730/172‑1 (MQ=18)                                                                                                          aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                          >  2:938002/1‑169 (MQ=34)                                                                                                          aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                          <  1:938005/169‑1 (MQ=35)                                                                                                            ccAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAag                                                                                                     <  1:729393/236‑1 (MQ=255)                                                                                                               gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCt                                                                                                                                                                                                       >  1:819283/1‑135 (MQ=34)                                                                                                                  gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTggg                                                                                                                                                                                                    <  1:327934/135‑1 (MQ=21)                                                                                                                  gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTggg                                                                                                                                                                                                    >  2:327931/1‑135 (MQ=25)                                                                                                                   cagggtgaCGGCCTGGTTGGGACCCTGCGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCtt                                                                                                                                                                <  1:759403/163‑1 (MQ=11)                                                                                                                       cTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCg                                                                                                                                            >  1:348877/1‑186 (MQ=35)                                                                                                                                     gTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAag                                                                                                     >  2:362651/1‑211 (MQ=255)                                                                                                                                         cccTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCg                                                                        >  1:922327/1‑236 (MQ=255)                                                                                                                                           cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc                                                                                                                                                                               >  1:51921/1‑131 (MQ=34)                                                                                                                                           cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGt                                                                      >  1:30218/1‑236 (MQ=255)                                                                                                                                               gCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCAGTCGGTGAc                                                                   >  1:320684/1‑235 (MQ=255)                                                                                                                                                caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTCGTGGCCAAGCTGACTGCGAcc                                                                              >  1:475745/1‑223 (MQ=255)                                                                                                                                                caacaaTGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGAcc                                                                  >  1:617701/1‑235 (MQ=255)                                                                                                                                                      tGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGccg                                                                                                                                      >  1:247919/1‑159 (MQ=21)                                                                                                                                                        cccAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                         <  1:180141/236‑1 (MQ=255)                                                                                                                                                              tCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATc                                                   <  1:270630/236‑1 (MQ=255)                                                                                                                                                                aCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                 <  1:638276/236‑1 (MQ=255)                                                                                                                                                                aCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                 <  1:486791/236‑1 (MQ=255)                                                                                                                                                                    tCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTa                                              >  1:806948/1‑235 (MQ=255)                                                                                                                                                                     cAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGGTCCCCTGGATGCATCTC‑CGTACCCAGGTGCCTGCTAGATCAGCCTGAGTATGAACTCGGCGGGGGAAGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCg                                                                        >  1:544336/1‑208 (MQ=14)                                                                                                                                                                     cAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGcc                                                                                             >  1:325093/1‑187 (MQ=255)                                                                                                                                                                      aaGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                         <  1:56306/222‑1 (MQ=255)                                                                                                                                                                         gCCGGTGAAACCAGCGCGCCGTACACCCACGCGGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCa                                                  <  2:249612/226‑1 (MQ=255)                                                                                                                                                                               gAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCt                                  >  1:447850/1‑236 (MQ=255)                                                                                                                                                                                aaaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCg                                                                                                                                            >  1:493834/1‑129 (MQ=35)                                                                                                                                                                                 aaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGa                                >  1:689795/1‑236 (MQ=255)                                                                                                                                                                                 aaCCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCAAGATCAGCCTGGGTATCAACTcggcgga                                                                                                                                                                                     >  2:290944/1‑86 (MQ=18)                                                                                                                                                                                  aCCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCAc                                                 <  2:617698/218‑1 (MQ=255)                                                                                                                                                                                        gcgcCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGt                                                         <  1:793627/204‑1 (MQ=255)                                                                                                                                                                                                 acCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGGGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACCGCGACCCCGTCGGTGACCGAAGGCGGCga                                                       >  1:266747/1‑197 (MQ=255)                                                                                                                                                                                                     aCGCTGCGCAGGGTG‑ATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCg            <  1:540590/236‑1 (MQ=255)                                                                                                                                                                                                        cTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACcaa                            <  1:84603/217‑1 (MQ=255)                                                                                                                                                                                                           cgcAGTTTG‑ATGACGTCTACGACGGTGCGGGTCAGATCAGCCTGGGTCTCAACTCGGCGGTGGACGCCACTGGTGCGACCTTc                                                                                                                                                               <  1:318293/83‑1 (MQ=11)                                                                                                                                                                                                           cgcAGGGTG‑ATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTc                                                                                                                              <  1:326416/116‑1 (MQ=18)                                                                                                                                                                                                              aGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGt                                                                                                                                                                        >  1:112427/1‑71 (MQ=34)                                                                                                                                                                                                               gggTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAAcc  <  1:362654/236‑1 (MQ=255)                                                                                                                                                                                                               gggTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAAcc  <  1:793400/236‑1 (MQ=255)                                                                                                                                                                                                                ggTG‑ATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCGTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAAcc  <  1:9795/235‑1 (MQ=255)                                                                                                                                                                                                                ggTG‑ATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAAcc  <  1:384586/235‑1 (MQ=255)                                                                                                                                                                                                                ggTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAAcc  >  1:875688/1‑235 (MQ=255)                                                                                                                                                                                                                              aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCg            >  1:668468/1‑212 (MQ=255)                                                                                                                                                                                                                              aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCg            <  2:668465/212‑1 (MQ=255)                                                                                                                                                                                                                              aCAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACAACACCGACGAAGTAGTGGCCa                                                                                            >  1:365445/1‑132 (MQ=255)                                                                                                                                                                                                                              aCAACGACGCTGGCCAGATAAGCCTGGGTATCAACTCAGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGTTGGGCGGTGCCGCTTCGGTCCAAGTGACCGACACCCTCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCgg                                                                       >  2:882635/1‑153 (MQ=11)                                                                                                                                                                                                                                aaCGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGAc                               >  1:660517/1‑191 (MQ=255)                                                                                                                                                                                                                                aaCGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCGTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACAcc                                          >  2:814780/1‑180 (MQ=255)                                                                                                                                                                                                                                aaCGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGCCGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACcaacaa                         <  1:200085/197‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                       GAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTG‑ATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTCCAAGTGACCGACACCACCGACGAAGTAGTGGCCAAGCTGACTGCGACCCCGTCGGTGACCGAAGGCGGTGAGATCACTTACACCATTACCCTGACCAACAAAGACGGTCTGCCGATTAACAACC  >  AE015451/211907‑212348

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