Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,551,282	G→T	9.0%	intergenic (+41/+85)	groL → / ← pykA	chaperonin GroEL/pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,551,282	0	G	T	9.0%	197.8 / 9.3	67	intergenic (+41/+85)	groL/pykA	chaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base T (3/3);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

GTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551135‑1551430                                                                                                                                                     |                                                                                                                                                     gTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                     <  6:191255/149‑1 (MQ=255)    tgttgatCACCACCGAAGCCATGGTTGCCGTCCTTCCGGAAGTCAAGCCAGCTGCTGGCATGGCTGATATGGGTGGCATGGGTGGGATGGGCGGCATGATGTACAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCga                                                                                                                                                   <  1:84824/146‑1 (MQ=255)    tgatgaTCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                     >  7:13486/1‑146 (MQ=255)     gatgatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCgag                                                                                                                                                  >  4:56371/1‑148 (MQ=255)      atgatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCgag                                                                                                                                                  <  8:158965/147‑1 (MQ=255)      atgatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGaa                                                                                                                                                >  3:164235/1‑149 (MQ=255)        gatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAAcc                                                                                                                                              >  4:95018/1‑149 (MQ=255)         atCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCg                                                                                                                                             >  8:92178/1‑149 (MQ=255)          tCACCACCGAAGCCATGGTTGCCGACCTGTCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCg                                                                                                                                             <  7:342523/148‑1 (MQ=255)           caccacCGAAGCCATGGTTGCCGACCTACCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCggg                                                                                                                                           >  8:303529/1‑149 (MQ=255)              caccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCt                                                                                                                                        >  6:351247/1‑149 (MQ=255)               accGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtt                                                                                                                                       >  4:9569/1‑149 (MQ=255)                 cGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtttt                                                                                                                                     >  7:319275/1‑149 (MQ=255)                        aTGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGc                                                                                                                              >  8:55684/1‑149 (MQ=255)                        aTGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTGTTTATTGc                                                                                                                              >  3:45404/1‑149 (MQ=255)                                cGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCgggg                                                                                                                      >  3:310013/1‑149 (MQ=255)                                 gACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGc                                                                                                                     >  7:361274/1‑149 (MQ=255)                                 gACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGc                                                                                                                     <  8:13486/149‑1 (MQ=255)                                             gACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgcagca                                                                                                         >  6:21908/1‑149 (MQ=255)                                              aCAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        <  5:351247/149‑1 (MQ=255)                                                aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCggg                                                                                                                                           >  4:392277/1‑112 (MQ=255)                                                aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCggg                                                                                                                                           >  4:398614/1‑112 (MQ=255)                                                aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCggg                                                                                                                                           <  3:398614/112‑1 (MQ=255)                                                aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCggg                                                                                                                                           <  3:392277/112‑1 (MQ=255)                                                    cAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCt                                                                                                  >  5:253095/1‑149 (MQ=255)                                                        tgctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTgcgc                                                                                              >  1:20858/1‑149 (MQ=255)                                                             gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCGCCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTACTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  3:276528/148‑1 (MQ=255)                                                             gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  5:353931/148‑1 (MQ=255)                                                                   cTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtgt                                                                                   <  2:20858/149‑1 (MQ=255)                                                                      atatGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtgtg                                                                                  <  5:21908/147‑1 (MQ=255)                                                                       tatGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                               <  6:253095/149‑1 (MQ=255)                                                                        atgggtggcatgggtggcATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAATAAGCCCCGGTTCGCGAGAACCGGGGCTTGTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtgt                                                                              >  5:231797/1‑149 (MQ=255)                                                                             tggcatgggtggcatgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         <  4:164235/149‑1 (MQ=255)                                                                             tggcatgggtggcatgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         <  8:319275/149‑1 (MQ=255)                                                                             tggcatgggggccaggggCGCCATGATGAA‑CCCCAGCATTACCCCCCTCCCCCACCAAAAAGCCCCGGTGCGCGATAACCGGGGGTTTTTTATTGCCGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         <  1:106668/149‑1 (MQ=255)                                                                                catgggtggcatgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAg                                                                      <  4:310013/149‑1 (MQ=255)                                                                                  tgggtggcatgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         <  4:100597/144‑1 (MQ=255)                                                                                  tgggtggcatgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         >  3:100597/1‑144 (MQ=255)                                                                                  tggggggcatgggCGGCATTATGTA‑AGCCATCCTTACCCCCCTGCACCATCAAAAATCCCCGGTTCGCGACAACCGGGGCTTTTTTATTTCTGTCGGGGCCTCTTCGCAGCACAAGGCt                                                                                                  <  6:244272/119‑1 (MQ=255)                                                                                     gtggcatgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTGGTGCTGTGAAACGAGTGGAt                                                                 <  3:95018/149‑1 (MQ=255)                                                                                         catgggCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAc                                                             >  4:339669/1‑149 (MQ=255)                                                                                                ggcATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  4:216960/1‑134 (MQ=255)                                                                                                ggcATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  3:216960/134‑1 (MQ=255)                                                                                                 gcATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGc                                                     <  2:295962/149‑1 (MQ=255)                                                                                                     gatgTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGAGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                                 >  3:320081/1‑149 (MQ=255)                                                                                                       tgTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCATTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTAGTGCTGTGAAACGACTGGATCAACCGACCAGCGGATCa                                               >  4:131657/1‑149 (MQ=255)                                                                                                        gTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAc                                              >  1:281406/1‑149 (MQ=255)                                                                                                        gTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAc                                              >  8:51728/1‑149 (MQ=255)                                                                                                         tA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAc                                              <  7:51728/148‑1 (MQ=255)                                                                                                          a‑aGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCtt                                                                                            <  5:331213/101‑1 (MQ=255)                                                                                                          a‑aGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCtt                                                                                            >  6:331213/1‑101 (MQ=255)                                                                                                            agccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          >  7:43432/1‑102 (MQ=255)                                                                                                            agccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  8:43432/102‑1 (MQ=255)                                                                                                                        cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCt                               >  5:347045/1‑149 (MQ=255)                                                                                                                        cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCt                               <  3:9569/149‑1 (MQ=255)                                                                                                                          cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGAGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  4:320081/149‑1 (MQ=255)                                                                                                                           ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgc                                                                                                             >  3:153109/1‑68 (MQ=255)                                                                                                                           ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgc                                                                                                             <  4:153109/68‑1 (MQ=255)                                                                                                                               cACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATg                          <  8:361274/147‑1 (MQ=255)                                                                                                                                   aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTg                    >  3:161037/1‑149 (MQ=255)                                                                                                                                      aaaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTTCGCTTGAGCCTTGTGCTGTGAAACGAGTGTATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                  <  3:131657/148‑1 (MQ=255)                                                                                                                                       aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                  <  2:80794/147‑1 (MQ=255)                                                                                                                                       aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                  >  1:80794/1‑147 (MQ=255)                                                                                                                                           gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAt                      <  5:185820/139‑1 (MQ=255)                                                                                                                                           gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAt                      >  6:185820/1‑139 (MQ=255)                                                                                                                                           gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            >  7:29260/1‑149 (MQ=255)                                                                                                                                              ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGAt         >  7:28243/1‑149 (MQ=255)                                                                                                                                              ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGAt         >  8:60708/1‑149 (MQ=255)                                                                                                                                              ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGAt         >  5:381932/1‑149 (MQ=255)                                                                                                                                                    cgcgAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAATGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtgt    <  2:281406/148‑1 (MQ=255)                                                                                                                                                     gcgAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  >  8:151211/1‑149 (MQ=255)                                                                                                                                                     gcgAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  <  7:151211/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTA‑AGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551135‑1551430

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑