Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,809,179	A→C	7.4%	D77A (GAC→GCC)	eno →	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,809,179	0	A	C	7.4%	204.1 / 6.2	68	D77A (GAC→GCC)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base C (2/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GACAACGGCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACC  >  NC_002947/1809031‑1809325                                                                                                                                                     |                                                                                                                                                   gACAACGGCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTAGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGa                                                                                                                                                    >  8:150927/1‑149 (MQ=255)       ggCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCCACGGCCCGATCCGTGACctgct                                                                                                                                              >  4:240827/1‑149 (MQ=255)             atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAggt                                                                                                                                        >  6:219311/1‑148 (MQ=255)                ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGc                                                                                                                                       >  5:290985/1‑147 (MQ=255)                 gCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCacg                                                                                                                                    <  5:219311/149‑3 (MQ=255)                     cgcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATc                                                                                                                                <  1:167548/149‑1 (MQ=255)                     cgcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATc                                                                                                                                >  3:81773/1‑149 (MQ=255)                     cgcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATc                                                                                                                                <  1:271889/149‑1 (MQ=255)                       cTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                              <  4:81773/149‑3 (MQ=255)                        ttGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCtt                                                                                                                             >  7:187082/1‑149 (MQ=255)                          gcgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCtt                                                                                                                             <  3:34123/147‑1 (MQ=255)                                    cGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAgg                                                                                                                 >  7:172819/1‑149 (MQ=255)                                      gTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACctgctgctg                                                                                                                                          >  4:119683/1‑122 (MQ=255)                                       tGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGctc                                                                                                              <  8:187082/149‑1 (MQ=255)                                         cTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTg                                                                                                            <  6:290985/149‑1 (MQ=255)                                         cTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTg                                                                                                            >  4:174037/1‑149 (MQ=255)                                         cTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGAAGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCCTGACCTGCTGCTGGGCAAGGATCCTTCCGAACAGAAGGCTCTg                                                                                                            >  4:248849/1‑149 (MQ=255)                                            cTACCGTCTCGCGCGAAGCGCTGGAGCTGCGTTATGGCGAAAAGAGCCTTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCTATCCGTGACCTGCTGCTTGGCAAGGATCCTTCCGATCAGAATGctcttg                                                                                                           <  4:154196/147‑3 (MQ=255)                                              aCCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcg                                                                                                       >  2:96100/1‑149 (MQ=255)                                              aCCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcg                                                                                                       >  3:79111/1‑149 (MQ=255)                                               ccGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGAc                                                                                                         <  3:174037/146‑1 (MQ=255)                                               ccGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcgc                                                                                                      <  4:79111/149‑1 (MQ=255)                                                      gcgcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGa                                                                                               >  5:135565/1‑149 (MQ=255)                                                       cgcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGTCCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGAt                                                                                              <  1:202149/149‑1 (MQ=255)                                                        gcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCa                                                                                                  <  6:355967/144‑1 (MQ=255)                                                        gcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCa                                                                                                  >  5:355967/1‑144 (MQ=255)                                                         cgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATc                                                                                             >  3:310027/1‑148 (MQ=255)                                                            aGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                         >  1:71712/1‑149 (MQ=255)                                                              cgcTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTg                                                                                       <  5:148704/149‑1 (MQ=255)                                                                   gAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCAGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACgg                                                                                  >  6:66752/1‑149 (MQ=255)                                                                          gTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCgag                                                                           <  4:310027/149‑1 (MQ=255)                                                                                  gACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGc                                                                   <  7:150927/149‑1 (MQ=255)                                                                                    cAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGcc                                                                  >  1:9761/1‑148 (MQ=255)                                                                                     aaGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGcc                                                                  >  5:123002/1‑147 (MQ=255)                                                                                       gagCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCg                                                                             <  3:49355/134‑1 (MQ=255)                                                                                       gagCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCg                                                                             >  4:49355/1‑134 (MQ=255)                                                                                        agCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAg                                                                     <  1:145779/141‑1 (MQ=255)                                                                                        agCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAg                                                                     >  2:145779/1‑141 (MQ=255)                                                                                          ccGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                         >  6:80331/1‑119 (MQ=255)                                                                                          ccGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                         <  5:80331/119‑1 (MQ=255)                                                                                                    ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                         >  6:349688/1‑141 (MQ=255)                                                                                                    ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                         <  5:349688/141‑1 (MQ=255)                                                                                                    ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACg                                                  <  2:9761/148‑1 (MQ=255)                                                                                                         gggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAg                                                               >  3:220213/1‑130 (MQ=255)                                                                                                         gggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAg                                                               <  4:220213/130‑1 (MQ=255)                                                                                                          ggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGTCGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCt                                           <  6:123002/149‑1 (MQ=255)                                                                                                              tGCTGAAGGCCGTCGCCAACATTAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCt                                       >  7:283096/1‑149 (MQ=255)                                                                                                                 tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGc                                                                   >  4:142235/1‑118 (MQ=255)                                                                                                                 tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    >  7:84797/1‑149 (MQ=255)                                                                                                                 tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    >  1:23058/1‑149 (MQ=255)                                                                                                                 tGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    >  1:254770/1‑149 (MQ=255)                                                                                                                   aaGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGtt                                    <  5:66752/147‑1 (MQ=255)                                                                                                                       ccGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCt                               <  2:193804/148‑1 (MQ=255)                                                                                                                       ccGTCGCCAACATCAACGGCCCGATCCGAGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCt                               >  1:193804/1‑148 (MQ=255)                                                                                                                        cGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTgg                             <  3:240827/149‑1 (MQ=255)                                                                                                                         gTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGc                            <  8:84797/149‑1 (MQ=255)                                                                                                                         gTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGc                            <  2:254770/149‑1 (MQ=255)                                                                                                                               aaCATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTgccgc                      >  5:229198/1‑149 (MQ=255)                                                                                                                                 cATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTg                                          >  6:144776/1‑127 (MQ=255)                                                                                                                                 cATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTg                                          <  5:144776/127‑1 (MQ=255)                                                                                                                                  aTCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCa                    >  1:326293/1‑148 (MQ=255)                                                                                                                                     aaCGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAg                 <  8:283096/148‑1 (MQ=255)                                                                                                                                              aTCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCcca       <  6:229198/149‑1 (MQ=255)                                                                                                                                                ccGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGccc        <  2:83451/146‑1 (MQ=255)                                                                                                                                                ccGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGccc        >  1:83451/1‑146 (MQ=255)                                                                                                                                                 cGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccagg     <  3:248849/148‑1 (MQ=255)                                                                                                                                                   tGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggacc  <  6:135565/149‑1 (MQ=255)                                                                                                                                                   tGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggac   >  4:15695/1‑148 (MQ=255)                                                                                                                                                    gACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                         <  2:231942/93‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GACAACGGCATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACC  >  NC_002947/1809031‑1809325

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