Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,695,005	T→C	5.9%	L324S (TTG→TCG)	wspC →	biofilm formation methyltransferase WspC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,695,005	0	T	C	5.9%	200.6 / 2.2	67	L324S (TTG→TCG)	wspC	biofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GCCCAAACTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCG  >  NC_002947/1694857‑1695147                                                                                                                                                      |                                                                                                                                                 gCCCAAACTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGtt                                                                                                                                                  <  5:366861/149‑1 (MQ=255)    cAAACTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑t                                                                                                                                              >  8:333559/1‑149 (MQ=255)     aaaCTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑Tg                                                                                                                                             >  6:363093/1‑149 (MQ=255)     aaaCTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑Tg                                                                                                                                             <  3:63808/149‑1 (MQ=255)       aCTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGc                                                                                                                                           >  2:324064/1‑149 (MQ=255)           ccAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGTCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑Cg                                                                                                                                        <  2:254203/147‑1 (MQ=255)           ccAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑Cgg                                                                                                                                       >  1:79501/1‑148 (MQ=255)           ccAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGc                                                                                                                                      <  6:112911/149‑1 (MQ=255)           ccAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCAACTCGC‑TGGc‑ccc                                                                                                                                       <  2:178924/148‑3 (MQ=255)            cAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCa                                                                                                                                     <  1:66345/149‑1 (MQ=255)              aGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCAtt                                                                                                                                   <  7:333559/149‑1 (MQ=255)                    cgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCg                                                                                                                              <  4:85563/148‑1 (MQ=255)                     gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGc                                                                                                                            >  6:42110/1‑149 (MQ=255)                         acgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTgg                                                                                                                         <  4:276447/148‑1 (MQ=255)                         acgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGc                                                                                                                        >  3:292017/1‑149 (MQ=255)                          cgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGc                                                                                                                        >  1:57124/1‑148 (MQ=255)                           gacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGcc                                                                                                                       >  1:376790/1‑148 (MQ=255)                               ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCGGTTGCCTCGTTCTCGTCGTGTGCTGCCGGTGGCTGCCAGGCCAGCGCGAGACCGCGAGCA‑CAGCCACGTGGGGGACAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGc                                                                                                                  <  2:250129/149‑1 (MQ=255)                                  cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGGCCAGCCGGGCAATGCCg                                                                                                                <  8:136359/148‑1 (MQ=255)                                  cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                <  3:198722/148‑1 (MQ=255)                                  cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCCGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CCGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCCCCTGGCCAATGCCgg                                                                                                               >  4:21079/1‑149 (MQ=255)                                      cACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                            >  7:287867/1‑148 (MQ=255)                                         cGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGc                                                                                                        <  8:131187/149‑1 (MQ=255)                                           gTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGc                                                                                                        >  7:185898/1‑147 (MQ=255)                                           gTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                      <  6:341404/149‑1 (MQ=255)                                               aTCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAg                                                                                                  >  5:78927/1‑149 (MQ=255)                                                       cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCAGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAATAGCAACCCCCCGGGGGGCCGCGAGTTTG‑GGGC‑CCGCATTTCCCGGGTGGCCAATTCCCGTGGCCGCGCGCCGGCCCg                                                                                             >  1:314946/1‑147 (MQ=255)                                                        ccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                         >  5:321971/1‑149 (MQ=255)                                                            tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                     <  2:57124/149‑1 (MQ=255)                                                                 cGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCaa                                                                                >  6:127345/1‑149 (MQ=255)                                                                     gCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGc                                                                             >  4:204461/1‑148 (MQ=255)                                                                       cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                          >  8:235242/1‑149 (MQ=255)                                                                           ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                      <  5:110853/149‑1 (MQ=255)                                                                           ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                      <  2:376790/149‑1 (MQ=255)                                                                              gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                    <  5:141198/148‑1 (MQ=255)                                                                              gcggcgGGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                    <  3:288498/148‑1 (MQ=255)                                                                               cggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                  <  6:51878/149‑1 (MQ=255)                                                                                     gtgCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAccc                                                            <  6:321971/149‑1 (MQ=255)                                                                                           ccGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTc                                                      >  5:21321/1‑149 (MQ=255)                                                                                           ccGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTc                                                      >  6:2227/1‑149 (MQ=255)                                                                                            cGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCg                                                     <  4:380992/149‑1 (MQ=255)                                                                                               tGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGca                                                  <  4:292017/149‑1 (MQ=255)                                                                                                 gcggcCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGcac                                                 >  5:201354/1‑148 (MQ=255)                                                                                                 gcggcCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGcac                                                 <  6:201354/148‑1 (MQ=255)                                                                                                   ggccaggccaGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAgg                                              >  6:57519/1‑149 (MQ=255)                                                                                                    gccaggccagCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                             <  5:42110/149‑1 (MQ=255)                                                                                                      caggccagCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtgt                                           <  5:127345/149‑1 (MQ=255)                                                                                                       aggccagCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGCGCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                            >  1:295508/1‑148 (MQ=255)                                                                                                         gccagCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtac                                         >  4:45076/1‑148 (MQ=255)                                                                                                         gccagCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtac                                         >  7:216074/1‑148 (MQ=255)                                                                                                             gcgcGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActg                                    >  6:192088/1‑149 (MQ=255)                                                                                                             gcgcGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActg                                    <  7:214397/149‑1 (MQ=255)                                                                                                              cgcgAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                   <  5:192088/149‑1 (MQ=255)                                                                                                                 gAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctg                                >  6:298503/1‑149 (MQ=255)                                                                                                                 gAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctg                                <  8:194870/149‑1 (MQ=255)                                                                                                                     ccGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGCGTGCC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATCCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGt                             >  1:253861/1‑149 (MQ=255)                                                                                                                            ca‑caGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCATCGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTc                       <  7:266432/147‑1 (MQ=255)                                                                                                                            ca‑caGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCACCGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTc                       >  8:266432/1‑147 (MQ=255)                                                                                                                             a‑caGCCTCGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATGGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                    <  8:266283/149‑1 (MQ=255)                                                                                                                             a‑caGCCACGAGGGGGCCAGCGAGTTGG‑TGGC‑CGGCATTGCGCGGCTGGCCAATGCCGGTGCCAGGGAGCAGGCCGGTAGCGAATGCCAACGCTACCTGAGCGAGTACCCGCCGTCGGGACAGGTGTACTACTGGCTGGGTTTGCTCatt                    <  3:21079/149‑3 (MQ=255)                                                                                                                                    caatgccggCCAGCAACTCGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAcc              >  6:380959/8‑149 (MQ=255)                                                                                                                                            gCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa      <  5:363093/149‑1 (MQ=255)                                                                                                                                            gCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa      >  7:42465/1‑149 (MQ=255)                                                                                                                                             ccAGCAACTCGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                             >  2:164901/1‑109 (MQ=255)                                                                                                                                             ccAGCAACTCGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                             <  1:164901/109‑1 (MQ=255)                                                                                                                                               aGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc   >  8:61115/1‑149 (MQ=255)                                                                                                                                                gCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgcg  <  5:2227/149‑1 (MQ=255)                                                                                                                                                gCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgcg  >  2:123264/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCCCAAACTGCCAAGCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCA‑CAGCCACGAGGGGGCCAGCGAGTTGC‑TGGC‑CGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCG  >  NC_002947/1694857‑1695147

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