Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,528,705	T→C	5.9%	V349A (GTA→GCA)	ftsA →	cell division protein FtsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,528,705	0	T	C	5.9%	199.7 / 2.0	68	V349A (GTA→GCA)	ftsA	cell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GACGAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAA  >  NC_002947/1528563‑1528852                                                                                                                                               |                                                                                                                                                     gacgaGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCt                                                                                                                                                >  3:139588/1‑149 (MQ=255)            cACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGTGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGc                                                                                                                                        <  4:124575/146‑1 (MQ=255)               ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGACCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑Gca                                                                                                                                   >  7:283591/1‑147 (MQ=255)               ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑Gca                                                                                                                                   >  5:63800/1‑147 (MQ=255)               ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑Gca                                                                                                                                   >  4:166935/1‑147 (MQ=255)                      cAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACC‑ggcat                                                                                                                                 >  1:202072/1‑138 (MQ=255)                      cAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACC‑ggcat                                                                                                                                 <  2:202072/142‑5 (MQ=255)                       aGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCgg                                                                                                                         >  7:143388/1‑149 (MQ=255)                                tGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGc                                                                                                                >  6:254901/1‑149 (MQ=255)                                tGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGc                                                                                                                >  2:260215/1‑149 (MQ=255)                                   gccgcAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGAc                                                                                                             >  3:182994/1‑149 (MQ=255)                                      gcAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACgtg                                                                                                          <  6:63800/149‑1 (MQ=255)                                        aGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACgtggt                                                                                                        <  5:254901/149‑1 (MQ=255)                                               aCGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGAAGTAGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTAGGGGGCTGAGCGACGTGGTGCGCAAc                                                                                                 >  8:48648/1‑149 (MQ=255)                                                cGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCaa                                                                                                  <  8:309026/147‑1 (MQ=255)                                                 gAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAAccc                                                                                               >  7:45190/1‑149 (MQ=255)                                                    gACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            <  8:45190/149‑1 (MQ=255)                                                       cTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCaa                                                                                                  <  2:229454/140‑1 (MQ=255)                                                       cTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCaa                                                                                                  >  1:229454/1‑140 (MQ=255)                                                       cTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            <  5:48026/146‑1 (MQ=255)                                                       cTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            >  6:48026/1‑146 (MQ=255)                                                       cTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTa                                                                                         >  3:286290/1‑149 (MQ=255)                                                       cTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTa                                                                                         >  7:99224/1‑149 (MQ=255)                                                          gTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTc                                                                                      >  3:45916/1‑149 (MQ=255)                                                             cccGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCAc                                                                                   >  6:102079/1‑149 (MQ=255)                                                                gccgccATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACgtggt                                                                                                        <  5:107343/125‑1 (MQ=255)                                                                gccgccATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACgtggt                                                                                                        >  6:107343/1‑125 (MQ=255)                                                                gccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            >  4:127980/1‑137 (MQ=255)                                                                gccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            <  3:127980/137‑1 (MQ=255)                                                                 ccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAAcc                                                                                                <  4:282939/132‑1 (MQ=255)                                                                 ccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAAcc                                                                                                >  3:282939/1‑132 (MQ=255)                                                                  cgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGgtg                                                                              >  3:69537/1‑149 (MQ=255)                                                                   ggCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGgtgt                                                                             >  7:217586/1‑149 (MQ=255)                                                                   ggCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGgtgt                                                                             <  8:99224/149‑1 (MQ=255)                                                                   ggCATCGGCCTCCCCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGgtgt                                                                             <  8:283591/149‑1 (MQ=255)                                                                     cATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTgg                                                                           <  6:13674/149‑1 (MQ=255)                                                                          tCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTg                                                                      <  4:139588/149‑1 (MQ=255)                                                                           ccTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGc                                                                         <  8:217586/145‑1 (MQ=255)                                                                             tCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCt                                                                    >  1:236202/1‑148 (MQ=255)                                                                                  ggcggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            >  8:95985/1‑119 (MQ=255)                                                                                  ggcggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                            <  7:95985/119‑1 (MQ=255)                                                                                    cggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTTCGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACg                                                            <  6:137868/149‑1 (MQ=255)                                                                                     ggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACgg                                                           <  4:45916/149‑1 (MQ=255)                                                                                              aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                  <  8:16946/149‑1 (MQ=255)                                                                                              aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                  <  4:69537/149‑1 (MQ=255)                                                                                              aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                  <  1:260215/149‑1 (MQ=255)                                                                                              aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTCGCTCACCTACgg                                                           <  5:316709/140‑1 (MQ=255)                                                                                              aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTCGCTCACCTACgg                                                           >  6:316709/1‑140 (MQ=255)                                                                                                gATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAg                                                 >  5:95087/1‑148 (MQ=255)                                                                                                        gcgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAg                                        <  8:143388/149‑1 (MQ=255)                                                                                                           ccGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                     <  6:95087/149‑1 (MQ=255)                                                                                                            cGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAg                                        >  2:44947/1‑145 (MQ=255)                                                                                                            cGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAg                                        <  1:44947/145‑1 (MQ=255)                                                                                                             gTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGGGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACtt                                   <  7:48648/149‑1 (MQ=255)                                                                                                                gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGc                                >  6:131082/1‑149 (MQ=255)                                                                                                                 aaCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                               <  3:166935/149‑1 (MQ=255)                                                                                                                 aaCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                               <  5:102079/149‑1 (MQ=255)                                                                                                                      gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTgg                                                                           <  4:67126/100‑1 (MQ=255)                                                                                                                      gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTgg                                                                           >  3:67126/1‑100 (MQ=255)                                                                                                                       ccGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGAcc                         <  4:182994/149‑1 (MQ=255)                                                                                                                            gAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaa                    <  2:236202/149‑1 (MQ=255)                                                                                                                             aaaTCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                     <  4:282666/131‑1 (MQ=255)                                                                                                                                      acaTGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaacag           >  7:24772/1‑148 (MQ=255)                                                                                                                                        tgtgcGGCACGCCCAGGCGTACCGGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaaca                  >  5:34647/5‑140 (MQ=255)                                                                                                                                        tgtgcGGCACGCCCAGGCGTACCGGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaaca                  <  6:34647/136‑1 (MQ=255)                                                                                                                                          gCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctg                                                       >  2:201276/1‑100 (MQ=255)                                                                                                                                          gCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctg                                                       <  1:201276/100‑1 (MQ=255)                                                                                                                                              gTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCaa  >  2:80457/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     GACGAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCC‑GCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAA  >  NC_002947/1528563‑1528852

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