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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	212,794	0	C	T	8.1%	168.6 / 2.6	60	D6100D (GAC→GAT)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (3/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACA  >  AE015451/212572‑213011                                                                                                                                                                                                                               |                                                                                                                                                                                                                          aGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCa                                                                                                                                                                                                              <  1:99494/236‑1 (MQ=32)                    ggCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTcg                                                                                                                                                                                             <  1:725405/234‑1 (MQ=37)                           tGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTcg                                                                                                                                                                                             <  1:499911/227‑1 (MQ=34)                             gTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACg                                                                                                                                                                                   <  1:41587/235‑1 (MQ=33)                                          tCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGAc                                                                                                                                                                      >  1:696655/1‑235 (MQ=33)                                              cgccgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGc                                                                                                                                                                         <  1:285214/228‑1 (MQ=34)                                                 cgACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCgaga                                                                                                                                                              <  1:1106309/236‑1 (MQ=34)                                                       gACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCCGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAg                                                                                                                                                                                                             <  1:700012/183‑1 (MQ=14)                                                                     agaACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCggca                                                                                                                                            >  2:1102287/1‑232 (MQ=255)                                                                                 cGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTg                                                                                                                                                                             <  1:993070/189‑1 (MQ=35)                                                                                       tCACCGTTCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCgaga                                                                                                                                                              <  1:959514/198‑1 (MQ=9)                                                                                         aCCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTcg                                                                                                                                                                                             <  2:193772/165‑1 (MQ=35)                                                                                           cGTGCACTTCACCCCGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCa                                                                                                                                                                                                  <  2:369057/158‑1 (MQ=16)                                                                                                aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGcc                                                                                                                                          <  1:1083472/209‑1 (MQ=35)                                                                                                aCGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAg                                                                                                                                    >  2:966697/1‑215 (MQ=37)                                                                                                             gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACcg                                                                                                              >  2:517949/1‑224 (MQ=35)                                                                                                             gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGACAGAGCAGCCCGGGTATCAACTcggcg                                                                                                                                                                                          >  2:338138/1‑148 (MQ=255)                                                                                                               agggtcACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt                                                                                                                                                                                        <  1:743915/143‑1 (MQ=12)                                                                                                                    cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGtt                                                                                           <  1:517950/236‑1 (MQ=37)                                                                                                                         ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCt                                                                                                                                                                                                                        <  2:1072371/106‑1 (MQ=33)                                                                                                                             gTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGGATCAACTCGGCGGTGGACGCCACTGGGGCGACCTTCGAGAACCTGGAg                                                                                                                                                      >  2:314707/1‑168 (MQ=18)                                                                                                                              tCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTc                                                                                                      <  1:515057/215‑1 (MQ=37)                                                                                                                               cgggACCCTGAGGAACAACGCCCAGTTAACCATCAAGGCGGGCCAAAGCAGTGTGTCATACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCt                                                                                                                                                                                                                        <  1:598813/97‑1 (MQ=11)                                                                                                                                      cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGt                                                                                                                                             >  1:1111856/1‑168 (MQ=35)                                                                                                                                      cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCAACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                                                              <  1:940757/215‑1 (MQ=37)                                                                                                                                       tGAGCAACAACGTCCAGGTCACCATCAAGGCTGGCGAACCCAGCGCGCCGTACGCCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGAGACCTGCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGAc                                                                                        <  1:145524/220‑1 (MQ=255)                                                                                                                                           caacaaCGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAg                                                                                             >  1:30873/1‑211 (MQ=35)                                                                                                                                                   cccAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCt                                                                                                                                                    <  1:258098/148‑1 (MQ=37)                                                                                                                                                   cccAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGGAGGTGGac                                                                                                                              >  2:940756/1‑169 (MQ=2)                                                                                                                                                     cAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc                                                            <  2:696654/234‑1 (MQ=34)                                                                                                                                                     cAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcg                                                              <  2:699042/232‑1 (MQ=18)                                                                                                                                                      aGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCga                                                          >  1:143102/1‑235 (MQ=33)                                                                                                                                                      aGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCg                                                                    >  2:994368/1‑225 (MQ=16)                                                                                                                                                        gTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCGCGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAg                                                                 <  1:700914/226‑1 (MQ=21)                                                                                                                                                          cACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGGGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAgg                                                                <  1:1089604/225‑1 (MQ=21)                                                                                                                                                           aCCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTAc                                                    <  1:358291/236‑1 (MQ=35)                                                                                                                                                                cAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGTCGCCACTGGTGc                                                                                                                                                                         >  2:258097/1‑114 (MQ=255)                                                                                                                                                                     cTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGCACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGACGGCGGCGAGATTACCTACACCATc                                          >  2:104950/1‑236 (MQ=37)                                                                                                                                                                     cTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGGGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACtt                                                                                                                                                                                              >  1:303374/1‑87 (MQ=255)                                                                                                                                                                             aCCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGAc                                  <  1:994369/236‑1 (MQ=25)                                                                                                                                                                                    cgacGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGc                                                                                                                                                <  1:68496/116‑1 (MQ=25)                                                                                                                                                                                       cGTACACCCACGCGGCGCAGGGTGATGACGGCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAg                          <  1:759851/234‑1 (MQ=25)                                                                                                                                                                                          acacCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGAcacc                                                                                                                  >  2:695762/1‑143 (MQ=37)                                                                                                                                                                                           cacCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGcc                    <  1:942492/236‑1 (MQ=37)                                                                                                                                                                                              ccACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                 <  1:482589/236‑1 (MQ=39)                                                                                                                                                                                              ccACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGc                 <  1:22242/236‑1 (MQ=37)                                                                                                                                                                                                 cGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTcgg                                                                                                                                                                                            >  1:617106/1‑62 (MQ=35)                                                                                                                                                                                                    tGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGCGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGGGCCGCGAAGGGGGATGTCAc                                                                                                                         >  2:342390/1‑126 (MQ=255)                                                                                                                                                                                                        cAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGcc                >  1:15286/1‑227 (MQ=37)                                                                                                                                                                                                         aaggcgatgaCGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCTGCGGTGGACGTCGATGGTCCGACCTTCGAGAACCTGGATCTGGGCGGCGCCGCTAAGGTTGAAGTCACCGACACCACCGACGAAGTCGTGGcc                                                                                                >  2:1070113/6‑146 (MQ=11)                                                                                                                                                                                                         aGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTg                                                                                                                                                   >  2:478310/1‑95 (MQ=37)                                                                                                                                                                                                         aGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGTTGTCACCGACACCACcgaa                                                                                                            >  1:170336/1‑133 (MQ=21)                                                                                                                                                                                                            gTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAAcc     <  2:871970/234‑1 (MQ=40)                                                                                                                                                                                                              gatgaCGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATc                                          <  1:737941/195‑1 (MQ=37)                                                                                                                                                                                                              gatgaCGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATc                                          >  2:737940/1‑195 (MQ=37)                                                                                                                                                                                                              gatgaCGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACcacc  >  1:1104819/1‑234 (MQ=39)                                                                                                                                                                                                               atgaCGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGCGCCGCGTAgg                                                                                                                                   >  2:184547/1‑105 (MQ=255)                                                                                                                                                                                                                        aCAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACcaaca                             >  2:907135/1‑198 (MQ=34)                                                                                                                                                                                                                        aCAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAACTGGGCGGTACCGCGAAGGTGCAAGTGACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTa                                                 >  1:808846/1‑178 (MQ=11)                                                                                                                                                                                                                          aaCGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTCGCCAAGTTGACCGCCAcccc                                                                               >  1:1072779/1‑146 (MQ=21)                                                                                                                                                                                                                            cgatgCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGgcggcg                                                           <  1:451524/164‑1 (MQ=21)                                                                                                                                                                                                                               |                                                                                                                                                                                                                          AGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACA  >  AE015451/212572‑213011

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