Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,551,282	G→T	8.2%	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	8.2%	180.6 / 3.1	59	intergenic (+41/+85)	groL/pykA	chaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base T (2/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551135‑1551430                                                                                                                                                    |                                                                                                                                                      gTCTGATGATCACCACCGAAGCCTTGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGGGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                      <  6:184409/149‑1 (MQ=255)          tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCgg                                                                                                                                            >  3:80931/1‑149 (MQ=255)           caccacCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCggg                                                                                                                                           >  2:106922/1‑149 (MQ=255)            accaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGC‑GAACCgggg                                                                                                                                          <  7:335731/149‑1 (MQ=255)                  gAAGCCATGGTTGCCGCCCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCttttt                                                                                                                                    >  5:107134/1‑149 (MQ=255)                      ccATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTAtt                                                                                                                                >  6:143832/1‑149 (MQ=255)                       cATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTg                                                                                                                               >  8:171027/1‑149 (MQ=255)                         tGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGCGGCATGGGGGGCATGGGCGGCATGATGTAAGCCCGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCt                                                                                                                             <  6:309447/149‑1 (MQ=255)                               ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                       <  7:193828/149‑1 (MQ=255)                                     tGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCctat                                                                                                                 >  3:245991/1‑147 (MQ=255)                                       ccGGAAGACAATCAAGCTGCTGGCATGCCTGATATGGGGTGCATGCGGCGCATGGGCGGCATGATGTAAGCCAGCCTTACACCCCTGCACCATCAAAAAGCCCCGGTTCGCGC‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCctct                                                                                                                 <  4:213142/147‑1 (MQ=255)                                          gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                            <  4:29870/149‑1 (MQ=255)                                                   ccAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGc                                                                                                   >  2:203970/1‑149 (MQ=255)                                                   ccAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGc                                                                                                   <  5:143832/149‑1 (MQ=255)                                                     aGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGC‑GAACCgg                                                                                                                                            >  3:125074/1‑106 (MQ=255)                                                     aGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGC‑GAACCgg                                                                                                                                            <  4:125074/106‑1 (MQ=255)                                                            ggaaTGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  7:270344/146‑1 (MQ=255)                                                             gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGATCTACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  1:106922/149‑1 (MQ=255)                                                                tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCt                                                                                      >  1:354320/1‑149 (MQ=255)                                                                    tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                       >  3:234616/1‑112 (MQ=255)                                                                    tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                       <  4:234616/112‑1 (MQ=255)                                                                       tatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                               <  5:353085/149‑1 (MQ=255)                                                                       tatGGGTGGCATGGGAGGCATGCGCGGCATGAAGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCCGTTCGCGA‑CAACCGGGGCTGTTTTATAGCTCTCGGGGCCTCTTCGCACCACAAGGCTGCGCTTGCGCCTAGTGCtg                                                                               >  6:309075/1‑149 (MQ=255)                                                                            gtggcatgggtggcatggGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGaaa                                                                          <  8:238098/149‑1 (MQ=255)                                                                             tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          >  5:384164/1‑132 (MQ=255)                                                                             tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  6:384164/132‑1 (MQ=255)                                                                               gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         >  4:44734/1‑147 (MQ=255)                                                                                 atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  6:107134/149‑1 (MQ=255)                                                                                 atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  3:15462/149‑1 (MQ=255)                                                                                    ggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                   <  7:171027/148‑1 (MQ=255)                                                                                          atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAcc                                                            >  4:11258/1‑149 (MQ=255)                                                                                             ggcggcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGAc                                                         <  3:11258/149‑1 (MQ=255)                                                                                              gcggcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGAc                                                         >  6:359397/1‑148 (MQ=255)                                                                                                 gcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGc                                                     <  1:343234/149‑1 (MQ=255)                                                                                                     gatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGAGCAACCGACCAGCGGAt                                                 >  1:271962/1‑149 (MQ=255)                                                                                                      atgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                                  >  6:311451/1‑131 (MQ=255)                                                                                                      atgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                                  <  5:311451/131‑1 (MQ=255)                                                                                                         tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAcc                                             >  4:31387/1‑149 (MQ=255)                                                                                                          aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATc                                                                <  7:171353/129‑1 (MQ=255)                                                                                                          aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATc                                                                >  8:171353/1‑129 (MQ=255)                                                                                                           agccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGGCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGa                                           <  2:354320/149‑1 (MQ=255)                                                                                                            gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                          >  4:123518/1‑149 (MQ=255)                                                                                                            gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                          >  8:202925/1‑149 (MQ=255)                                                                                                            gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                          >  5:137097/1‑149 (MQ=255)                                                                                                             ccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGCCCAGCGGATCACCGACg                                         >  8:160462/1‑149 (MQ=255)                                                                                                                  cTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAg                                    <  3:31387/149‑1 (MQ=255)                                                                                                                         cccTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                                 <  4:60592/129‑1 (MQ=255)                                                                                                                         cccTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                                 >  3:60592/1‑129 (MQ=255)                                                                                                                         cccTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  6:137097/149‑1 (MQ=255)                                                                                                                           cTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCATGTGCTGTGAAACGAGAGGATCAACCGACCAGCGGATCACCGAAGTGCAGGAACTTCAt                           >  8:57035/1‑149 (MQ=255)                                                                                                                               aCCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCa                       >  3:354551/1‑149 (MQ=255)                                                                                                                               aCCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGAGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGCCCAGCGGATCACCGACGTGCAGGATCTTCATGCCa                       >  6:55255/1‑149 (MQ=255)                                                                                                                                  aTCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTg                    <  7:160462/149‑1 (MQ=255)                                                                                                                                      aaaaGCCCCGGATCTCGC‑GACCGGGGGCTTTTTTATTGATGTCGGGGCCTGTTCCCAGCACAAGGGTGCGCTTGAGCGTAGTGCTGTGAAACGAGTGGATCACCGGACGAGCGGATCACCGCCGTGCAGGATGTTCATGCCATTGGtcc                >  8:122467/1‑147 (MQ=255)                                                                                                                                         aGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             >  6:139177/1‑149 (MQ=255)                                                                                                                                           ccccGGTTCGCGA‑GAACCGGGGCTTTTTTAATGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            >  1:36965/1‑148 (MQ=255)                                                                                                                                            cccGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa          >  1:169530/1‑149 (MQ=255)                                                                                                                                              cGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg        <  5:359397/149‑1 (MQ=255)                                                                                                                                                gTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt      <  7:202925/149‑1 (MQ=255)                                                                                                                                                 ttCTCGC‑GAACCGGGGCTTTTTGA‑TGGTGCAGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtgt    <  6:374338/149‑1 (MQ=38)                                                                                                                                                    gcgA‑TAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  <  4:80931/149‑1 (MQ=255)                                                                                                                                                    gcgA‑GTACCGGGGCTTTTTTATTGCTGTCGGGGCGTCTTCGGAGGACAAGGCTGCGCTTGAGCTTTGTGCTGTGAAACGAGTGGATCAAGCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGtgtg   <  3:15958/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGA‑GAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551135‑1551430

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