Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,312,061	C→T	8.9%	R78C (CGC→TGC)	aaeA_2 →	p‑hydroxybenzoic acid efflux pump subunit AaeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,312,061	0	C	T	8.9%	177.6 / 4.5	56	R78C (CGC→TGC)	aaeA_2	p‑hydroxybenzoic acid efflux pump subunit AaeA
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base T (3/2);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACG  >  ACB143/1311926‑1312196                                                                                                                                        |                                                                                                                                        tGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAgc                                                                                                                                         >  2:34636/1‑136 (MQ=255) tGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACGTGTTGATGAGTATTGACCGCGAgc                                                                                                                                         >  2:4621/1‑136 (MQ=255)    aCGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCt                                                                                                                                      >  8:226603/1‑136 (MQ=255)        gTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCa                                                                                                                                  <  7:148476/136‑1 (MQ=255)          gATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGCCTGGTTGATGAGTATTGACAGCGAGCGCTTCCAgg                                                                                                                                >  6:307206/1‑136 (MQ=255)          gATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGTGCTTCCAgg                                                                                                                                >  7:308350/1‑136 (MQ=255)          gATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGTGCTTCCAgg                                                                                                                                >  3:177768/1‑136 (MQ=255)          gATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAgg                                                                                                                                >  7:14309/1‑136 (MQ=255)          gATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAgg                                                                                                                                >  8:124795/1‑136 (MQ=255)              cccGTGTGTGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGTCcg                                                                                                                             <  8:3444/135‑1 (MQ=255)              cccGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCcgc                                                                                                                            >  6:206357/1‑136 (MQ=255)                 gtgtGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCt                                                                                                                         >  5:355840/1‑136 (MQ=255)                  tgtgCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTc                                                                                                                        <  8:352759/136‑1 (MQ=255)                   gtgCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTc                                                                                                                        >  6:72216/1‑135 (MQ=255)                    tgCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGa                                                                                                                      <  5:227519/136‑1 (MQ=255)                     gtgcgcCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGCTCAAGGCTGGCGACTTGATGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGAt                                                                                                                     >  8:45781/3‑136 (MQ=255)                     gcgcgcCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCTAGGCCGCGCTCGAt                                                                                                                     >  5:119074/1‑136 (MQ=255)                     gcgcgcCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGAt                                                                                                                     >  1:232210/1‑136 (MQ=255)                     gcgcacCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGAt                                                                                                                     >  7:344196/1‑136 (MQ=255)                           cGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGc                                                                                                                >  2:4800/1‑135 (MQ=255)                              tgtgGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAgc                                                                                                            <  3:255889/136‑1 (MQ=255)                              tgtgGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAgc                                                                                                            <  5:206357/136‑1 (MQ=255)                                                  aTGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACAcgcg                                                                                         <  3:245789/135‑1 (MQ=255)                                                  aTGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACAcgc                                                                                          <  5:310240/134‑1 (MQ=255)                                                       tCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAg                                                                                    <  7:12892/135‑1 (MQ=255)                                                          ggCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCaa                                                                                 >  8:367187/1‑135 (MQ=255)                                                          ggCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCaa                                                                                 >  7:150970/1‑135 (MQ=255)                                                          ggCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAAc                                                                                >  4:236377/1‑136 (MQ=255)                                                          ggCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGCCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCaa                                                                                 >  1:363393/1‑135 (MQ=255)                                                                    gcgAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGCCCGAGACACGCGCCCAGCAACTGCGCTTgcg                                                                      >  2:183581/1‑136 (MQ=255)                                                                    gcgAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTgcg                                                                      >  5:154464/1‑136 (MQ=255)                                                                                  cATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGTGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAagca                                                         <  8:308350/135‑1 (MQ=255)                                                                                   aTGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAagcagc                                                       <  5:72216/136‑1 (MQ=255)                                                                                       cAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCgg                                                   >  2:358846/1‑136 (MQ=255)                                                                                       cAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCgg                                                   >  5:258929/1‑136 (MQ=255)                                                                                       cAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGAGAAGCAGCCCgg                                                   >  2:76559/1‑136 (MQ=255)                                                                                        aaCCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGc                                                  <  8:252830/136‑1 (MQ=255)                                                                                             gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGTGCTTCTAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                             >  1:251450/1‑136 (MQ=255)                                                                                             gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGacac                                                                                            >  4:20860/1‑89 (MQ=255)                                                                                             gcagGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGCTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACt                                             >  2:134855/1‑136 (MQ=255)                                                                                             gcagGTCAAGGCTGGCGAATTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGacac                                                                                            <  3:20860/89‑1 (MQ=255)                                                                                               atgTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCTCTTCCAGGCCGCGCTCGATCAGGCGAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGc                                           <  5:257969/134‑1 (MQ=255)                                                                                               agGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGTGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGc                                           <  3:82367/136‑1 (MQ=255)                                                                                               agGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGc                                           <  7:226603/136‑1 (MQ=255)                                                                                               agGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGc                                           <  4:95133/136‑1 (MQ=255)                                                                                                              aCTTGTTGATGAGTATTGACCTCGAGCGCTTCCAGGCAGAGCTCGATCAGGCCAGCGCGGTGACCGAGGCACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                            <  3:54033/136‑1 (MQ=255)                                                                                                              aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                            >  8:284606/1‑136 (MQ=255)                                                                                                              aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                            <  8:71066/136‑1 (MQ=255)                                                                                                              aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                            >  5:54431/1‑136 (MQ=255)                                                                                                              aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                            <  5:16130/136‑1 (MQ=255)                                                                                                              aCTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGc                            >  1:61784/1‑136 (MQ=255)                                                                                                                  gttgATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGa                          <  5:45883/134‑1 (MQ=255)                                                                                                                      atgaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGt                     <  4:252567/135‑1 (MQ=255)                                                                                                                        gaGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGcc                  >  1:218531/1‑136 (MQ=255)                                                                                                                           tATTGACCGCGAGCGCTTCCAGGCCGCGCTCGGTCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAg               >  3:104618/1‑136 (MQ=255)                                                                                                                             ttGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGtt             >  1:340660/1‑136 (MQ=255)                                                                                                                                aCCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCg          >  1:120233/1‑136 (MQ=255)                                                                                                                                  cgcgAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTg        <  8:150970/136‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGCTCAGGCCAGCGCGGTGACCGAGCCACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   <  4:48196/136‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGa                          >  6:176642/1‑113 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGa                          <  5:176642/113‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   <  7:61517/136‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   >  7:322165/1‑136 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   >  7:314554/1‑136 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   <  2:165747/136‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   <  5:12549/136‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   <  5:307206/136‑1 (MQ=255)                                                                                                                                       gcgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAAc   >  6:45770/1‑136 (MQ=255)                                                                                                                                        cgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACg  >  4:81253/1‑136 (MQ=255)                                                                                                                                        cgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACg  >  4:300395/1‑136 (MQ=255)                                                                                                                                        cgcTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCATGGGGCCAGAGGCGATCAGTGCAGAGTTGCGTGAAAACg  >  2:131729/1‑136 (MQ=255)                                                                                                                                        |                                                                                                                                        TGGACGCGTGATGCCCGTGTGCGCGCCGATGTGGTGGTGATTGCGCCTGATGTGTCCGGCTGGGTTCGCGAACTCAAGGTCCATGACAACCAGCAGGTCAAGGCTGGCGACTTGTTGATGAGTATTGACCGCGAGCGCTTCCAGGCCGCGCTCGATCAGGCCAGCGCGGTGACCGAGACACGCGCCCAGCAACTGCGCTTGCGCGAGCGTGAAGCAGCCCGGCGTACTGCCCTGGGGCCAGAGGCGATCAGTGCCGAGTTGCGTGAAAACG  >  ACB143/1311926‑1312196

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