Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,682,856	T→A	6.2%	L459M (TTG→ATG)	recJ →	ssDNA‑specific exonuclease RecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,682,856	0	T	A	6.2%	188.1 / 2.6	65	L459M (TTG→ATG)	recJ	ssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base A (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CGCCACCCGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCA  >  NC_002947/1682709‑1683004                                                                                                                                                    |                                                                                                                                                     cgcCACCCGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCtgtt                                                                                                                                                     <  4:27997/149‑1 (MQ=255) cgcCACCCGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCtgtt                                                                                                                                                     <  8:196188/149‑1 (MQ=255)    cACCCGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTc                                                                                                                                                  >  3:204164/1‑149 (MQ=255)    cACCCGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTc                                                                                                                                                  >  7:163639/1‑149 (MQ=255)        cGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGAc                                                                                                                                              <  4:31436/149‑1 (MQ=255)         gCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACg                                                                                                                                             >  6:166825/1‑149 (MQ=255)            gCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACgg                                                                                                                                            >  1:161294/1‑147 (MQ=255)             cTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAg                                                                                                                                         >  1:183062/1‑149 (MQ=255)             cTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAg                                                                                                                                         >  6:90189/1‑149 (MQ=255)                aTCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCt                                                                                                                                      >  7:208684/1‑149 (MQ=255)                    gCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCg                                                                                                                                  >  5:349016/1‑149 (MQ=255)                    gCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCg                                                                                                                                  <  5:260254/149‑1 (MQ=255)                     cAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCgg                                                                                                                                 >  5:85565/1‑149 (MQ=255)                       aGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTa                                                                                                                               <  1:284165/149‑1 (MQ=255)                        gTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCt                                                                                                                                      <  4:6147/141‑1 (MQ=255)                        gTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCt                                                                                                                                      >  3:6147/1‑141 (MQ=255)                         ttCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTaga                                                                                                                             <  7:72272/149‑1 (MQ=255)                            ggcggcCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAgagga                                                                                                                          <  2:161294/149‑1 (MQ=255)                            ggcggcCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAgagga                                                                                                                          >  4:13477/1‑149 (MQ=255)                                           gCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCt                                                                                                           <  3:13477/149‑1 (MQ=255)                                                 gggTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGcc                                                                                                       >  7:106724/1‑147 (MQ=255)                                                  ggTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCATTCGACGAAGAGGTGCGACGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAg                                                                                                    >  1:72038/1‑149 (MQ=255)                                                       tCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCt                                                                                               <  5:274832/149‑1 (MQ=255)                                                         gTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTg                                                                                              <  6:349016/148‑1 (MQ=255)                                                             ccAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGt                                                                                                                        <  4:249451/118‑1 (MQ=255)                                                             ccAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGt                                                                                                                        >  3:249451/1‑118 (MQ=255)                                                                   gCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCgg                                                                                   >  2:212352/1‑149 (MQ=255)                                                                    cTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGc                                                                                  >  5:322285/1‑149 (MQ=255)                                                                     tAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGcc                                                                                 >  7:216524/1‑149 (MQ=255)                                                                     tAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGcc                                                                                 >  2:216140/1‑149 (MQ=255)                                                                       aTTTCTCGGCGTTTGCCGAGGCGTGCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGcc                                                                                 <  7:205734/147‑1 (MQ=255)                                                                           ctcGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATgg                                                                            <  7:301250/148‑1 (MQ=255)                                                                                cGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAg                                                                      <  6:85565/149‑1 (MQ=255)                                                                                     gCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAGGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCCATGGGCCGGCatg                                                                  <  2:228855/148‑3 (MQ=255)                                                                                       cGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTcc                                                               >  4:239221/1‑149 (MQ=255)                                                                                       cGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTcc                                                               >  5:109659/1‑149 (MQ=255)                                                                                       cGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTCGCCAAGGCCCTGCGCCATGACGGCCCATGGGGCCAGCACTTcc                                                               >  6:115818/1‑149 (MQ=255)                                                                                         aGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCt                                                             >  8:219678/1‑149 (MQ=255)                                                                                         aGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCt                                                             <  1:212352/149‑1 (MQ=255)                                                                                          ggCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTg                                                            <  7:219678/149‑1 (MQ=255)                                                                                           gCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGa                                                           <  5:166825/149‑1 (MQ=255)                                                                                           gCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGa                                                           <  8:106724/149‑1 (MQ=255)                                                                                                 gacgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCt                                                     <  8:208684/149‑1 (MQ=255)                                                                                                   cgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGt                                                   <  2:183062/149‑1 (MQ=255)                                                                                                      agagGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTc                                                 >  7:342839/1‑148 (MQ=255)                                                                                                        agGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCAt                                              <  2:133599/149‑1 (MQ=255)                                                                                                           tGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGc                                           <  8:163639/149‑1 (MQ=255)                                                                                                                  tCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTc                                     >  5:375057/1‑148 (MQ=255)                                                                                                                      cTGCGCGAAGAAGACCTGACTGGCCGCCTGATGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGtt                                >  4:328496/1‑149 (MQ=255)                                                                                                                      cTGCGCGAAGAAGACCTGACTGGCCGCCTGATGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGtt                                >  4:333848/1‑149 (MQ=255)                                                                                                                         cgcgAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGGCGCTGTTCCATGGCGTGTTCCAGTtggt                             <  1:216140/149‑1 (MQ=255)                                                                                                                          gcgAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTtggtg                            <  5:90189/149‑1 (MQ=255)                                                                                                                            gaagaaGACGTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGa                          <  4:204164/149‑1 (MQ=255)                                                                                                                            gaagaaGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGa                          <  6:322285/149‑1 (MQ=255)                                                                                                                                aagaCCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCa                                   <  7:210709/136‑1 (MQ=255)                                                                                                                                aagaCCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCa                                   >  8:210709/1‑136 (MQ=255)                                                                                                                                 agaCCTGACTGGCCGCCTGATGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGagcagc                     <  3:333848/149‑1 (MQ=255)                                                                                                                                 agaCCTGACTGGCCGCCTGATGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGagcagc                     <  3:328496/149‑1 (MQ=255)                                                                                                                                 agaACTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGagcagc                     <  8:342839/149‑1 (MQ=255)                                                                                                                                    ccTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCgtg                  <  8:216524/149‑1 (MQ=255)                                                                                                                                      tGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGtgtg                >  8:195036/1‑149 (MQ=255)                                                                                                                                      tGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGtgtg                >  6:363003/1‑149 (MQ=255)                                                                                                                                          tGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGa                          <  8:137734/135‑1 (MQ=255)                                                                                                                                          tGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGa                          >  7:137734/1‑135 (MQ=255)                                                                                                                                            gccgccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTg            <  2:72038/147‑1 (MQ=255)                                                                                                                                            gccgccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGc          >  1:95464/1‑149 (MQ=255)                                                                                                                                            gccgccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTCCCACCGCGACCTGGCCAAGGCCCTGCGCCCTGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTg            <  4:222422/147‑1 (MQ=255)                                                                                                                                                    ttgtCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCa  <  2:95464/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGCCACCCGCAGCTGATCAGCAAGTTCGGCGGCCACGCCATGGCTGCCGGGTTGTCGTTGCCAGAGGCTAATTTCTCGGCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCCATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCA  >  NC_002947/1682709‑1683004

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