Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	670,841	Δ1 bp	6.0%	coding (2283/2355 nt)	pbpC →	peptidoglycan glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	670,841	0	C	.	6.0%	200.6 / 6.7	67	coding (2283/2355 nt)	pbpC	peptidoglycan glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base . (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTGGTGC  >  NC_002947/670694‑670975                                                                                                                                                     |                                                                                                                                       tCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCg                                                                                                                                       >  4:236267/1‑149 (MQ=255)   tGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCggg                                                                                                                                     >  8:234957/1‑149 (MQ=255)     cgccgcCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGGGGTGGGGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCggg                                                                                                                                     <  8:146767/147‑1 (MQ=255)      gccgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAg                                                                                                                                  <  3:236267/149‑1 (MQ=255)        cgccCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGc                                                                                                                                <  4:84211/149‑1 (MQ=255)          cccGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTCTCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCg                                                                                                                              <  7:121052/149‑1 (MQ=255)             ggCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGaa                                                                                                                            >  8:123191/1‑148 (MQ=255)               ccaccaGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGCGCAGGCCGAACt                                                                                                                          >  8:307419/1‑148 (MQ=255)                   caGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAgcg                                                                                                                     >  6:224690/1‑149 (MQ=255)                      cACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGT‑‑GGCAGGCCGACCTCCGCGCGCTg                                                                                                                >  4:325584/1‑149 (MQ=255)                      cACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCAACAGGT‑‑GGCAGGCCGAACTCAGCGCGCTg                                                                                                                >  4:330858/1‑149 (MQ=255)                         tGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTg                                                                                                                >  6:66213/1‑148 (MQ=255)                          gAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGCTTGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTg                                                                                                                >  1:274590/1‑148 (MQ=255)                               gTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGaaa                                                                                                         <  7:234957/149‑1 (MQ=255)                                ttGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAg                                                                                                        >  4:2008/1‑149 (MQ=255)                                  gCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGc                                                                                                       <  6:182781/148‑1 (MQ=255)                                   cAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCgg                                                                                                     >  1:249692/1‑149 (MQ=255)                                   cAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCg                                                                                                      <  5:305660/148‑1 (MQ=255)                                        gCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCg                                                                                                      <  3:191790/143‑1 (MQ=255)                                        gCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCg                                                                                                      >  4:191790/1‑143 (MQ=255)                                         cAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGc                                                                                                                                >  6:220391/1‑116 (MQ=255)                                         cAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGc                                                                                                                                <  5:220391/116‑1 (MQ=255)                                         cAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGa                                                                                                <  2:249692/148‑1 (MQ=255)                                         cAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGAc                                                                                               >  8:68004/1‑149 (MQ=255)                                                  gCCTTGGGTGGTGGTGGGCGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGt                                                                                      >  7:49838/1‑149 (MQ=255)                                                   ccTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTg                                                                                     <  5:338298/149‑1 (MQ=255)                                                          tggtggtggACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGATGGCGTTcc                                                                              <  8:3580/149‑1 (MQ=255)                                                           ggtggtggACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                                                             <  1:69376/149‑1 (MQ=255)                                                           ggtggtggACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCa                                                                             <  6:354640/149‑1 (MQ=255)                                                            gtggtggACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAg                                                                            <  7:171455/149‑1 (MQ=255)                                                              ggtggACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGt                                                                          >  5:173166/1‑149 (MQ=255)                                                                  gACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGt                                                                      <  5:66213/149‑1 (MQ=255)                                                                  gACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGt                                                                      <  2:79112/149‑1 (MQ=255)                                                                    cGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGa                                                                    <  6:293451/149‑1 (MQ=255)                                                                       cgc‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTa                                                                 <  7:68004/149‑1 (MQ=255)                                                                         c‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAgg                                                               <  5:297842/149‑1 (MQ=255)                                                                           tggtggTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGgt                                                              >  2:284113/1‑149 (MQ=255)                                                                                gtTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTc                                                          <  7:123191/148‑1 (MQ=255)                                                                                 ttCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcga                                                        <  7:280824/149‑1 (MQ=255)                                                                                   ccTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTcgac                                                       <  8:106940/148‑1 (MQ=255)                                                                                       aaCGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACa                                                   <  5:224690/148‑1 (MQ=255)                                                                                             cAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCgg                                              >  4:52073/1‑147 (MQ=255)                                                                                             cAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCgg                                              <  3:52073/147‑1 (MQ=255)                                                                                              aGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGcc                                           >  7:326840/1‑149 (MQ=255)                                                                                              aGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGc                                            >  3:328929/1‑148 (MQ=255)                                                                                                ccGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTATGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAg                                         <  7:307419/149‑1 (MQ=255)                                                                                                   ttGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCCGGTGGCGATCCAGGTCAGTGAGTAGGTGTTCGACGACAGACGGGCCAGCAc                                      >  2:105740/1‑149 (MQ=255)                                                                                                        tGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACg                                 >  6:23658/1‑149 (MQ=255)                                                                                                           aaCCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCGACAGGTCGGGCAGGCCGAACTCAGCGCCCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCtt                              <  6:8196/149‑1 (MQ=255)                                                                                                           aaCCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCtt                              <  6:13400/149‑1 (MQ=255)                                                                                                            aCCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTg                             >  6:317190/1‑149 (MQ=255)                                                                                                             cccAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGT‑‑GGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTTGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCaa                          <  3:330858/149‑1 (MQ=255)                                                                                                             cccAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGG‑‑GGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCaa                          <  3:325584/149‑1 (MQ=255)                                                                                                              ccAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGc                            >  5:45004/1‑148 (MQ=255)                                                                                                               cAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCaa                          <  4:328929/149‑1 (MQ=255)                                                                                                                     cAGGACAGCTTGTTGCTGCGTTTGCCACAGGTTGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGt                                                                                 <  1:12794/88‑1 (MQ=255)                                                                                                                     cAGGACAGCTTGTTGCTGCGTTTGCCACAGGTTGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGt                                                                                 >  2:12794/1‑88 (MQ=255)                                                                                                                      aGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACgg                   <  6:137584/149‑1 (MQ=255)                                                                                                                        gACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGgt                 <  8:49838/149‑1 (MQ=255)                                                                                                                        gACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGgt                 >  1:258024/1‑149 (MQ=255)                                                                                                                                gttgCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCgg         >  5:47127/1‑149 (MQ=255)                                                                                                                                   gctgcGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTg      >  6:221019/1‑149 (MQ=255)                                                                                                                                    ctgcGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTgg     >  5:358868/1‑149 (MQ=255)                                                                                                                                    ctgcGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTg      <  1:284113/148‑1 (MQ=255)                                                                                                                                     tgcGTTTGCCTCAGGGCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCTGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTGGt    <  3:2008/149‑1 (MQ=255)                                                                                                                                       cGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTGGTGc  >  5:28984/1‑149 (MQ=255)                                                                                                                                                    tCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGgt                 >  4:128764/1‑121 (MQ=255)                                                                                                                                                     |                                                                                                                                       TCTGCGCCGCCCGGCCACCAGCACTGAGCCGTTGCAATTGCAGGTATCTGCCTTGGGTGGTGGTGGACGGCGC‑TGGTGGTTCCTCAACGGCCAGCCGTTGGGTGAAACCCAGGGGCAGGACAGCTTGTTGCTGCGTTTGCCACAGGTCGGGCAGGCCGAACTCAGCGCGCTGGACGAAAGCGGTCAGACGGCGCGGGTGGCGTTCCAGGTCAGTGAGTAGGTGTTCGACGACAGCCGGGCCAGCACCTACGCTTGCAAGTGACGGGTGTGACCGGCTGGTGC  >  NC_002947/670694‑670975

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