Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	203,422	T→C	9.0%	A2976A (GCT→GCC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	203,422	0	T	C	9.0%	145.2 / 2.9	53	A2976A (GCT→GCC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base C (2/3);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACAC  >  AE015451/203198‑203652                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       aaGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                                             <  1:375020/236‑1 (MQ=33)  aGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                                                                                                                                                            <  1:469824/236‑1 (MQ=33)                  gCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAccc                                                                                                                                                                                                               <  2:696382/233‑1 (MQ=35)                                   tCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                                                                                                                                             >  2:756597/1‑202 (MQ=32)                                          gTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCCAAACCAGCGCGCCGTACACCCACg                                                                                                                                                                                                            >  1:31695/1‑212 (MQ=21)                                                       gaGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGcc                                                                                                                                                                        <  1:276511/235‑1 (MQ=40)                                                        aGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCa                                                                                                                                                                       <  1:66443/235‑1 (MQ=38)                                                                  cccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAcc                                                                                                                                                                                                                <  2:640971/184‑1 (MQ=37)                                                                                 cGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCTCAGGGt                                                                                                                                                                                                 >  1:770266/1‑184 (MQ=11)                                                                                   aCTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCGACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGTGAAACCAGCGCGCCGTAGACCCAcccgt                                                                                                                                                                                                         >  1:760527/1‑170 (MQ=12)                                                                                                 aCCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                        <  2:650083/161‑1 (MQ=37)                                                                                                 aCCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTg                                                                                                                              <  1:755531/235‑1 (MQ=37)                                                                                                     tCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTcg                                                                                                                                                 <  1:365768/212‑1 (MQ=37)                                                                                                     tCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTcg                                                                                                                                                 >  2:365766/1‑212 (MQ=38)                                                                                                      cacaGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCACGGtgct                                                                                                                                                                                              >  2:224406/1‑164 (MQ=37)                                                                                                       acaGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtg                                                                                                                                                                                                <  2:248328/163‑1 (MQ=38)                                                                                                           ccGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTgg                                                                                                                                >  2:536054/1‑223 (MQ=35)                                                                                                                            gCTGAGACCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTgcgc                                                                                                                                                                                                      <  2:147593/136‑1 (MQ=37)                                                                                                                                  aaCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                                                                                                                                              >  1:720653/1‑106 (MQ=33)                                                                                                                                   aCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGAATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGc                                                                                                    >  1:13358/1‑227 (MQ=17)                                                                                                                                          aCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTgg                                                                                    <  1:195969/236‑1 (MQ=35)                                                                                                                                                 ttCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCAc                                                                             <  1:536056/236‑1 (MQ=37)                                                                                                                                                           aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGAGATGACGTTTACAACGACGCTGGCCAGATCAGcc                                                                                                                                                               >  1:137211/1‑144 (MQ=17)                                                                                                                                                             gccAACCCGCCGCTAGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCTCAGGGTGATGACGTTTACAACGACGc                                                                                                                                                                             <  2:137209/126‑1 (MQ=255)                                                                                                                                                             gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTg                                                                                                                              >  1:153489/1‑175 (MQ=31)                                                                                                                                                             gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAGGTGGAt                                                                                  <  2:494015/219‑1 (MQ=18)                                                                                                                                                             gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGTTGGGCGGTGCCGCGAAGGTGGAt                                                                                  >  1:494017/1‑219 (MQ=18)                                                                                                                                                                aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGcc                                                                                                                                                                                                                       <  2:474476/83‑1 (MQ=18)                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAccc                                                                                                                                                                                                               <  2:104892/91‑1 (MQ=33)                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAccc                                                                                                                                                                                                               >  1:104893/1‑91 (MQ=33)                                                                                                                                                                aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAACCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGAGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACATGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcgacga                                                              >  1:761347/1‑236 (MQ=255)                                                                                                                                                                  ccAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGt                                                                                                                               <  2:140979/169‑1 (MQ=31)                                                                                                                                                                     gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtgctg                                                                                                                                                                                             >  2:170738/1‑104 (MQ=17)                                                                                                                                                                        gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGaaa                                                                                                                                                                                                                                 <  1:772409/65‑1 (MQ=35)                                                                                                                                                                        gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                        <  2:747776/90‑1 (MQ=33)                                                                                                                                                                        cttGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTg                                                                                     <  1:105582/203‑1 (MQ=14)                                                                                                                                                                               cACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCa                                                                                                                                                                                                              >  1:501014/1‑77 (MQ=35)                                                                                                                                                                                          cGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                                         >  2:350249/1‑55 (MQ=35)                                                                                                                                                                                                 cTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAgcgcg                                                                                                                                                                                                                         >  1:90146/1‑48 (MQ=38)                                                                                                                                                                                                 cTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGcc                                                                                                                                                                                                                       >  2:756663/1‑50 (MQ=38)                                                                                                                                                                                                 cTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGttc                                                                                                                             >  2:738921/1‑138 (MQ=17)                                                                                                                                                                                                     gCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCAc                          >  1:476909/1‑235 (MQ=37)                                                                                                                                                                                                       aacaaTGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAg                                                 >  2:383638/1‑210 (MQ=38)                                                                                                                                                                                                       aacaaTGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTcggcggg                                                                                                                                            >  2:251156/1‑118 (MQ=17)                                                                                                                                                                                                                 aGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCga              >  1:555552/1‑235 (MQ=35)                                                                                                                                                                                                                 aGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCg                               >  2:273027/1‑218 (MQ=14)                                                                                                                                                                                                                    tCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATc          >  1:572002/1‑236 (MQ=34)                                                                                                                                                                                                                      aCCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc        <  1:273029/236‑1 (MQ=14)                                                                                                                                                                                                                       ccATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc        <  1:488610/235‑1 (MQ=14)                                                                                                                                                                                                                           cAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGGCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggccg                 >  2:431837/1‑220 (MQ=21)                                                                                                                                                                                                                             aGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                <  1:69734/221‑1 (MQ=34)                                                                                                                                                                                                                              ggCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTacac  >  1:250995/1‑234 (MQ=35)                                                                                                                                                                                                                                cTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGGGCCGCGAAGGTGGATGTCAcc                                                                            >  1:135823/1‑158 (MQ=25)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       AAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACAC  >  AE015451/203198‑203652

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