Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,709	Δ1 bp	8.1%	coding (191/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,709	0	T	.	8.1%	189.8 / 10.7	62	coding (191/831 nt)	PP_3358	hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base . (2/3);  total (29/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 = 8.34e-01

CCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAG  >  NC_002947/3798563‑3798857                                                                                                                                                    |                                                                                                                                                      ccATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGACCATGCCGGCG‑GTCCAGGACTc                                                                                                                                                      <  5:235303/148‑1 (MQ=255)        gATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGC‑CGGCat                                                                                                                                             <  6:84667/149‑2 (MQ=255)               ggTTTGGCTTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCAc                                                                                                                                             <  2:238442/142‑1 (MQ=255)               ggTTTGGCTTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCAc                                                                                                                                             <  2:233356/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                                                                                                                                       <  5:178554/138‑1 (MQ=255)                         aCAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCa                                                                                                                                       >  6:178554/1‑138 (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                                                                         <  2:307059/147‑1 (MQ=255)                                                                              tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                                         >  1:307059/1‑147 (MQ=255)                                                                              tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGtt                                                                       >  3:277316/1‑149 (MQ=255)                                                                               ggaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGC‑CGGc                                                                                                                                               >  3:175031/1‑76 (MQ=37)                                                                               ggaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGC‑CGGc                                                                                                                                               <  4:175031/76‑1 (MQ=37)                                                                                    aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGt                                                                 <  4:277316/149‑1 (MQ=255)                                                                                      ttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTcc                                                                                              >  7:325681/1‑118 (MQ=255)                                                                                      ttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTcc                                                                                              <  8:325681/118‑1 (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                         >  4:322748/1‑148 (MQ=255)                                                                                                                             gTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTc                         <  3:322748/148‑1 (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                                      >  2:121924/1‑128 (MQ=255)                                                                                                                                    ccGGCGAGTCCTGGAC‑CGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                                      <  1:121924/128‑1 (MQ=255)                                                                                                                                     cggcg‑gTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGc                 >  4:242462/1‑148 (MQ=255)                                                                                                                                         g‑gTCCAGGACTCGC‑CGGCACCGGTCAATCCCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGc            <  3:18568/149‑1 (MQ=255)                                                                                                                                             ccAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCg         >  7:181166/1‑149 (MQ=255)                                                                                                                                                  aCTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTcc    <  8:181166/149‑1 (MQ=255)                                                                                                                                                   cTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTc     >  4:163411/1‑147 (MQ=255)                                                                                                                                                    tCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAg  >  2:201545/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAG  >  NC_002947/3798563‑3798857

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