Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	195,724	A→G	9.6%	V410V (GTA→GTG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	195,724	0	A	G	9.6%	187.4 / 7.2	61	V410V (GTA→GTG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base G (2/4);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGAC  >  AE015451/195515‑195948                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 tACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                                        <  1:317725/236‑1 (MQ=12)                   gTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                       >  1:627840/1‑235 (MQ=255)                    acacagTCACCATCGAGAACGCCACAGGCGCTAATTTCGAGCAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACATCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCTACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                                       <  1:1068851/229‑3 (MQ=11)                          tCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                                       <  2:582209/228‑3 (MQ=2)                                 cGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCt                                                                                                                                                                                                                       >  1:19047/1‑189 (MQ=25)                                        gCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAc                                                                                                                                                                                                                                >  1:827003/1‑173 (MQ=12)                                            ccGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                            <  1:538606/173‑1 (MQ=255)                                             aggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGCCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                       <  1:349008/208‑1 (MQ=255)                                               gcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                                          <  1:1091451/236‑1 (MQ=255)                                                            gCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGTTCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCACCCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                            >  1:954116/1‑157 (MQ=255)                                                                         aaTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                            >  1:1078198/1‑160 (MQ=38)                                                                             cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACAATTACCGTTGa                                                                                                                                                                                            <  1:712047/172‑1 (MQ=255)                                                                             cGACCCCGGCTCAGACCACCATGAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGTGCCCGTCCGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                                                    <  1:517154/148‑1 (MQ=255)                                                                               aCCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                                       <  2:834135/175‑1 (MQ=255)                                                                                ccccGGCTCAGACCCCGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                                        <  1:439296/157‑1 (MQ=255)                                                                                ccccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                              <  1:246302/151‑1 (MQ=38)                                                                                                               cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTc                                                                                                                                        >  2:501623/1‑190 (MQ=38)                                                                                                               cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCAc                                                                                                           >  1:105553/1‑219 (MQ=39)                                                                                                                               aCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACACCGGCCAATGAc                                                                                                                                           >  1:94591/1‑171 (MQ=38)                                                                                                                                cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                                       <  2:870510/126‑3 (MQ=2)                                                                                                                                cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                                       >  1:870511/1‑124 (MQ=2)                                                                                                                                      cGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGt                                                                                                                                                                                                                     <  2:264333/90‑1 (MQ=255)                                                                                                                                               cgacGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                                        <  2:557762/91‑1 (MQ=255)                                                                                                                                               cGTGGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                            <  2:348698/106‑1 (MQ=255)                                                                                                                                                  cGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                              <  2:652725/101‑1 (MQ=255)                                                                                                                                                   ggTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCGGACGCCGGTAACCGTGAGCCTGTCCAACGGCCa                                                                                                                                                                                                           <  2:225819/87‑1 (MQ=255)                                                                                                                                                    gtgaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                              <  2:48221/144‑1 (MQ=21)                                                                                                                                                           aaGGCGGCGTGATCACTTACACCGTGAGCCTGTGCAAATCTGCCAAGACGCCGGTAACCGTGAAACTGTGCAAGGGCTAAACGATTAGAGATGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGATCGACGGTCAGCGTCa                                                                                                            >  1:769229/1‑174 (MQ=11)                                                                                                                                                              gcggcgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                                                >  1:48444/1‑71 (MQ=255)                                                                                                                                                                  cgTAATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                                                    <  2:563188/63‑1 (MQ=255)                                                                                                                                                                       ccACGTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTggcg                                                                                       <  1:980770/182‑1 (MQ=255)                                                                                                                                                                         aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACg                                 <  1:501624/235‑1 (MQ=255)                                                                                                                                                                               acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTggcgtca                                                                                    <  2:1080577/178‑4 (MQ=255)                                                                                                                                                                                    gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAAc                                                                                                                                                                                                        >  2:1028684/1‑57 (MQ=255)                                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                               >  2:1037477/1‑96 (MQ=255)                                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAcc                                                                                                          >  1:177855/1‑149 (MQ=37)                                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAcc                                                                                                          <  2:177855/149‑1 (MQ=37)                                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCggc                                                                                        >  2:895494/1‑167 (MQ=255)                                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACAATCGAGAGCGCCACTGGCGGCAACTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTc                             >  1:987862/1‑226 (MQ=255)                                                                                                                                                                                      cccTGAGCAACCCTGCCCAGACGCAGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCATTGACGACTACaacaac                                                                                                                               >  1:302452/1‑128 (MQ=255)                                                                                                                                                                                        cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                                                    >  2:913500/1‑41 (MQ=255)                                                                                                                                                                                        cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                                                    <  1:913501/41‑1 (MQ=255)                                                                                                                                                                                        cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCTACGGCAGCGTCGATTTCCAGACCCCGGCCAATTACACCTACAACAGCgc                                                                                                                             >  1:1076058/1‑127 (MQ=11)                                                                                                                                                                                        cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCgg                                                   >  2:1123014/1‑202 (MQ=255)                                                                                                                                                                                        cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACGggcggc                                                                                     >  2:468264/1‑168 (MQ=255)                                                                                                                                                                                         cgAGCAACCCTGCCCAGACGCCGGTAACCGTGAGCCGGTCCAACGGCCAAACGATTTCCGtt                                                                                                                                                                                              <  2:380739/61‑1 (MQ=14)                                                                                                                                                                                            gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                                                      >  1:397553/1‑99 (MQ=255)                                                                                                                                                                                              aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaaca                                                                                                                                 <  2:641975/118‑1 (MQ=12)                                                                                                                                                                                              aaCCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaaca                                                                                                                                 >  1:641976/1‑118 (MQ=255)                                                                                                                                                                                              aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCaa                                                                                                                                                                                                          >  1:13454/1‑45 (MQ=255)                                                                                                                                                                                              aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGGCAGCGTCACCATCGAGAgcgc                                                                                               >  1:343465/1‑152 (MQ=38)                                                                                                                                                                                              aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGAttt                                                                                                                                                              >  1:890804/1‑89 (MQ=255)                                                                                                                                                                                                cccTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                                                         <  2:409169/156‑1 (MQ=40)                                                                                                                                                                                                  aggCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACg                                 <  1:1033401/208‑1 (MQ=255)                                                                                                                                                                                                   tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                              <  1:798411/100‑1 (MQ=255)                                                                                                                                                                                                   tGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCAGGGCAGc                                                                                                                                                                      >  1:273203/1‑76 (MQ=255)                                                                                                                                                                                                        aGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACaccaccacctgtctcttata  >  1:576693/1‑224 (MQ=255)                                                                                                                                                                                                             ccGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACCGGCGGCAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGAcc       >  1:110639/1‑225 (MQ=255)                                                                                                                                                                                                              cGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCg                                                                                                                                                                                        >  2:224730/1‑47 (MQ=255)                                                                                                                                                                                                              cGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCg                            <  2:576692/203‑1 (MQ=255)                                                                                                                                                                                                               ggTAACAGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAg                                                                                                                                                                                          >  1:928432/1‑44 (MQ=255)                                                                                                                                                                                                                  aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGGTCGACGGTCAGCGTCACCATCGAGAACGCCCCCGGCGGCa                                                                                    >  2:326820/1‑143 (MQ=11)                                                                                                                                                                                                                  aaCCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGagag                                                                                                  >  2:81974/1‑129 (MQ=38)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 TACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGAC  >  AE015451/195515‑195948

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