Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,288,779	G→T	5.6%	A18S (GCC→TCC)	PP_1126 →	hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,288,779	0	G	T	5.6%	291.1 / 5.7	89	A18S (GCC→TCC)	PP_1126	hydrolase
Reads supporting (aligned to +/- strand):  ref base G (35/49);  new base T (2/3);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

CCTGCGGCCGCCTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCTTCGT  >  NC_002947/1288632‑1288927                                                                                                                                                    |                                                                                                                                                     ccTGCGGCCGCCTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTgcgc                                                                                                                                                     >  3:397743/1‑149 (MQ=255)  cTGCGGCCGCCTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGcc                                                                                                                                                    >  3:335822/1‑149 (MQ=255)     cGGCCGCCTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGCGCAAGTGCTCCGgc                                                                                                                                                 <  3:16411/149‑1 (MQ=255)     cGGCCGCCTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGgc                                                                                                                                                 <  5:42142/149‑1 (MQ=255)        ccgccTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCa                                                                                                                                              <  4:330574/149‑1 (MQ=255)        ccgccTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCa                                                                                                                                              >  3:330574/1‑149 (MQ=255)          gccTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCTCCGTTCCGGCGGGAAATTTCCTGATCGCCGGCGCAAGTGCTCCGgcgatca                                                                                                                                            <  8:51157/149‑5 (MQ=255)          gccTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGCGCAAGTGCTCCGgcgatca                                                                                                                                            <  5:259919/149‑5 (MQ=255)            cTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGt                                                                                                                                          >  4:329324/1‑149 (MQ=255)                  cGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTa                                                                                                                                    <  2:198679/149‑1 (MQ=255)                       cGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGCCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGt                                                                                                                               >  2:53082/1‑149 (MQ=255)                        gAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCTCGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTc                                                                                                                              >  1:306378/1‑149 (MQ=255)                        gAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGATACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTc                                                                                                                              <  6:130753/149‑1 (MQ=255)                           cAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTg                                                                                                                           >  7:320168/1‑149 (MQ=255)                            aGGTTCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTgg                                                                                                                          <  1:122906/149‑1 (MQ=255)                            aGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTgg                                                                                                                          <  7:368780/149‑1 (MQ=255)                             ggATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTggg                                                                                                                         >  3:36905/1‑149 (MQ=255)                               aTCGCGGCGTTCTCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGcc                                                                                                                       <  2:405331/149‑1 (MQ=255)                                  gcggcgTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCAt                                                                                                                    >  4:126721/1‑149 (MQ=255)                                   cggcgTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCAt                                                                                                                    >  7:78370/1‑148 (MQ=255)                                    ggcgTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCa                                                                                                                  <  3:258308/149‑1 (MQ=255)                                    ggcgTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCa                                                                                                                  <  3:126721/149‑1 (MQ=255)                                      cgTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAgg                                                                                                                <  2:111961/149‑1 (MQ=255)                                       gTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAggg                                                                                                               <  4:399600/149‑1 (MQ=255)                                       gTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAggg                                                                                                               <  8:73724/149‑1 (MQ=255)                                          aTCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTcc                                                                                                            >  6:333709/1‑149 (MQ=255)                                          aTCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTc                                                                                                             <  5:114242/148‑1 (MQ=255)                                             aCCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTg                                                                                                          >  7:181561/1‑148 (MQ=255)                                              ccTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGaa                                                                                                        <  7:87902/149‑1 (MQ=255)                                               cttgcttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAg                                                                                                       <  2:204565/149‑1 (MQ=255)                                               cttgcttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAg                                                                                                       >  5:355100/1‑149 (MQ=255)                                               cttgcttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAg                                                                                                       <  8:95286/149‑1 (MQ=255)                                                   cttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCGATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGga                                                                                                     >  5:77036/1‑147 (MQ=255)                                                   cttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGgag                                                                                                    <  2:113364/148‑1 (MQ=255)                                                   cttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGgag                                                                                                    <  6:152479/148‑1 (MQ=255)                                                    ttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGagag                                                                                                  >  4:154589/1‑149 (MQ=255)                                                    ttgACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGagag                                                                                                  >  6:378263/1‑149 (MQ=255)                                                      gACCGGCGCTTGGTCACGCAGCGCTACGGCAAGTCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGc                                                                                                 <  1:53082/148‑1 (MQ=255)                                                      gACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTGAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGATCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGc                                                                                                 <  4:385133/148‑1 (MQ=255)                                                      gACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGc                                                                                                 <  4:335822/148‑1 (MQ=255)                                                        ccGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCt                                                                                               >  1:428740/1‑148 (MQ=255)                                                         cGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCt                                                                                               >  1:210259/1‑147 (MQ=255)                                                         cGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCt                                                                                               >  1:335800/1‑147 (MQ=255)                                                         cGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCa                                                                                             <  8:136979/149‑1 (MQ=255)                                                         cGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCa                                                                                             <  6:259611/149‑1 (MQ=255)                                                           gcgcTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGc                                                                                           >  4:311958/1‑149 (MQ=255)                                                             gcTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCTGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCt                                                                                         >  6:186607/1‑149 (MQ=255)                                                                tGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCTGGCGCATTGTCCATTACTCGACTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTAt                                                                                      >  5:375553/1‑149 (MQ=255)                                                                 ggTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTa                                                                                       <  2:428740/147‑1 (MQ=255)                                                                 ggTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTa                                                                                       <  5:136904/147‑1 (MQ=255)                                                                  gTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGa                                                                                    <  4:22444/149‑1 (MQ=255)                                                                     aCGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTgg                                                                                                                          <  1:24794/108‑1 (MQ=255)                                                                     aCGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTgg                                                                                                                          >  2:24794/1‑108 (MQ=255)                                                                       gCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATccgc                                                                               >  4:11215/1‑149 (MQ=255)                                                                             cTACGGCAAGGCTATTCTCGATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCAcc                                                                         <  6:77036/149‑1 (MQ=255)                                                                                cGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCt                                                                       <  8:181561/148‑1 (MQ=255)                                                                                cGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTACTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTg                                                                      <  3:329324/149‑1 (MQ=255)                                                                                  gCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGcc                                                                    <  2:268958/149‑1 (MQ=255)                                                                                      ggCTATTCTCAATGCGCTGCCGCTGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCccgcc                                                                <  5:186607/149‑1 (MQ=255)                                                                                        gtATGCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGt                                                              <  4:273670/148‑1 (MQ=255)                                                                                           ttCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCTGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTa                                                             <  7:300782/147‑1 (MQ=255)                                                                                                aaTGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCt                                                      >  8:244257/1‑149 (MQ=255)                                                                                                 aTGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCtt                                                     <  2:88172/149‑1 (MQ=255)                                                                                                 aTGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGACGCACCCTGCCCGCCGTATTCGCCtt                                                     <  3:154589/149‑1 (MQ=255)                                                                                                  tGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCttg                                                    >  1:42530/1‑149 (MQ=255)                                                                                                    cgcTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCttgtt                                                  <  5:18334/149‑1 (MQ=255)                                                                                                         ccgccgTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCa                                             >  4:286520/1‑149 (MQ=255)                                                                                                            ccgTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGtt                                          <  3:286520/149‑1 (MQ=255)                                                                                                             cgTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTc                                         <  4:417683/149‑1 (MQ=255)                                                                                                                      gggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTAGCCTTGTTGTGCAGTTCGCCCTTGCt                                >  4:360113/1‑149 (MQ=255)                                                                                                                       ggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCccgcc                                                                >  1:300299/1‑116 (MQ=255)                                                                                                                       ggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCccgcc                                                                <  2:300299/116‑1 (MQ=255)                                                                                                                       ggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTg                               <  2:306378/149‑1 (MQ=255)                                                                                                                       ggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTg                               <  8:320168/149‑1 (MQ=255)                                                                                                                       ggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTg                               >  3:379911/1‑149 (MQ=255)                                                                                                                       ggAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTg                               >  3:402866/1‑149 (MQ=255)                                                                                                                         aaaTTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGc                             >  4:335979/1‑149 (MQ=255)                                                                                                                              tCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGgc                        >  7:326965/1‑149 (MQ=255)                                                                                                                               ccTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGgc                        <  4:36905/148‑1 (MQ=255)                                                                                                                               ccTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGgc                        <  5:333709/148‑1 (MQ=255)                                                                                                                                cTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGgca                       <  4:397743/148‑1 (MQ=255)                                                                                                                                   aTCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGgcagca                    <  6:355100/148‑1 (MQ=255)                                                                                                                                      gCCGGCGCAAGTGCTCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGGCCTTGCTGGCCGGGCAGCAAACg                >  5:74554/1‑149 (MQ=255)                                                                                                                                      gCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACg                <  4:363574/149‑1 (MQ=255)                                                                                                                                       ccGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCCCCCTGACCGCCGGATTAGCCTTGTTGTGAAGTTAGCCCTTGCTGTCCGGGCAGCAAACGc               >  8:384413/1‑149 (MQ=255)                                                                                                                                        cGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCt              <  5:378263/149‑1 (MQ=255)                                                                                                                                         ggAGCACTTGCGCTGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTg             <  6:375553/149‑1 (MQ=255)                                                                                                                                          gAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCATCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCCGCAAACggggg            >  6:294026/1‑145 (MQ=255)                                                                                                                                            gcaagtGCTCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGGCTTATGATCCGCCGCACCCTGCCCGACGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTc          >  2:157004/6‑149 (MQ=255)                                                                                                                                                ttGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCt      >  3:182989/1‑149 (MQ=255)                                                                                                                                                   cgcCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGGCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCTTCg   <  7:384413/149‑1 (MQ=255)                                                                                                                                                    gcCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCTTCGt  <  4:379911/149‑1 (MQ=255)                                                                                                                                                    gcCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCTTCGt  <  3:311958/149‑1 (MQ=255)                                                                                                                                                    gcCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCTTCGt  >  2:300613/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCTGCGGCCGCCTTCTGCGTACCGAGCAGGATCGCGGCGTTATCACCTTGCTTGACCGGCGCTTGGTCACGCAGCGCTACGGCAAGGCTATTCTCAATGCGCTGCCGCCGTTCCGGCGGGAAATTTCCTGATCGCCGGAGCACTTGCGCCGGCGCATTGTCCATTACTCGTCTTGGGCCCATCAGGGTCCTGAAGGAGAGCCTCAGCCTTATGATCCGCCGCACCCTGCCCGCCGTATTCGCCTTGTTGTGCAGTTCGCCCTTGCTGGCCGGGCAGCAAACGCTGTTCAGCTTCGT  >  NC_002947/1288632‑1288927

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