Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,704	A→T	8.0%	W66R (TGG→AGG)	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,704	0	A	T	8.0%	177.9 / 4.7	62	W66R (TGG→AGG)	PP_3358	hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (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.71e-01

CGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCG  >  NC_002947/3798557‑3798850                                                                                                                                                     |                                                                                                                                                    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                                                                                                                                                     >  6:159140/1‑130 (MQ=255)                   tCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCa                                                                                                                                                     <  5:159140/130‑1 (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:316659/148‑1 (MQ=255)                           gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCa                                                                                                                          <  5:310985/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                                                                                        >  1:153273/1‑149 (MQ=255)                                                            gAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                        <  6:231324/149‑1 (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                                                                                                                                        <  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                                                         <  7:305689/132‑1 (MQ=255)                                                                                                            aCCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                         >  8:305689/1‑132 (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                               >  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)                                                                                                                                                     |                                                                                                                                                    CGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGC‑CGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCG  >  NC_002947/3798557‑3798850

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