Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,799,421	G→A	5.9%	A94A (GCG→GCA)	PP_3359 →	hydroxycinnamic acid degradation regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,799,421	0	G	A	5.9%	200.7 / 2.1	68	A94A (GCG→GCA)	PP_3359	hydroxycinnamic acid degradation regulator
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (2/2);  total (35/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.32e-01

CGGCCCATGACCTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAACCTCTCGGATGACG  >  NC_002947/3799276‑3799563                                                                                                                                                  |                                                                                                                                                cGGCCCATGACCTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCCCCTGCGTTGGCAATCGACCCCCCGCCAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTgg                                                                                                                                             >  1:300282/1‑150 (MQ=255)     ccATGACCTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCc‑gt                                                                                                                                        >  1:578854/1‑150 (MQ=255)     ccATGACCTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCc‑gt                                                                                                                                        >  1:1324585/1‑150 (MQ=255)        tGACCTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTgcg                                                                                                                                     <  1:479920‑M2/150‑4 (MQ=255)            cTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCc‑gtgcgcgtg                                                                                                                                 >  2:1518591/1‑150 (MQ=255)                        cAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCg                                                                                                                     <  1:513246/150‑1 (MQ=255)                         aaTTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCgg                                                                                                                    >  1:1181240/1‑150 (MQ=255)                             ttCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACaa                                                                                                                <  2:629968/150‑1 (MQ=255)                             ttCCGCCCTGCCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACaa                                                                                                                <  2:901558/150‑1 (MQ=255)                                cGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAgcg                                                                                                             <  1:1342852/150‑1 (MQ=255)                                    cTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCACGGCCA‑GCgcgcgctcttccagtttgttcagcattgg                                                                                                           <  1:899664‑M2/148‑30 (MQ=255)                                    cTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCAC‑GGCCAGCgcgcgctcttccagtttgttcagcattgg                                                                                                           >  2:899664‑M2/1‑119 (MQ=255)                                     tGGCGATCATAAAGATGAGCCCCGGGCCGAGACCTGCCCACCTGGCGCGTGCGTTGTCAATCGCCCCCCCGCAAGTAGTGCCAATGCTTAACAGACTGGAAGAGCGCGCGCGGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGt                                                                                                        <  2:300282/150‑1 (MQ=255)                                           tcatcGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAAGGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGa                                                                                                  >  1:108089/1‑150 (MQ=255)                                              tcGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCCATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGAttt                                                                                               >  1:31589/1‑150 (MQ=255)                                              tcGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGAttt                                                                                               >  2:1366626/1‑150 (MQ=255)                                              tcGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGAttt                                                                                               >  1:381576/1‑150 (MQ=255)                                              tcGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGAttt                                                                                               >  1:180439/1‑150 (MQ=255)                                                  ccAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCt                                                                                           <  2:381576/150‑1 (MQ=255)                                                  ccAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCt                                                                                           >  2:970330/1‑150 (MQ=255)                                                  ccAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCt                                                                                           <  1:1366626/150‑1 (MQ=255)                                                   cAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTg                                                                                                                                              <  2:534900/99‑1 (MQ=255)                                                   cAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTg                                                                                                                                              >  1:534900/1‑99 (MQ=255)                                                     gAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAg                                                                                        <  2:1181240/150‑1 (MQ=255)                                                     gAACACCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAg                                                                                        >  2:533979/1‑150 (MQ=255)                                                                  aTGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAg                                                                           <  2:578854/150‑1 (MQ=255)                                                                   tGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAgg                                                                          >  1:1516186/1‑150 (MQ=255)                                                                    gCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGc                                                                         >  2:223235/1‑150 (MQ=255)                                                                       ggCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGcc                                                                        <  1:1389159/148‑1 (MQ=255)                                                                       ggCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGcc                                                                        >  2:1389159/1‑148 (MQ=255)                                                                        gCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCt                                                                      <  1:751690/149‑1 (MQ=255)                                                                        gCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCt                                                                      >  2:751690/1‑149 (MQ=255)                                                                                  cgcgTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGt                                                           <  2:1069192/150‑1 (MQ=255)                                                                                      tgcgtTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCACGGCCA‑GCgcgcgctcttccagtttgtt                                                                                                                    >  2:885419‑M2/1‑69 (MQ=255)                                                                                      tgcgtTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCACGGCCA‑GCgcgcgctcttccagtttgtt                                                                                                                    <  1:885419‑M2/89‑21 (MQ=255)                                                                                       gcgtTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAACTCGGGCGAGGCCCTGTTGAAGGAGTTGAAg                                                      <  1:1311802/150‑1 (MQ=255)                                                                                       gcgtTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAg                                                      <  1:1143742/150‑1 (MQ=255)                                                                                       gcgtTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAg                                                      <  1:1143743/150‑1 (MQ=255)                                                                                          ttGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGcac                                                   <  2:689619/150‑1 (MQ=255)                                                                                           tGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGcaca                                                  >  1:1214575/1‑150 (MQ=255)                                                                                              cAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATcg                                               >  1:1497519/1‑150 (MQ=255)                                                                                               aaTCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAATCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGTAGTTGAAGGACATcgc                                              >  1:516169/1‑150 (MQ=255)                                                                                                      cccccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAgcgc                                                                                                            >  2:1521951/1‑81 (MQ=255)                                                                                                      cccccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAgcgc                                                                                                            <  1:1521951/81‑1 (MQ=255)                                                                                                      cccccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCt                                                                                           <  2:1582983/98‑1 (MQ=255)                                                                                                      cccccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCt                                                                                           >  1:1582983/1‑98 (MQ=255)                                                                                                        cccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGc                                                     >  2:1526545/1‑134 (MQ=255)                                                                                                        cccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGc                                                     <  1:1526545/134‑1 (MQ=255)                                                                                                         ccGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGc                                    >  2:1152890/1‑150 (MQ=255)                                                                                                                  gTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAg                           >  1:841873/1‑150 (MQ=255)                                                                                                                    gCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGc                         <  2:841873/150‑1 (MQ=255)                                                                                                                       aaTGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGa                                        >  1:280180/1‑132 (MQ=255)                                                                                                                       aaTGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGa                                        <  2:280180/132‑1 (MQ=255)                                                                                                                       aaTGCGGACCAAACTGGAAGAGCGCTCGCTGGCC‑GTTCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGATCAAATCGGGCTAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGAc                      <  2:702414/150‑1 (MQ=255)                                                                                                                          gCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATc                   <  2:1214575/150‑1 (MQ=255)                                                                                                                           cTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTAt                                                                                                      >  1:1072982/1‑66 (MQ=255)                                                                                                                           cTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTAt                                                                                                      <  2:1072982/66‑1 (MQ=255)                                                                                                                            tGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGTCGACAGCGATCACCAGGCGACATCCa                 >  2:1362272/1‑150 (MQ=255)                                                                                                                              aaCAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAAc               <  1:970330/150‑1 (MQ=255)                                                                                                                              aaCAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAAc               >  2:617887/1‑150 (MQ=255)                                                                                                                                 aaaCTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAACctc            >  1:840160/1‑150 (MQ=255)                                                                                                                                  aaCTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAACctct           >  1:1398778/1‑150 (MQ=255)                                                                                                                                   aCTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACa                                      >  1:871868/1‑122 (MQ=255)                                                                                                                                   aCTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACa                                      <  2:871868/122‑1 (MQ=255)                                                                                                                                   aCTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCGGAGCAAATCGAGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCTATCACCAGGCGACATCCAACCtctc          <  2:180439/150‑1 (MQ=255)                                                                                                                                       gAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAACCTCTCGGAt      <  2:108089/150‑1 (MQ=255)                                                                                                                                       gAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAACCTCTCGGAt      >  2:125922/1‑150 (MQ=255)                                                                                                                                           agCGCGCGCTGGCC‑GTgcgcgtgcgctcttccagtttgttcagcattggcacgacttgcggggggtcgattgccaacgcacgcgccaggtcggcctgcatcagcccggggttctggtcgatgatcgccagggcggaaaattgggccgggc  <  1:1234058‑M2/150‑135 (MQ=255)                                                                                                                                                  |                                                                                                                                                CGGCCCATGACCTGCGCCCGGCCCAATTTTCCGCCCTGGCGATCATCGACCAGAACCCCGGGCTGATGCAGGCCGACCTGGCGCGTGCGTTGGCAATCGACCCCCCGCAAGTCGTGCCAATGCTGAACAAACTGGAAGAGCGCGCGCTGGCC‑GTGCGCGTGCGGTGCAAACCGGACAAGCGCTCGTATGGGATTTTCCTGAGCAAATCGGGCGAGGCCCTGTTGAAGGAGTTGAAGCACATCGCCGCCGACAGCGATCACCAGGCGACATCCAACCTCTCGGATGACG  >  NC_002947/3799276‑3799563

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