Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	206,662	G→A	11.0%	T4056T (ACG→ACA)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	206,662	0	G	A	11.0%	141.3 / 8.5	58	T4056T (ACG→ACA)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (20/31);  new base A (2/4);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTT  >  AE015451/206460‑206895                                                                                                                                                                                                           |                                                                                                                                                                                                                                          tGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCa                                                                                                                                                                                                            >  1:480922/1‑234 (MQ=35)      gCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGAc                                                                                                                                                                                                     <  1:112038/236‑1 (MQ=33)        cccGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACtt                                                                                                                                                                                                                              >  1:517119/1‑209 (MQ=32)                    cTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGt                                                                                                                                                                                                   <  1:228408/224‑1 (MQ=34)                                cAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGc                                                                                                                                                                              >  2:266413/1‑233 (MQ=34)                                    gACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGc                                                                                                                                                                       >  1:651072/1‑236 (MQ=34)                                      cccTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAAcc                                                                                                                                                                                                        <  1:121423/201‑1 (MQ=34)                                         tGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACg                                                                                                                                                                  <  1:566846/236‑1 (MQ=34)                                                aTGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                            <  1:266414/235‑1 (MQ=34)                                                 tGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCgg                                                                                                                                                                                                                           <  2:427877/171‑1 (MQ=34)                                                   tCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGCAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACg                                                                                                                                                                                                               >  2:121423/1‑181 (MQ=25)                                                   tCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGCGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCg                                                                                                                                                                            >  1:178037/1‑216 (MQ=21)                                                           aaCCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCgaga                                                                                                                                                                                                                  <  2:582561/170‑1 (MQ=35)                                                                   gAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCCTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                                        <  1:411606/236‑1 (MQ=18)                                                                     ggTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGc                                                                                                                                      >  1:184638/1‑236 (MQ=33)                                                                          ggACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                                  >  1:618769/1‑235 (MQ=35)                                                                            aCGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                               >  2:283144/1‑220 (MQ=34)                                                                                    aaGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                                        <  2:618768/235‑1 (MQ=33)                                                                                          gAGCAACGGACGCTGGACAAGCCCCCGGTCAGCACTACCGTCGCCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGcc                                                                                                                                                                             <  1:150378/176‑1 (MQ=16)                                                                                                            aaGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAAcc                                                                                                                                                                                                        <  2:232815/131‑1 (MQ=34)                                                                                                                 cccGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                                                                                                                             <  1:662506/121‑1 (MQ=35)                                                                                                                 cccGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTg                                                                                                                                                 >  2:534494/1‑181 (MQ=33)                                                                                                                 cccGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTGCGGTGGCCGAGAACGTCaaa                                                                                                                                                                                                          >  2:633529/1‑124 (MQ=25)                                                                                                                      tCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                                     <  1:266710/236‑1 (MQ=33)                                                                                                                       cAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                                     <  2:364700/235‑1 (MQ=32)                                                                                                                         gCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGa                                                                                                                                                                 <  2:536064/157‑1 (MQ=34)                                                                                                                                     ccGACGCGCCTGTCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTg                                                                                                                                                                         <  1:534702/137‑1 (MQ=21)                                                                                                                                     ccGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAAAGTCAAACCGACGTTCACCGTACACg                                                                                                                                                                                       >  2:410574/1‑123 (MQ=37)                                                                                                                                       gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCTCGCCGTACACCCACGCTGCGCAGGGTGATGACGt                                                                    >  1:171354/1‑236 (MQ=17)                                                                                                                                          gaGCCTGGCACCCCAGGCCATCCGGGCGGCCGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                                                                                        <  1:248041/181‑1 (MQ=18)                                                                                                                                           aGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTa                                                                 <  1:437826/235‑1 (MQ=32)                                                                                                                                            gCCTGGCACACCAGGCAATCCGGGCGGCAGAAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCCCCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAgc                                                                                                       <  1:181319/196‑1 (MQ=18)                                                                                                                                               tGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTgg                                                                                                                 <  2:559629/183‑1 (MQ=33)                                                                                                                                                  cACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTgcg                                                                                  >  1:162511/1‑211 (MQ=32)                                                                                                                                                    ccccAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGCCGTGACCCTGAGCAACAACGCCCAgg                                                                                                                                 >  2:117492/1‑162 (MQ=18)                                                                                                                                                        aGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                   <  1:534495/236‑1 (MQ=255)                                                                                                                                                          gCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCAc                                                                                                                                                                   >  1:348419/1‑122 (MQ=37)                                                                                                                                                          gCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCAc                                                                                                                                                                   <  2:348418/122‑1 (MQ=37)                                                                                                                                                          gCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGcc                                                                                                                   >  2:724200/1‑169 (MQ=25)                                                                                                                                                                     ggcAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTCTCCGAGAACGTCAATCCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                                                                                       <  2:392761/123‑1 (MQ=12)                                                                                                                                                                         gcaACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCa                                                                                                                                   <  1:480127/139‑1 (MQ=17)                                                                                                                                                                          caACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                                  <  1:480880/139‑1 (MQ=32)                                                                                                                                                                          caACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTAAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAgc                                                                                                       <  2:34595/166‑1 (MQ=14)                                                                                                                                                                           aaCGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCa                                 <  1:298633/235‑1 (MQ=31)                                                                                                                                                                           aaCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGc                                                                                                                                                                              <  2:594568/94‑1 (MQ=21)                                                                                                                                                                           aaCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGCCAAACCGACGTTCACCGTACACGTCAACCAGc                                                                                                                                                                              >  1:594569/1‑94 (MQ=255)                                                                                                                                                                            aCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCa                                                                                                         >  2:285406/1‑162 (MQ=17)                                                                                                                                                                            aCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCa                                                                                                         <  1:285407/162‑1 (MQ=14)                                                                                                                                                                                     aCCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACTGTACACGTCAACCAGc                                                                                                                                                                              <  2:581849/84‑1 (MQ=18)                                                                                                                                                                                     aCCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                       <  1:340681/235‑1 (MQ=9)                                                                                                                                                                                       cTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGa                    <  1:6262/236‑1 (MQ=255)                                                                                                                                                                                         gggcAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGAc                   <  1:85494/232‑1 (MQ=9)                                                                                                                                                                                                  cACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTg         >  1:648488/1‑236 (MQ=31)                                                                                                                                                                                                     cATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                       <  2:480921/219‑1 (MQ=31)                                                                                                                                                                                                       tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCa                                                                                                                                                                                                >  1:565240/1‑48 (MQ=40)                                                                                                                                                                                                       tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCa                                                                                                                                                                                                <  2:565239/48‑1 (MQ=40)                                                                                                                                                                                                       tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGc        >  2:715362/1‑232 (MQ=31)                                                                                                                                                                                                         aCGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCtt  <  1:628276/236‑1 (MQ=14)                                                                                                                                                                                                           |                                                                                                                                                                                                                                          TGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTT  >  AE015451/206460‑206895

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: ‑