Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,494,825	G→A	6.2%	P463P (CCC→CCT)	puuD ←	uricase/urate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,494,825	0	G	A	6.2%	244.2 / 2.6	80	P463P (CCC→CCT)	puuD	uricase/urate oxidase
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base A (2/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGCTGCCGATC  >  NC_002947/3494684‑3494972                                                                                                                                              |                                                                                                                                                    tgtgTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgtgt                                                                                                                                              <  6:394045/149‑5 (MQ=255)       aTCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAaca                                                                                                                                        >  3:230165/1‑149 (MQ=255)       aTCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGCTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTACGTCACGGCCTTGTACACGCCTGCCTAGCCGGGCACTTCGAaca                                                                                                                                        >  3:235308/1‑149 (MQ=255)         cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAacac                                                                                                                                       <  3:164560/148‑1 (MQ=255)         cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAacac                                                                                                                                       >  4:164560/1‑148 (MQ=255)         cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgt                                                                                                                                                <  8:264541/139‑3 (MQ=255)         cccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgt                                                                                                                                                >  7:264541/1‑137 (MQ=255)          ccTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCt                                                                                                                                     <  7:387165/149‑1 (MQ=255)           cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                                    >  2:191020/1‑149 (MQ=255)           cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                                    >  4:275091/1‑149 (MQ=255)           cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                                    >  4:320322/1‑149 (MQ=255)           cTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTc                                                                                                                                    >  4:325970/1‑149 (MQ=255)            tGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCg                                                                                                                                   >  7:300906/1‑149 (MQ=255)                   aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCGACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCggg                                                                                                                            >  5:359245/1‑149 (MQ=255)                   aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGgcgc                                                                                                                               >  3:361080/1‑146 (MQ=255)                   aaCGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGgcgc                                                                                                                               <  4:361080/146‑1 (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                                                                                                                             <  4:398012/125‑1 (MQ=255)                                          ccTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCgg                                                                                                                             >  3:398012/1‑125 (MQ=255)                                              tgtcAGCTCGGCAACCAGCCGCAGCGAGGCTGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGTTCACGGCGTTGTTCACGCCTGGGTTGCCGTGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                                  <  4:235308/145‑1 (MQ=255)                                              tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttgg                                                                                                 >  5:193584/1‑149 (MQ=255)                                              tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                                  >  1:383153/1‑148 (MQ=255)                                              tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                                  <  2:383153/148‑1 (MQ=255)                                              tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttg                                                                                                  <  4:230165/148‑1 (MQ=255)                                              tgGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggcctt                                                                                                   <  3:275091/147‑1 (MQ=255)                                                gCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGc                                                                                                                           <  2:362408/121‑1 (MQ=255)                                                gCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGc                                                                                                                           >  1:362408/1‑121 (MQ=255)                                                gCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttggc                                                                                                <  3:393558/148‑1 (MQ=255)                                                  aGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGcctc                                                                                             >  3:31028/1‑149 (MQ=255)                                                  aGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGcctc                                                                                             >  3:250995/1‑149 (MQ=255)                                                    cTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgtgt                                                                                                                                              <  1:391273/98‑5 (MQ=255)                                                    cTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCCGGCAAGCCAGGcgtgt                                                                                                                                              >  2:391273/1‑94 (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:33026/1‑149 (MQ=255)                                                          aaCCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCctggctgg                                                                                     <  3:294482/149‑1 (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                                                                                                 >  8:82685/1‑101 (MQ=255)                                                                                              gTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAActtggccttgg                                                                                                 <  7:82685/101‑1 (MQ=255)                                                                                              gTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCAc                                                  <  3:147606/148‑1 (MQ=255)                                                                                                    cGGGCGCAAGTAGGTCACGGCCTGGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGa                                           <  8:300906/149‑1 (MQ=255)                                                                                                           aaGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAgcg                                    >  4:285337/1‑149 (MQ=255)                                                                                                            aGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGgccagcca                                                   <  8:279891/133‑1 (MQ=255)                                                                                                            aGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGgccagcca                                                   >  7:279891/1‑133 (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)                                                                                                                                       gCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGCTg        >  5:27247/1‑149 (MQ=255)                                                                                                                                         ttGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCg                                <  2:131731/123‑1 (MQ=255)                                                                                                                                         ttGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCg                                >  1:131731/1‑123 (MQ=255)                                                                                                                                            ccGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGCTGCCGAt   <  4:250995/149‑1 (MQ=255)                                                                                                                                             cGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCTCGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATgc                      <  1:289729/129‑1 (MQ=255)                                                                                                                                             cGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATgc                      >  2:289729/1‑129 (MQ=255)                                                                                                                                             cGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGCTGCCGATc  <  8:198598/149‑1 (MQ=255)                                                                                                                                             cGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGCTGCCGATc  <  4:31028/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    TGTGTCATCCCTGTACACAACGGCCTTTCAGCCGCGCGCGCCCTGTGGCAGCTCGGCAACCAGCCGCAGCGAGGCGGTCAGCTCATCCAGCTGGTAGTGCGGGCGCAAGTAGGTCACGGCCTTGTACACGCCTGGCTTGCCGGGCACTTCGAACACCTCGGCGCGGGCTTCGCGCAGCGGGAACTTGGCCTTGGCCTCCTGGCTGGCGGTATCGTCGAGCAGCACGTAGCTGGCCAGCCACTGGTTGAGGAAGCGCTCGACATCCATGCGCGAGGCGAAGCTGCCGATC  >  NC_002947/3494684‑3494972

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