Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,711,169	A→G	7.1%	A174A (GCT→GCC)	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,169	0	A	G	7.1%	159.9 / 2.8	56	A174A (GCT→GCC)	paaC	ring 1,2‑phenylacetyl‑CoA epoxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base G (2/2);  total (26/30)
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.76e-01

CATCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTGGCG  >  NC_002947/3711021‑3711313                                                                                                                                                              |                                                                                                                                                 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                                                                                                                                   >  8:322984/1‑121 (MQ=255)                                            gACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACg                                                                                                                                   <  7:322984/121‑1 (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                                                                                         >  8:118956/1‑149 (MQ=255)                                                          tCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTgg                                                                                         <  5:320866/149‑1 (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               >  4:324381/1‑147 (MQ=255)                                                                                                                                               ttCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGActtcctt               >  1:130292/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)                                                                                                                                                              |                                                                                                                                                 CATCCAGGTAGAAGTGGCTGGCGGCCTTGGGCAGCGGCAGTTCGACCGAAGCAAACGTCGCGGCGACCTGCTTCAGCCAGGTGGCGCCGACGGTCGCTGGGTCGGCGGCGATAC‑‑CG‑‑‑‑‑CAT‑‑TCCGCCAGACGGACTTCGTTGTCGCTGCCAGCAGCGAGTTCGACGGTGAAGCGCCACAGTGCCGGGACGGCGGCGAGCATGCGCTGGCGGCTGGCATCGGTACCACCGCCCAGACGTTGTACCCATTCGCTGGACCGGCGCAGGTGGTAGGTGACTTCCTTCAGGGCCTTGGCG  >  NC_002947/3711021‑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: ‑