Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,699:1	+A	8.2%	coding (201/831 nt)	PP_3358 ←	hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,798,699	1	.	A	8.2%	233.7 / 14.9	61	coding (201/831 nt)	PP_3358	hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base . (26/30);  new base A (2/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

ACCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGC  >  NC_002947/3798551‑3798847                                                                                                                                                      |                                                                                                                                                      aCCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCcggcg                                                                                                                                                        <  4:293607/149‑1 (MQ=255) aCCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCcggcg                                                                                                                                                        >  1:180467/1‑149 (MQ=255)   cAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGTGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCcggcg‑g                                                                                                                                                      >  8:302051/1‑148 (MQ=255)       cGTTGACCATGGCGATGGTCGGTTAGTCATACACACGCAGCAACTTCCCTTCCCATTGCGAGGCTGCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCCTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCTGCG‑GTCCAg                                                                                                                                                 <  8:130157/149‑1 (MQ=255)       aaccgACCATCGCCATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAg                                                                                                                                                 >  5:223157/5‑149 (MQ=255)             ccATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGACCATGCCGGCG‑GTCCAGGACTc                                                                                                                                            <  5:235303/148‑1 (MQ=255)                    gATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGC‑CGGCat                                                                                                                                   <  6:84667/149‑2 (MQ=255)                         tCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCa                                                                                                                                                  <  5:159140/130‑1 (MQ=255)                         tCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCa                                                                                                                                                  >  6:159140/1‑130 (MQ=255)                           ggTTTGGCTTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCAc                                                                                                                                   <  2:233356/142‑1 (MQ=255)                           ggTTTGGCTTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCAc                                                                                                                                   <  2:238442/142‑1 (MQ=255)                           ggTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCAc                                                                                                                                   >  1:233356/1‑142 (MQ=255)                           ggTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCAc                                                                                                                                   >  1:238442/1‑142 (MQ=255)                                 gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCa                                                                                                                       <  5:310985/148‑1 (MQ=255)                                 gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCa                                                                                                                       <  5:316659/148‑1 (MQ=255)                                  cATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGc                                                                                                                     <  3:137326/149‑1 (MQ=255)                                   aTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGc‑c                                                                                                                                        >  8:60421/1‑129 (MQ=255)                                   aTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGc‑c                                                                                                                                        <  7:60421/129‑1 (MQ=255)                                     aCAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCa                                                                                                                             >  6:178554/1‑138 (MQ=255)                                     aCAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCa                                                                                                                             <  5:178554/138‑1 (MQ=255)                                       aGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGc                                                                                                                 >  8:101353/1‑148 (MQ=255)                                        gACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGc                                                                                                                 >  5:231324/1‑147 (MQ=255)                                                   ttCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGAGTTCCGGGCCGGCGTCCACATCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCagcgtcagcgtc                                                                                                    >  2:234897/1‑149 (MQ=255)                                                      cattgccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTg                                                                                                 >  2:170020/1‑149 (MQ=255)                                                           ccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTc                                                                                              >  6:333178/1‑147 (MQ=255)                                                           ccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCaa                                                                                            <  1:170020/149‑1 (MQ=255)                                                                 cGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GGCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGtt                                                                                       <  8:176577/148‑1 (MQ=255)                                                                  gAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                     <  6:231324/149‑1 (MQ=255)                                                                  gAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                     >  1:153273/1‑149 (MQ=255)                                                                    ggCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCa                                                                                   <  2:180467/149‑1 (MQ=255)                                                                      cTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGc                                                                                 <  2:153273/149‑1 (MQ=255)                                                                       ttCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCa                                                                                >  2:350932/1‑149 (MQ=255)                                                                       ttCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGATTCCAGCa                                                                                >  3:57623/1‑149 (MQ=255)                                                                         cgcgACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCa                                                                                <  1:233957/147‑1 (MQ=255)                                                                          gcgACGAATCTTTTCCTGGAGGATTTCCGGGCCGGTGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGt                                                                             <  7:302051/149‑1 (MQ=255)                                                                           cgacgaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTc                                                                            <  7:101353/149‑1 (MQ=255)                                                                              cgaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         <  4:57623/149‑1 (MQ=255)                                                                                          tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               >  1:307059/1‑147 (MQ=255)                                                                                          tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               <  2:307059/147‑1 (MQ=255)                                                                                          tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGtt                                                             >  3:277316/1‑149 (MQ=255)                                                                                           ggaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGC‑CGGc                                                                                                                                     <  4:175031/76‑1 (MQ=37)                                                                                           ggaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGC‑CGGc                                                                                                                                     >  3:175031/1‑76 (MQ=37)                                                                                                aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGt                                                       <  4:277316/149‑1 (MQ=255)                                                                                                  ttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTcc                                                                                    <  8:325681/118‑1 (MQ=255)                                                                                                  ttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTcc                                                                                    >  7:325681/1‑118 (MQ=255)                                                                                                   tCCGGGCCGGCGTCCTCCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTgg                                                     <  6:223157/148‑1 (MQ=255)                                                                                                         ccGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCAt                                              <  5:333178/149‑1 (MQ=255)                                                                                                            gCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGc                                           <  1:350932/149‑1 (MQ=255)                                                                                                               tCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCAtt                                        >  5:338533/1‑149 (MQ=255)                                                                                                                  aCCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                      >  8:305689/1‑132 (MQ=255)                                                                                                                  aCCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                      <  7:305689/132‑1 (MQ=255)                                                                                                                         ggAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCACGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTc                               >  1:250190/1‑149 (MQ=255)                                                                                                                               aCTCCTTCAGGTCCATGCCGGCG‑GTCCACGCCTCGC‑CGGCACCGGTCAATACCAGCACGCCCGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTACCGGTACAGGGTGGGGCTCATGGCATTGCCTTTTTCCGGGCGa                        >  4:18568/1‑149 (MQ=255)                                                                                                                                cTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGAt                       <  5:136057/149‑1 (MQ=255)                                                                                                                                   cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGAtt                      <  6:338533/147‑1 (MQ=255)                                                                                                                                         gTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                            >  7:171877/1‑135 (MQ=255)                                                                                                                                         gTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTc               <  3:322748/148‑1 (MQ=255)                                                                                                                                         gTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTc               >  4:322748/1‑148 (MQ=255)                                                                                                                                         gTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCa              >  7:72563/1‑149 (MQ=255)                                                                                                                                         gTCCATGCCGGCG‑GTCCACGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                            <  8:171877/135‑1 (MQ=255)                                                                                                                                           ccATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAcc            >  7:256246/1‑149 (MQ=255)                                                                                                                                             aTGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAcc            >  2:294037/1‑147 (MQ=255)                                                                                                                                                ccGGCGAGTCCTGGAC‑CGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                            <  1:121924/128‑1 (MQ=255)                                                                                                                                                ccGGCGAGTCCTGGAC‑CGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                            >  2:121924/1‑128 (MQ=255)                                                                                                                                                 cggcg‑gTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGc       >  4:242462/1‑148 (MQ=255)                                                                                                                                                     g‑gTCCAGGACTCGC‑CGGCACCGGTCAATCCCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGc  <  3:18568/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      ACCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGC  >  NC_002947/3798551‑3798847

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