Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,018,622	G→T	6.3%	G87V (GGT→GTT)	bcsZ →	cellulose synthase endo‑1,4‑D‑glucanase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,018,622	0	G	T	6.3%	203.5 / 3.8	63	G87V (GGT→GTT)	bcsZ	cellulose synthase endo‑1,4‑D‑glucanase subunit
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base T (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

TGGATGGGCGGGTCATCGACCCGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGTTGTGGGACCAGCCGGCC  >  NC_002947/3018475‑3018770                                                                                                                                                    |                                                                                                                                                       tgGATGGGCGGGTCATCGACCCGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAggg                                                                                                                                                        <  3:29873/148‑1 (MQ=255)           ggTCATCGACCCGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGccc                                                                                                                                             <  8:198650/149‑1 (MQ=255)                  gACCCGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATc                                                                                                                                      <  6:364209/148‑1 (MQ=255)                  gACCCGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATc                                                                                                                                      <  6:299966/148‑1 (MQ=255)                    cccGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATctgc                                                                                                                                   <  8:169425/149‑1 (MQ=255)                        aGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATctgcctgc                                                                                                                               <  1:159017/149‑1 (MQ=255)                          tGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTgcctgcct                                                                                                                             <  7:244411/149‑1 (MQ=255)                                gcgcCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCtgt                                                                                                                       >  1:33141/1‑149 (MQ=255)                                  gcCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTgg                                                                                                                     <  5:278531/149‑1 (MQ=255)                                   ccTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTggg                                                                                                                    >  7:11463/1‑149 (MQ=255)                                        tCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGt                                                                                                               >  4:278413/1‑149 (MQ=255)                                            caccTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACg                                                                                                           >  6:102566/1‑149 (MQ=255)                                             accTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCt                                                                                                                                                  <  5:327889/110‑1 (MQ=255)                                             accTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCt                                                                                                                                                  >  6:327889/1‑110 (MQ=255)                                              ccTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCtgt                                                                                                                       <  4:230811/135‑1 (MQ=255)                                              ccTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCtgt                                                                                                                       >  3:230811/1‑135 (MQ=255)                                               cTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACggg                                                                                                         >  6:122883/1‑148 (MQ=255)                                                      gggCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGgcagcag                                                                                                   >  1:43435/1‑147 (MQ=255)                                                              ttACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAg                                                                                          <  8:11463/148‑1 (MQ=255)                                                                 cgcgCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGc                                                                                      >  1:95787/1‑149 (MQ=255)                                                                  gcgcTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCt                                                                                     >  3:272557/1‑149 (MQ=255)                                                                   cgcTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTg                                                                                    <  5:102566/149‑1 (MQ=255)                                                                    gcTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGTTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTgg                                                                                   >  7:126111/1‑149 (MQ=255)                                                                    gcTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGTTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTg                                                                                    >  7:130879/1‑148 (MQ=255)                                                                        ttcttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGAc                                                                                 <  2:95787/147‑1 (MQ=255)                                                                        ttcttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACg                                                                                <  4:360/148‑1 (MQ=255)                                                                          cttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGt                                                                                               >  2:344819/1‑131 (MQ=255)                                                                          cttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGt                                                                                               <  1:344819/131‑1 (MQ=255)                                                                           ttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATctgcc                                                                                                                                  <  1:184197/95‑1 (MQ=255)                                                                           ttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATctgcc                                                                                                                                  >  2:184197/1‑95 (MQ=255)                                                                           ttcGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGAGGCAGGTGCTGGACGCcaa                                                                            >  6:203190/1‑149 (MQ=255)                                                                               cccTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCcaa                                                                            <  6:334773/145‑1 (MQ=255)                                                                               cccTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCcaa                                                                            >  5:334773/1‑145 (MQ=255)                                                                               cccTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCaacaac                                                                        >  2:89726/1‑149 (MQ=255)                                                                                ccTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGgcagcagca                                                                                                 <  4:136406/123‑1 (MQ=255)                                                                                ccTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGgcagcagca                                                                                                 >  3:136406/1‑123 (MQ=255)                                                                                ccTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACg                                                                       >  4:149737/1‑149 (MQ=255)                                                                                    gtaggtaACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCa                                                                    <  2:33141/148‑1 (MQ=255)                                                                                     taggtaACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGTTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGc                                                                  <  8:126111/149‑1 (MQ=255)                                                                                     taggtaACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGTTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGc                                                                  <  8:130879/149‑1 (MQ=255)                                                                                          aaCGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACATCGCCGAGGGTGACCTGGCCC‑GGCATATGCCTGCCAGGCTGTGGGGCCGGGACGGGCAGCAGCAGTGGAAGGTGCTGGACGCCCACAAAGCAAGCGATGc                                                             >  8:49804/1‑149 (MQ=255)                                                                                           aCGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGc                                                             <  5:122883/148‑1 (MQ=255)                                                                                             gACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCg                                                                 <  4:310583/142‑1 (MQ=255)                                                                                             gACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCg                                                                 >  3:310583/1‑142 (MQ=255)                                                                                                 ggCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCGGCAGCAGTGGCAGGTGCTGGACGCCAACAACGc                                                                      <  5:302436/133‑1 (MQ=255)                                                                                                 ggCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCGGCAGCAGTGGCAGGTGCTGGACGCCAACAACGc                                                                      >  6:302436/1‑133 (MQ=255)                                                                                                 ggCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCtg                                                      <  3:278413/149‑1 (MQ=255)                                                                                                  gCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCtgt                                                     <  5:203190/149‑1 (MQ=255)                                                                                                  gCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCtgt                                                     <  3:149737/149‑1 (MQ=255)                                                                                                            gcgcAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTa                                           >  2:14848/1‑149 (MQ=255)                                                                                                              gcAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACa                                         <  1:89726/149‑1 (MQ=255)                                                                                                              gcAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACa                                         <  2:43435/149‑1 (MQ=255)                                                                                                               cAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCCGGGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTa                                           >  1:261898/1‑147 (MQ=255)                                                                                                                    gctgcGTTGGACCAGCAACAACATCGCCAAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCAGTGGGGACGTGACGTGCAGCAGCAGTCGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCctgc                                   >  1:343330/1‑149 (MQ=255)                                                                                                                       gcGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCctgctg                                 >  2:246666/1‑148 (MQ=255)                                                                                                                                  gCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACa                                         >  8:21551/1‑129 (MQ=255)                                                                                                                                  gCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACa                                         <  7:21551/129‑1 (MQ=255)                                                                                                                                  gCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCgg                    >  1:44156/1‑149 (MQ=255)                                                                                                                                   caacaaCCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGTCGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGt                   >  1:274890/1‑149 (MQ=255)                                                                                                                                    aacaacCTCGCCGAGGGTGACCTGGCCC‑GGCATCTTCCTGCCTGGCTGTGGGCCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACTCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTT‑CCAGCAggctgtag                  <  7:49804/149‑7 (MQ=255)                                                                                                                                         ccTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGTTGTggg             <  2:44156/149‑1 (MQ=255)                                                                                                                                              ccGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCctgct                                  >  6:356681/1‑124 (MQ=255)                                                                                                                                              ccGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCctgct                                  <  5:356681/124‑1 (MQ=255)                                                                                                                                               cGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGTTGTGGGACCAGc       >  5:53887/1‑149 (MQ=255)                                                                                                                                                  gggTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTGCGAAGCAGGGCGGTTGTGGGACCAGCCg     <  1:14848/149‑1 (MQ=255)                                                                                                                                                   ggTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGTTGTGGGACCAGCCGGc   >  3:204612/1‑149 (MQ=255)                                                                                                                                                    gTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGTTGTGGGACCAGCCGGcc  <  2:274890/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TGGATGGGCGGGTCATCGACCCGAGTGACGAGCGCCTGATCACCACCTCGGAAGGGCAGAGTTACGCGCTGTTCTTCGCCCTGGTAGGTAACGACCGGCAAACCTTCGCGCAGTTGCTGCGTTGGACCAGCAACAACCTCGCCGAGGGTGACCTGGCCC‑GGCATCTGCCTGCCTGGCTGTGGGGCCGTGACGGGCAGCAGCAGTGGCAGGTGCTGGACGCCAACAACGCCAGCGATGCCGATCTGTGGATCGCCTACAGCCTGCTG‑GAAGCAGGGCGGTTGTGGGACCAGCCGGCC  >  NC_002947/3018475‑3018770

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