Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,494,816	T→C	6.9%	P466P (CCA→CCG)	puuD ←	uricase/urate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,494,816	0	T	C	6.9%	209.0 / 3.6	72	P466P (CCA→CCG)	puuD	uricase/urate oxidase
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base C (2/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGC  >  NC_002947/3494673‑3494964                                                                                                                                                |                                                                                                                                                     cTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCtt                                                                                                                                                 <  6:2882/149‑1 (MQ=255)            tgtgTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgtgt                                                                                                                                      <  6:394045/149‑5 (MQ=255)                  aTCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAaca                                                                                                                                >  3:230165/1‑149 (MQ=255)                  aTCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGCTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTACGTCACGGCCTTGTACACGCCTGCCTAGCCGGGCACTTCGAaca                                                                                                                                >  3:235308/1‑149 (MQ=255)                    cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAacac                                                                                                                               >  4:164560/1‑148 (MQ=255)                    cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAacac                                                                                                                               <  3:164560/148‑1 (MQ=255)                    cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgt                                                                                                                                        <  8:264541/139‑3 (MQ=255)                    cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgt                                                                                                                                        >  7:264541/1‑137 (MQ=255)                     ccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCt                                                                                                                             <  7:387165/149‑1 (MQ=255)                      cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCt                                                                                                                                                      <  6:51500/123‑1 (MQ=255)                      cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                            >  4:320322/1‑149 (MQ=255)                      cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                            >  4:325970/1‑149 (MQ=255)                      cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                            >  4:275091/1‑149 (MQ=255)                      cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                            >  2:191020/1‑149 (MQ=255)                       tGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCg                                                                                                                           >  7:300906/1‑149 (MQ=255)                              aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCGACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCggg                                                                                                                    >  5:359245/1‑149 (MQ=255)                              aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGgcgc                                                                                                                       <  4:361080/146‑1 (MQ=255)                              aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGgcgc                                                                                                                       >  3:361080/1‑146 (MQ=255)                              aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCggg                                                                                                                    <  5:5961/149‑1 (MQ=255)                                  gCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTc                                                                                                                <  3:24623/149‑1 (MQ=255)                                   ccTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTcg                                                                                                               >  6:276125/1‑149 (MQ=255)                                                 gcgcCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActt                                                                                                 >  4:393558/1‑149 (MQ=255)                                                  cgcCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtg                                                                                                <  6:67517/149‑1 (MQ=255)                                                   gcCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtgg                                                                                               >  4:294482/1‑149 (MQ=255)                                                     ccTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCgg                                                                                                                     >  3:398012/1‑125 (MQ=255)                                                     ccTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCgg                                                                                                                     <  4:398012/125‑1 (MQ=255)                                                         tgtcAGCTCGGCAACCAGCCGCAGCGAGGCTGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGTTCACGGCGTTGTTCACGCCTGGGTTGCCGTGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                          <  4:235308/145‑1 (MQ=255)                                                         tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttgg                                                                                         >  5:193584/1‑149 (MQ=255)                                                         tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                          <  4:230165/148‑1 (MQ=255)                                                         tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                          <  2:383153/148‑1 (MQ=255)                                                         tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                          >  1:383153/1‑148 (MQ=255)                                                         tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggcctt                                                                                           <  3:275091/147‑1 (MQ=255)                                                           gCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGc                                                                                                                   >  1:362408/1‑121 (MQ=255)                                                           gCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGc                                                                                                                   <  2:362408/121‑1 (MQ=255)                                                           gCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttggc                                                                                        <  3:393558/148‑1 (MQ=255)                                                             aGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGcctc                                                                                     >  3:250995/1‑149 (MQ=255)                                                             aGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGcctc                                                                                     >  3:31028/1‑149 (MQ=255)                                                               cTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgtgt                                                                                                                                      >  2:391273/1‑94 (MQ=255)                                                               cTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgtgt                                                                                                                                      <  1:391273/98‑5 (MQ=255)                                                                 cGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTAGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctgg                                                                                 >  5:78157/1‑149 (MQ=255)                                                                  ggCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctgg                                                                                 <  6:193584/148‑1 (MQ=255)                                                                   gCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctggct                                                                               >  3:139784/1‑149 (MQ=255)                                                                     aaCCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctggctgg                                                                             <  3:294482/149‑1 (MQ=255)                                                                     aaCCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctggctgg                                                                             >  3:33026/1‑149 (MQ=255)                                                                     aaCCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCatggctgg                                                                             >  5:9029/1‑149 (MQ=255)                                                                       ccAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCg                                                                           <  4:139784/149‑1 (MQ=255)                                                                        cAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCgg                                                                          >  2:170799/1‑149 (MQ=255)                                                                         aGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCgg                                                                          >  2:318661/1‑148 (MQ=255)                                                                           ccGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTAt                                                                       <  5:255653/149‑1 (MQ=255)                                                                               aGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTAtcgtc                                                                   <  3:325970/149‑1 (MQ=255)                                                                               aGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTAtcgtc                                                                   <  3:154076/149‑1 (MQ=255)                                                                               aGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTAtcgtc                                                                   <  3:320322/149‑1 (MQ=255)                                                                                    gtcGGTCAGCTGATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCGTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTTGCCTCCTGGCTGGCGGTATCGTCGagca                                                              <  7:374457/147‑1 (MQ=255)                                                                                        gTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACTGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCAc                                                          <  1:134195/149‑1 (MQ=255)                                                                                             cTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGc                                                     >  4:147606/1‑149 (MQ=255)                                                                                                aTCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTgg                                                  >  3:176833/1‑149 (MQ=255)                                                                                                  ccAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctg                                                                                  <  8:360285/115‑1 (MQ=255)                                                                                                         gTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttgg                                                                                         <  7:82685/101‑1 (MQ=255)                                                                                                         gTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttgg                                                                                         >  8:82685/1‑101 (MQ=255)                                                                                                         gTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCAc                                          <  3:147606/148‑1 (MQ=255)                                                                                                               cGGGCGCAAGTAGGTCACGGCCTGGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGa                                   <  8:300906/149‑1 (MQ=255)                                                                                                                      aaGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAgcg                            >  4:285337/1‑149 (MQ=255)                                                                                                                       aGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGgccagcca                                           >  7:279891/1‑133 (MQ=255)                                                                                                                       aGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGgccagcca                                           <  8:279891/133‑1 (MQ=255)                                                                                                                               aCGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATc                   >  4:304709/1‑149 (MQ=255)                                                                                                                                cGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATc                   >  3:368906/1‑148 (MQ=255)                                                                                                                                     ttGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATgcg             >  6:267642/1‑149 (MQ=255)                                                                                                                                        tACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATgcgc            >  7:198598/1‑147 (MQ=255)                                                                                                                                         acacGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAg         >  3:223014/1‑149 (MQ=255)                                                                                                                                         acacGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAg         >  3:228156/1‑149 (MQ=255)                                                                                                                                          cacGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAgg        <  3:228688/149‑1 (MQ=255)                                                                                                                                           acGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGc       >  2:92807/1‑149 (MQ=255)                                                                                                                                               cTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAg   <  3:304709/149‑1 (MQ=255)                                                                                                                                                tGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGc  >  2:177958/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CTGGTTACTCATGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGC  >  NC_002947/3494673‑3494964

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