Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,551,282	G→T	8.8%	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.8%	177.9 / 4.9	63	intergenic (+41/+85)	groL/pykA	chaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base T (2/4);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

GTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551135‑1551430                                                                                                                                                    |                                                                                                                                                      gTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                      <  7:101563/149‑1 (MQ=255)          tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACC‑gg                                                                                                                                            <  6:72513/149‑1 (MQ=255)           caccacCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑ggg                                                                                                                                           >  8:229145/1‑149 (MQ=255)            accaccGAAGCCATGGTCGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACC‑gggg                                                                                                                                          <  3:82175/149‑1 (MQ=255)              caccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgcg                                                                                                                                                     >  6:129816/1‑137 (MQ=255)              caccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgcg                                                                                                                                                     <  5:129816/137‑1 (MQ=255)              caccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCt                                                                                                                                        >  8:305155/1‑149 (MQ=255)                     gCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑gg                                                                                                                                            <  6:100703/138‑1 (MQ=255)                     gCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑gg                                                                                                                                            >  5:100703/1‑138 (MQ=255)                      ccATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACC‑GGGGCTTTTTGAtg                                                                                                                                <  1:163061/149‑2 (MQ=255)                       cATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTg                                                                                                                               <  7:305155/149‑1 (MQ=255)                         tGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCt                                                                                                                             <  3:40864/149‑1 (MQ=255)                               ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCgg                                                                                                                        <  7:303464/148‑1 (MQ=255)                               ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATCGGTGGCTTGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCgg                                                                                                                        <  1:295464/148‑1 (MQ=255)                                     tGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCctct                                                                                                                 >  7:225748/1‑149 (MQ=255)                                      gCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCtt                                                                                                                >  2:106754/1‑149 (MQ=255)                                       ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTc                                                                                                               <  6:163897/149‑1 (MQ=255)                                       ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTc                                                                                                               <  7:267353/149‑1 (MQ=255)                                          gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCg                                                                                                              >  4:340048/1‑147 (MQ=255)                                                aaGCCAGGTGCTGGCATGCCTGATATGGGTGGGATGGGTTGCTTGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTGTCGCGAACC‑GGGGCtttt                                                                                                                                     <  6:310631/118‑1 (MQ=255)                                                aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  7:125935/1‑147 (MQ=255)                                                             gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAg                                                                                         <  8:225748/149‑1 (MQ=255)                                                                  ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGACGGGGCCTCTTCGCACCACAAGGCTGCGCTTGAGCCTtg                                                                                    >  6:256721/1‑149 (MQ=255)                                                                       tatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                               <  5:110830/149‑1 (MQ=255)                                                                         tgggtggcatgggtggcaTGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtgtg                                                                             <  6:199333/149‑1 (MQ=255)                                                                          gggtggcatgggtggcatGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtgtg                                                                             >  7:99481/1‑148 (MQ=255)                                                                             tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         >  5:58534/1‑149 (MQ=255)                                                                                catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAg                                                                      >  5:251694/1‑149 (MQ=255)                                                                                  tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTg                                                                    >  7:302605/1‑149 (MQ=255)                                                                                  tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGCTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGCGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTg                                                                    <  8:283123/149‑1 (MQ=255)                                                                                     gtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGAt                                                                 <  3:340048/149‑1 (MQ=255)                                                                                      tggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGAt                                                                 >  1:64228/1‑148 (MQ=255)                                                                                        gcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCaa                                                              >  2:340213/1‑149 (MQ=255)                                                                                          atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAcc                                                            >  1:339829/1‑149 (MQ=255)                                                                                          atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAcc                                                            <  8:302605/149‑1 (MQ=255)                                                                                          atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCACGGTTAGAGAGAACC‑GGGGCATTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAATGCTGCGCTTGAGCCTAGTGCTGTGAAACGAGTGGATCAAcc                                                            >  6:308648/1‑149 (MQ=255)                                                                                           tgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAcc                                                            >  3:92365/1‑148 (MQ=255)                                                                                                    tgatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGa                                                  <  2:339829/149‑1 (MQ=255)                                                                                                     gatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  5:62573/129‑1 (MQ=255)                                                                                                     gatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  6:62573/1‑129 (MQ=255)                                                                                                      atgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCGGTGAAACGACGGGATCAACCGCCCCGCGGAt                                                 >  8:101532/1‑148 (MQ=255)                                                                                                      atgTAAGCAAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGAGGTGAAACGAGTGGATCAACAGACCAGGGGAt                                                 >  8:106697/1‑148 (MQ=255)                                                                                                          aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCg                                            >  3:191943/1‑149 (MQ=255)                                                                                                            gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                          <  7:229145/149‑1 (MQ=255)                                                                                                              cagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGt                                        <  6:251694/149‑1 (MQ=255)                                                                                                               agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCTGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGt                                        >  1:261169/1‑149 (MQ=255)                                                                                                                     aCCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGa                                  <  2:207823/148‑1 (MQ=255)                                                                                                                      ccccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGCGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCACCCGACCAGCGGATCACCGACGTGCAGGATc                                <  8:125935/149‑1 (MQ=255)                                                                                                                         cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  2:64228/149‑1 (MQ=255)                                                                                                                         cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGACGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  5:256721/149‑1 (MQ=255)                                                                                                                          ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCctct                                                                                                                 >  4:256805/1‑64 (MQ=255)                                                                                                                          ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCctct                                                                                                                 <  3:256805/64‑1 (MQ=255)                                                                                                                          ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  8:26657/108‑1 (MQ=255)                                                                                                                          ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  7:26657/1‑108 (MQ=255)                                                                                                                          ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCa                            >  5:106802/1‑149 (MQ=255)                                                                                                                           cTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                           >  8:149042/1‑149 (MQ=255)                                                                                                                               aCCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCa                       <  6:106802/149‑1 (MQ=255)                                                                                                                                ccATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAt                      <  7:149042/149‑1 (MQ=255)                                                                                                                                      aaaaGCCCCGGTTCTCGCGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc                >  2:165577/1‑149 (MQ=255)                                                                                                                                         aGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGTAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             <  1:106754/149‑1 (MQ=255)                                                                                                                                          gCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            >  7:236666/1‑149 (MQ=255)                                                                                                                                          gCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            >  3:87584/1‑149 (MQ=255)                                                                                                                                            cccGGTTCTCGCGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCATCATGCCATTGGTGCCACCGa          >  7:309304/1‑149 (MQ=255)                                                                                                                                              cGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGCGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg        >  5:306292/1‑149 (MQ=255)                                                                                                                                                gTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt      <  4:87584/149‑1 (MQ=255)                                                                                                                                                  tCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtgt    >  3:333940/1‑149 (MQ=255)                                                                                                                                                    gcgAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  <  8:99481/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GTCTGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACC‑GGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  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: ‑