Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,711,168	C→T	10.9%	G175S (GGC→AGC)	paaC ←	ring 1,2‑phenylacetyl‑CoA epoxidase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,711,168	0	C	T	10.9%	147.9 / 7.4	55	G175S (GGC→AGC)	paaC	ring 1,2‑phenylacetyl‑CoA epoxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (3/3);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CCATCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTGGCG  >  NC_002947/3711020‑3711313                                                                                                                                                              |                                                                                                                                                  ccATCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGcc                                                                                                                                                   <  5:8493/149‑1 (MQ=255)  cATCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCca                                                                                                                                                  <  6:74067/149‑1 (MQ=255)    tCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCcag                                                                                                                                                 >  5:149819/1‑148 (MQ=255)            gAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCTGGCAGCGaca                                                                                                                                         <  7:87413/148‑3 (MQ=255)               gTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGa                                                                                                                                     <  4:129880/149‑1 (MQ=255)                  gctggcGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACg                                                                                                                                   <  6:149819/148‑1 (MQ=255)                    tggcGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGt                                                                                                                                 >  3:289234/1‑148 (MQ=255)                        ggcCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAgc                                                                                                                            <  8:138307/149‑1 (MQ=255)                        ggcCGTGGGCAGCGGCAGTTCGACCGAAGCAAACGGCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAgc                                                                                                                            <  6:204790/149‑1 (MQ=255)                                   gCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGcc                                                                                                                 >  2:353978/1‑149 (MQ=255)                                    cGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCg                                                                                                                >  6:320866/1‑149 (MQ=255)                                         gTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGAcg                                                                                                           >  3:364517/1‑149 (MQ=255)                                             gACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACg                                                                                                                                   <  7:322984/121‑1 (MQ=255)                                             gACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACg                                                                                                                                   >  8:322984/1‑121 (MQ=255)                                               ccGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGgcggcg                                                                                                     <  4:364517/149‑1 (MQ=255)                                                cGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACggcggc                                                                                                      >  6:332959/1‑147 (MQ=255)                                                cGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGa                                                                                                    <  3:115574/149‑1 (MQ=255)                                                  aaGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGGTGCTAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAg                                                                                                   >  1:302968/1‑148 (MQ=255)                                                     cAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGTGAGCATg                                                                                               >  5:150492/1‑149 (MQ=255)                                                       aaCGTCGCGGCGACCTGCTTCAGCCAGGTGGCGACGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATgcg                                                                                             >  2:260075/1‑149 (MQ=255)                                                           tCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTgg                                                                                         <  5:320866/149‑1 (MQ=255)                                                           tCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTgg                                                                                         >  8:118956/1‑149 (MQ=255)                                                             gcggcgACCTGCTTCAGCCAGGTGGCGCGGACGGGCGCTGGGTGGGCGGGGATACCACGACCTTCATTATAACCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCAt                                                                                                <  1:119514/149‑1 (MQ=255)                                                             gcggcgACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTggc                                                                                        <  8:115683/148‑1 (MQ=255)                                                                     cTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCca                                                                                                                                                  <  7:180889/82‑1 (MQ=255)                                                                     cTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCca                                                                                                                                                  >  8:180889/1‑82 (MQ=255)                                                                       gCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCg                                                                             >  8:221913/1‑149 (MQ=255)                                                                           cAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGt                                                                           <  8:96320/147‑1 (MQ=255)                                                                            aGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTac                                                                         <  1:353978/148‑1 (MQ=255)                                                                              ccAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCgg                                                                            <  8:147679/143‑1 (MQ=255)                                                                              ccAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCgg                                                                            >  7:147679/1‑143 (MQ=255)                                                                                aGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCg                                                                    >  3:188193/1‑149 (MQ=255)                                                                                  gTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGcc                                                                  >  4:182844/1‑149 (MQ=255)                                                                                    ggCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGcc                                                                  >  3:313061/1‑147 (MQ=255)                                                                                      cgcCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGa                                                              <  2:308756/149‑1 (MQ=255)                                                                                       gcCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGAc                                                             <  5:332959/149‑1 (MQ=255)                                                                                        ccGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACg                                                            >  1:72978/1‑149 (MQ=255)                                                                                          gACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGtt                                                          >  4:295495/1‑149 (MQ=255)                                                                                             ggTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTa                                                       <  3:182844/149‑1 (MQ=255)                                                                                                                aTAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCTAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGgcg                                    <  2:302968/149‑1 (MQ=255)                                                                                                                 tAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGgcg                                    >  6:278257/1‑148 (MQ=255)                                                                                                                              aT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAggtgg                             <  7:221913/149‑1 (MQ=255)                                                                                                                               t‑‑tCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAggtggt                            <  5:107715/149‑1 (MQ=255)                                                                                                                                    cGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAg                          <  4:289234/148‑1 (MQ=255)                                                                                                                                     gCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGTGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGt                        <  6:150492/149‑1 (MQ=255)                                                                                                                                     gCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGt                        <  7:118956/149‑1 (MQ=255)                                                                                                                                         gACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGAct                    <  4:188193/149‑1 (MQ=255)                                                                                                                                          aCGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAgg                                <  8:86586/136‑1 (MQ=255)                                                                                                                                          aCGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAgg                                >  7:86586/1‑136 (MQ=255)                                                                                                                                            ggACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGActtcc                 >  1:368142/1‑149 (MQ=255)                                                                                                                                               cTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGActtccttc              <  2:72978/149‑1 (MQ=255)                                                                                                                                                ttCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGActtcctt               >  1:130292/1‑147 (MQ=255)                                                                                                                                                ttCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGActtcctt               >  4:324381/1‑147 (MQ=255)                                                                                                                                                   gaactCGCTGCTGGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAgg           <  5:68347/144‑1 (MQ=255)                                                                                                                                                   gaactCGCTGCTGGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAgg           >  6:68347/5‑148 (MQ=255)                                                                                                                                                   gTTGTCGCTGCCAGCAGCGAGTTCGACGGTGTAGCGCGACTGTGGCGGGACGGGGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTTGAGCGGCGCAGGTGGTAGGTGACTTGCTTCAgg           <  4:291123/148‑1 (MQ=255)                                                                                                                                                     actCGCTGCTGGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGcc        >  6:168525/3‑149 (MQ=255)                                                                                                                                                       tCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTCCCCATTCGCTGGACCGGCGCAGGTGGTTGGTGACTTCCTTCAGGGCCtt      >  4:272415/1‑149 (MQ=255)                                                                                                                                                          cTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTggc   <  1:260075/149‑1 (MQ=255)                                                                                                                                                           tGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTggcg  >  4:282987/1‑149 (MQ=255)                                                                                                                                                             ccAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTggcg  <  2:130292/147‑1 (MQ=255)                                                                                                                                                              |                                                                                                                                                  CCATCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTGGCG  >  NC_002947/3711020‑3711313

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