Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	53,025	C→A	6.7%	P279P (CCC→CCA)	CJ019_00057 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	53,025	0	C	A	6.7%	242.3 / 5.4	74	P279P (CCC→CCA)	CJ019_00057	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/39);  new base A (2/3);  total (32/42)
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.03e-01

CGGCCCGGGAAAGCGCGCTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGACCGCCGCTGCGGCGCGCAA  >  ACB122/52896‑53153                                                                                                                                  |                                                                                                                                 cGGCCCGGGAAAGCGCGCTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGc                                                                                                                             <  6:368569/135‑1 (MQ=255) cGGCCCGGGAAAGCGCGCTACATTCTTGCGTACCCTGCCAGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGc                                                                                                                             <  2:113705/135‑1 (MQ=255)  ggCCCGGGAAAGCGCGCTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCg                                                                                                                            <  1:175243/135‑1 (MQ=255)      cGGGAAAGCGCGCTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGGGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATa                                                                                                                        <  2:304580/135‑1 (MQ=255)       gggAAAGCGCGCTACGTTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAg                                                                                                                       >  2:148703/1‑135 (MQ=255)                 gcTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTgg                                                                                                             >  7:98206/1‑135 (MQ=255)                 gcTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTgg                                                                                                             >  5:94716/1‑135 (MQ=255)                       ttCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGAc                                                                                                       <  8:262724/135‑1 (MQ=255)                       ttCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGAc                                                                                                       <  4:113399/135‑1 (MQ=255)                          ttGCGTACCCTGCCCGATTGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGt                                                                                                     <  4:219633/134‑1 (MQ=255)                          ttGCGTACCCTGCCCGATGGCTGAAGCGAAACGCGCTTCGAGGCCTTCGACATGCTGCACGTGGTTCCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGt                                                                                                     <  2:325683/134‑1 (MQ=255)                          ttGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGt                                                                                                     <  8:8744/134‑1 (MQ=255)                          ttGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGt                                                                                                     <  7:368744/134‑1 (MQ=255)                          ttGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGt                                                                                                     <  7:258837/134‑1 (MQ=255)                          ttGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGt                                                                                                     <  6:44811/134‑1 (MQ=255)                              gTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGTGGCCTTCGACATGCTGCACGTGGGACCGCCAGAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATc                                                                                                <  4:32738/135‑1 (MQ=255)                              gTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATc                                                                                                <  1:148703/135‑1 (MQ=255)                              gTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATc                                                                                                <  6:94716/135‑1 (MQ=255)                              gTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCACCGGCTGGGTTGACGt                                                                                                     <  3:272018/130‑1 (MQ=255)                              gTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCACCGGCTGGGTTGACGt                                                                                                     >  4:272018/1‑130 (MQ=255)                              gTAACCTGCCCGATGGCAGCATCTAAAAGAGCATCTAGGCCTTCGACATGCTGCACGTGGTCCCGACACAGGTGGCGACAGACTTCATTCGCGAAATCACATTTGAGGATAGCTCATGCTTGGTTGACTTCGATc                                                                                                <  4:28150/135‑1 (MQ=255)                                aCCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATcc                                                                                               >  1:56428/1‑134 (MQ=255)                                  ccTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGc                                                                                             >  6:141878/1‑134 (MQ=255)                                  ccTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGc                                                                                             >  4:118323/1‑134 (MQ=255)                                     gCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACCGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATAc                                                                                         >  3:257568/1‑135 (MQ=255)                                       ccGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATAccc                                                                                       >  3:207312/1‑135 (MQ=255)                                       ccGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATAccc                                                                                       >  2:224292/1‑135 (MQ=255)                                       ccGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATAccc                                                                                       >  1:223263/1‑135 (MQ=255)                                        cGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCt                                                                                      <  1:237100/135‑1 (MQ=255)                                        cGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCAGACCCt                                                                                      >  5:2086/1‑135 (MQ=255)                                         gATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTg                                                                                     >  3:331323/1‑135 (MQ=255)                                          aTGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTgc                                                                                    <  4:178501/135‑1 (MQ=255)                                               agcagcGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCAt                                                                               <  7:115336/135‑1 (MQ=255)                                               agcagcGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCAt                                                                               <  7:199700/135‑1 (MQ=255)                                                  agcGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGc                                                                            <  3:160189/135‑1 (MQ=255)                                                   gcGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTGCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGc                                                                            >  3:277676/1‑134 (MQ=255)                                                        aCGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATac                                                                      <  3:118323/135‑1 (MQ=255)                                                         cgcgCATCGAGTCCTTCGAAATGCTGCACGTGGTCACTCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATaca                                                                     <  8:140094/135‑1 (MQ=255)                                                          gcgcATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATacac                                                                    <  2:56428/135‑1 (MQ=255)                                                          gcgcATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATacac                                                                    <  1:234375/135‑1 (MQ=255)                                                              aTCGAGGCCTTCGACATGCTGCACGGGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAAc                                                                <  4:282980/135‑1 (MQ=255)                                                               tCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACg                                                               <  4:231418/135‑1 (MQ=255)                                                                cGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGt                                                              <  8:323048/135‑1 (MQ=255)                                                                             aTGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGt                                                 <  7:281951/135‑1 (MQ=255)                                                                             aTGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGt                                                 <  7:216444/135‑1 (MQ=255)                                                                              tgctgcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTg                                                <  5:141878/135‑1 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAACTACACCCACGTTCATGCCCTAGGTGACGt                                            >  8:139372/1‑135 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGt                                            >  3:186305/1‑135 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGt                                            >  8:293896/1‑135 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCATTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACg                                                               >  6:338577/1‑116 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCATTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACg                                                               >  6:333003/1‑116 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCATTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACg                                                               <  5:338577/116‑1 (MQ=255)                                                                                  gcACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCATTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACg                                                               <  5:333003/116‑1 (MQ=255)                                                                                   cACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGc                                                                            >  5:364389/1‑102 (MQ=255)                                                                                   cACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGc                                                                            <  6:364389/102‑1 (MQ=255)                                                                                      gTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGc                                         >  4:166030/1‑134 (MQ=255)                                                                                        ggTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTaa                                      <  5:162518/135‑1 (MQ=255)                                                                                            ccGCCACAGGTGGCGCCAGCCTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAaca                                    >  4:286903/1‑133 (MQ=255)                                                                                              gCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCa                                 >  4:88084/1‑134 (MQ=255)                                                                                                cacaGGTGGCGCCAGAATTCATTAGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCa                              >  7:15321/1‑135 (MQ=255)                                                                                                     gTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTg                                           <  8:278382/117‑1 (MQ=255)                                                                                                     gTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTg                                           >  7:278382/1‑117 (MQ=255)                                                                                                     gTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGc                          >  3:292362/1‑134 (MQ=255)                                                                                                       ggCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCaa                                                                          <  7:353879/84‑1 (MQ=255)                                                                                                       ggCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCaa                                                                          >  8:353879/1‑84 (MQ=255)                                                                                                       ggCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCaa                       <  4:207312/135‑1 (MQ=255)                                                                                                           ccAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGAcc                   >  3:352764/1‑135 (MQ=255)                                                                                                           ccAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGAcc                   >  7:48960/1‑135 (MQ=255)                                                                                                           ccAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGAcc                   >  5:315033/1‑135 (MQ=255)                                                                                                             aGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGAccgc                 >  7:157435/1‑135 (MQ=255)                                                                                                                    aTTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGACCGCCGCTg            <  4:186305/133‑1 (MQ=255)                                                                                                                      tCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGACAAGACCGCCGCTgcggc        >  6:160018/1‑135 (MQ=255)                                                                                                                          gAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGACCGCCGCTGCGgcgcgc    <  4:277676/135‑1 (MQ=255)                                                                                                                            aaGCCCCTTGGCGGATAGCTCCGGCTGGGGTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGGGACGTGGCTAACACCAGCAATGCCAAGACCGCCGCTGCGGCGCGCaa  <  4:300597/135‑1 (MQ=255)                                                                                                                                  |                                                                                                                                 CGGCCCGGGAAAGCGCGCTACATTCTTGCGTACCCTGCCCGATGGCAGCAGCGAAACGCGCATCGAGGCCTTCGACATGCTGCACGTGGTACCGCCACAGGTGGCGCCAGACTTCATTCGCGAAAGCCCCTTGGCGGATAGCTCCGGCTGGGTTGACGTCGATCCGCATACCCTGCGCCATCGCAAATACACCAACGTTCATGCCCTAGGTGACGTGGCTAACACCAGCAATGCCAAGACCGCCGCTGCGGCGCGCAA  >  ACB122/52896‑53153

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