Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,874,820	A→C	7.5%	S451A (TCG→GCG) ‡	ilvA‑II ←	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,874,820	0	A	C	7.5%	192.9 / 5.2	64	S451A (TCG→GCG) ‡	ilvA‑II	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base C (3/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGT  >  NC_002947/5874673‑5874966                                                                                                                                                    |                                                                                                                                                    tGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                     <  2:26503/148‑1 (MQ=255)      gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGtt                                                                                                                                               >  6:295937/1‑149 (MQ=255)       ccGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                              >  7:374961/1‑149 (MQ=255)       ccGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                              >  2:124515/1‑149 (MQ=255)        cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                              >  2:27892/1‑148 (MQ=255)         ggggTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCa                                                                                                                                            <  6:37067/149‑1 (MQ=255)           gTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAgc                                                                                                                                          >  8:176398/1‑149 (MQ=255)              ttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAgc‑gc                                                                                                                                       <  1:153673/148‑1 (MQ=255)              ttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAgc‑gc                                                                                                                                       <  1:148994/148‑1 (MQ=255)              ttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAgc‑gc                                                                                                                                       <  4:6013/148‑1 (MQ=255)              ttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGc‑gcc                                                                                                                                      >  7:278793/1‑149 (MQ=255)                 ggTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGc‑gcccgc                                                                                                                                   >  1:99557/1‑149 (MQ=255)                 ggTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGc‑gcccg                                                                                                                                    <  1:236753/148‑1 (MQ=255)                       gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCa                                                                                                                                <  1:118663/146‑1 (MQ=255)                             gTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGa                                                                                                                       <  6:144859/149‑1 (MQ=255)                             gTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGGCACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTCCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGa                                                                                                                       <  7:388761/149‑1 (MQ=255)                              tACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAg                                                                                                                      <  5:295937/149‑1 (MQ=255)                                 gggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGaa                                                                                                                   <  2:75292/149‑1 (MQ=255)                                   gTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                 >  5:426509/1‑149 (MQ=255)                                    tAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGtt                                                                                                                <  7:176398/149‑1 (MQ=255)                                       ccGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                             >  7:29425/1‑149 (MQ=255)                                       ccGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                             <  5:331220/149‑1 (MQ=255)                                        cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                             >  2:206955/1‑148 (MQ=255)                                               ggCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                           <  2:392854/143‑1 (MQ=255)                                               ggCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAg                                                                                                          <  6:426509/144‑1 (MQ=255)                                                  cAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCttgtt                                                                                                                         >  1:383020/1‑126 (MQ=255)                                                  cAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCttgtt                                                                                                                         <  2:383020/126‑1 (MQ=255)                                                          ccGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTc                                                                                           <  8:29425/148‑1 (MQ=255)                                                                 cAGGGGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCCCGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCACTCGGGGAAtt                                                                                   <  6:274833/149‑1 (MQ=255)                                                                     tGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCg                                                                                 >  6:32170/1‑147 (MQ=255)                                                                         gTTCGTCTTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAgcgc                                                                           <  1:123043/149‑1 (MQ=255)                                                                         gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAgcg                                                                            >  7:365725/1‑148 (MQ=255)                                                                          ttCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgca                                                                          <  2:99557/149‑1 (MQ=255)                                                                            cGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgcag                                                                         <  3:326866/148‑1 (MQ=255)                                                                                ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTc                                                                                                                                              >  3:50728/1‑76 (MQ=255)                                                                                ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTc                                                                                                                                              <  4:50728/76‑1 (MQ=255)                                                                                 tccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAt                                                                   >  7:25841/1‑149 (MQ=255)                                                                                   cGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAttt                                                                 >  5:291719/1‑149 (MQ=255)                                                                                    ggCACCTGCAGCCCCGCCACCACGCTCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtc                                                                <  4:360313/149‑1 (MQ=255)                                                                                    ggCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtc                                                                <  1:55399/149‑1 (MQ=255)                                                                                      cACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgt                                                              <  1:27892/149‑1 (MQ=255)                                                                                                 ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCg                                                   <  8:374961/149‑1 (MQ=255)                                                                                                      accagGCGCCCGTCCGCCGCGCCATGGTTGCGGTAGGGTAACATCGACATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGGCCTGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                              <  1:124515/149‑1 (MQ=255)                                                                                                         acGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTg                                           <  1:206955/149‑1 (MQ=255)                                                                                                         acGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTg                                           >  1:282897/1‑149 (MQ=255)                                                                                                           gcgcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTgg                                          >  1:331922/1‑149 (MQ=255)                                                                                                           gcgcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGgc                                         <  7:147149/149‑1 (MQ=255)                                                                                                                cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                              >  8:202184/1‑139 (MQ=255)                                                                                                                cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                              <  7:202184/139‑1 (MQ=255)                                                                                                                cGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGgccgcc                                     >  5:438029/1‑148 (MQ=255)                                                                                                                        gcgcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCatat                             >  8:313607/1‑148 (MQ=255)                                                                                                                        gcgcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCatat                             >  1:169744/1‑148 (MQ=255)                                                                                                                         cgcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGGTCCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGa                           <  1:309977/149‑1 (MQ=255)                                                                                                                          gcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGAc                          <  7:397162/149‑1 (MQ=255)                                                                                                                              tGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAat                      <  2:169744/149‑1 (MQ=255)                                                                                                                                gTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAatat                    >  4:216183/1‑149 (MQ=255)                                                                                                                                    cGGTAGTGGAACATCGAGATGCTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAACCGCAGCACCATTTCGTCGCCGGCACCGGCCGAGTGGCCGCCGACCATATTACGAATATGCAg                >  5:232438/1‑149 (MQ=255)                                                                                                                                        aGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACGATATGACGAATATGCAGCTTg            >  6:412053/1‑149 (MQ=255)                                                                                                                                          tGGAACATCTCGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGCACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGc          >  1:242523/1‑149 (MQ=255)                                                                                                                                          tGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGc          >  7:129022/1‑149 (MQ=255)                                                                                                                                            gAACATCTCGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGcc         >  4:335522/1‑148 (MQ=255)                                                                                                                                             aaCATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCaa       >  2:103866/1‑149 (MQ=255)                                                                                                                                                aTCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCTGCCGAGTGGCCGCCGACCATATGACGAATATGc                  <  3:81902/135‑1 (MQ=255)                                                                                                                                                aTCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTc    >  2:206827/1‑149 (MQ=255)                                                                                                                                                  cGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATTCAGCTTGGCCAATTCGt  >  3:297369/1‑149 (MQ=255)                                                                                                                                                   gagaTGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGt  <  6:438029/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGC‑GCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGT  >  NC_002947/5874673‑5874966

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