Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,681,307	G→A	24.6%	A159A (GCG→GCA)	gloB →	hydroxyacylglutathione hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,681,307	0	G	A	24.6%	103.5 / 33.9	61	A159A (GCG→GCA)	gloB	hydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base A (6/9);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CCGGCCACACCCTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTGCGCGCCGAAAATCGCAT  >  NC_002947/4681159‑4681450                                                                                                                                                     |                                                                                                                                                ccGGCCACACCCTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTggcg                                                                                                                                                 >  8:109492/1‑149 (MQ=255)   ggCCACACCCTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTggcggc                                                                                                                                               <  4:321867/149‑1 (MQ=255)        cacCCTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGc                                                                                                                                          <  4:4925/149‑1 (MQ=255)         acCCTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGc                                                                                                                                          >  3:16578/1‑148 (MQ=255)           ccTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGc                                                                                                                                          >  6:111716/1‑146 (MQ=255)               gggCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCagag                                                                                                                                     <  6:318285/147‑1 (MQ=255)                ggCATATCTGCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAg                                                                                                                                  <  8:298611/149‑1 (MQ=255)                 gCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGa                                                                                                                                 <  7:109492/149‑1 (MQ=255)                     atCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGacac                                                                                                                              >  7:162657/1‑148 (MQ=255)                      tCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACaa                                                                                                                            <  8:211904/149‑1 (MQ=255)                        gCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAgt                                                                                                                          <  4:16578/149‑1 (MQ=255)                             cttcAGCGACCAGCCAGCAACGCCCATCCTGTGCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTa                                                                                                                       <  4:159136/147‑1 (MQ=255)                                cAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCAGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTgct                                                                                                                  >  3:254229/1‑148 (MQ=255)                                 aGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTgcgc                                                                                                                 >  6:14688/1‑149 (MQ=255)                                   cGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGcc                                                                                                                >  6:289142/1‑148 (MQ=255)                                    gACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCa                                                                                                              <  5:14688/149‑1 (MQ=255)                                    gACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCaga                                                                                                                                      <  6:206283/125‑1 (MQ=255)                                    gACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCaga                                                                                                                                      >  5:206283/1‑125 (MQ=255)                                         gccagcAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGaa                                                                                                          >  2:51382/1‑148 (MQ=255)                                          ccagcAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCAt                                                                                                             <  1:212224/144‑1 (MQ=255)                                          ccagcAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCAt                                                                                                             >  2:212224/1‑144 (MQ=255)                                          ccagcAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATa                                                                                                        >  3:192565/1‑149 (MQ=255)                                           cagcAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATac                                                                                                       <  1:351752/149‑1 (MQ=255)                                                  gCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCa                                                                                                <  5:289142/149‑1 (MQ=255)                                                    ccATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAg                                                                                               >  6:263735/1‑148 (MQ=255)                                                    ccATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAg                                                                                               >  2:108663/1‑148 (MQ=255)                                                          tGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTg                                                                                        >  5:175368/1‑149 (MQ=255)                                                          tGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCAATGAATACACCCTCAGCAACCTg                                                                                        >  4:226056/1‑149 (MQ=255)                                                              cAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTgcgc                                                                                     <  3:232227/148‑1 (MQ=255)                                                               aGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCAt                                                                                                             >  3:296102/1‑123 (MQ=255)                                                               aGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCAt                                                                                                             <  4:296102/123‑1 (MQ=255)                                                                gTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTc                                                                                  >  7:279186/1‑149 (MQ=255)                                                                gTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCATCAGCACCCTGCGCTTc                                                                                  >  8:248051/1‑149 (MQ=255)                                                                  ggCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCg                                                                                 <  7:59523/148‑1 (MQ=255)                                                                       cacCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAgg                                                                           <  1:108663/149‑1 (MQ=255)                                                                       cacCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAgg                                                                           <  4:19790/149‑1 (MQ=255)                                                                          ccTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTgcgc                                                                                     >  5:200079/1‑136 (MQ=255)                                                                          ccTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTgcgc                                                                                     <  6:200079/136‑1 (MQ=255)                                                                          ccTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCtg                                                                        >  5:20065/1‑149 (MQ=255)                                                                              ttCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGATCAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTgg                                                                     <  4:254229/148‑1 (MQ=255)                                                                                 gccgccGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAAc                                                                  <  6:20065/148‑1 (MQ=255)                                                                                  ccgccgGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCg                                                                <  5:263735/149‑1 (MQ=255)                                                                                   cgccgGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCATCAACCTGCGCTTCGCCAAGGCTGTGGAACCGa                                                               <  7:248051/149‑1 (MQ=255)                                                                                         ttGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAAcc                                                         <  5:19078/149‑1 (MQ=255)                                                                                            cGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGc                                                      <  8:279186/149‑1 (MQ=255)                                                                                              ggCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCAc                                                    <  6:175368/149‑1 (MQ=255)                                                                                              ggCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCAc                                                    <  1:51382/149‑1 (MQ=255)                                                                                               gcgcATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCa                                                     >  1:151275/1‑147 (MQ=255)                                                                                               gcgcATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCa                                                     <  2:151275/147‑1 (MQ=255)                                                                                                      ttCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGc                                             <  5:111716/148‑1 (MQ=255)                                                                                                                cacCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCAGCCCTGCAAGAGCAGACACAATTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGa                                   <  4:192565/148‑1 (MQ=255)                                                                                                                     gAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGt                             <  3:226056/149‑1 (MQ=255)                                                                                                                           aTGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAg                                              <  4:224306/126‑1 (MQ=255)                                                                                                                           aTGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAg                                              >  3:224306/1‑126 (MQ=255)                                                                                                                             gCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGc                                             <  2:77636/125‑1 (MQ=255)                                                                                                                             gCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGc                                             >  1:77636/1‑125 (MQ=255)                                                                                                                             gCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTAc                          <  2:154750/144‑1 (MQ=255)                                                                                                                             gCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTAc                          >  1:154750/1‑144 (MQ=255)                                                                                                                               aGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTg                   >  2:341804/1‑149 (MQ=255)                                                                                                                                 cccgcccTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTgcg                 >  4:146481/1‑149 (MQ=255)                                                                                                                                      ccTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTGCGCGcc             >  8:62791/1‑148 (MQ=255)                                                                                                                                      ccTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCAAGAATACACCCTCAGCAACCTGCGCTTCGCCAAGCCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTGCGCGcc             >  6:342373/1‑148 (MQ=255)                                                                                                                                             cgccTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTGCGCGCCGAAAATCg     >  1:68069/1‑149 (MQ=255)                                                                                                                                                  ggcagcCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTGCGCGCCGAAAATCGCAt  >  3:27467/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCGGCCACACCCTGGGGCATATCGCCTTCTTCAGCGACCAGCCAGCAACGCCCATCCTGTTCAGTGGCGACACCCTGTTCGCCGCCGGTTGCGGGCGCATGTTCGAAGGTACACCCGAACAGATGCAGCCCGCCCTCGCCCGCCTGGCGGCCCTGCCAGAGCAGACACAAGTGTACTGCGCCCATGAATACACCCTCAGCAACCTGCGCTTCGCCAAGGCTGTGGAACCGACCAACCCGCACGTTCAGCAACGCTTCGAGGACGTTACCCGTTTGCGCGCCGAAAATCGCAT  >  NC_002947/4681159‑4681450

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