Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,033,502	G→C	8.2%	T401S (ACC→AGC) ‡	PP_3556 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,033,502	0	G	C	8.2%	177.0 / 6.4	61	T401S (ACC→AGC) ‡	PP_3556	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/22);  new base C (3/2);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGGAGAA  >  NC_002947/4033354‑4033650                                                                                                                                                     |                                                                                                                                                     cgccCGGCAAGGGCTGACGCCCTCCCCGCCTTACCAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCA‑CACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAggg                                                                                                                                                     <  8:256767/149‑2 (MQ=255)  gccCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGt                                                                                                                                                     >  7:112316/1‑149 (MQ=255)  gccCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGt                                                                                                                                                     >  3:110027/1‑149 (MQ=255)    ccGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCg                                                                                                                                                   >  1:268934/1‑149 (MQ=255)            ggCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAgg                                                                                                                                           <  4:185750/149‑1 (MQ=255)                 aCGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTc                                                                                                                                        <  1:148711/147‑1 (MQ=255)                 aCGCCCTCCCCGCCTTACAAGGTCATAGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGcccgaacgt                                                                                                                                        <  2:268934/147‑7 (MQ=255)                    cccTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTc                                                                                                                                   <  8:112316/149‑1 (MQ=255)                     ccTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCAACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTcc                                                                                                                                  >  7:81128/1‑149 (MQ=255)                           cGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCTGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGcc                                                                                                                            >  6:137491/1‑149 (MQ=255)                             ccTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCgg                                                                                                                          <  4:345164/149‑1 (MQ=255)                               ttACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGc                                                                                                                         >  7:20840/1‑148 (MQ=255)                                 aCAAGGTCATCGCTGCCAACCAGCCGAAAGCAATCAACGGCAGGTTCTAATCGATGAAGGTGGGCACAACTCTGTCACAGAAGTGGTGGTGCTTGCAATCCACTTTCAGGCCCGAGGTCTGCCCCAGCGTCTAGTACGAGGCCGGCGaa                                                                                                                      >  2:262591/1‑149 (MQ=255)                                     ggTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTg                                                                                                                  >  6:1049/1‑149 (MQ=255)                                      gTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGc                                                                                                                 <  8:138503/149‑1 (MQ=255)                                         aTCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATc                                                                                                              >  5:334004/1‑149 (MQ=255)                                              tGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGccc                                                                                                          >  4:279784/1‑148 (MQ=255)                                                    ccAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCCCCCAACGCAc                                                                                                   >  6:316329/1‑149 (MQ=255)                                                     cAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCAcc                                                                                                  >  1:258320/1‑149 (MQ=255)                                                      aGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCg                                                                                                 <  4:79935/149‑1 (MQ=255)                                                        ccGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCg                                                                                                                           <  8:137151/121‑1 (MQ=255)                                                        ccGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCgg                                                                                                >  7:27447/1‑148 (MQ=255)                                                          gAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAg                                                                                             >  6:308105/1‑149 (MQ=255)                                                              gCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTAcc                                                                                         <  8:20840/149‑1 (MQ=255)                                                              gCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGACCTCGGGcctgaac                                                                                                                                          <  2:310583/100‑6 (MQ=255)                                                              gCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGACCTCGGGcctgaac                                                                                                                                          >  1:310583/1‑95 (MQ=255)                                                                cATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCga                                                                                       >  3:111069/1‑149 (MQ=255)                                                                 aTCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCgat                                                                                      >  7:170097/1‑149 (MQ=255)                                                                 aTCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCgat                                                                                      >  1:148298/1‑149 (MQ=255)                                                                    aaCGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCgatga                                                                                    <  4:111069/148‑1 (MQ=255)                                                                         cAGGTTGTAATGGATGAAGGGGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGcaca                                                                              <  6:383259/149‑1 (MQ=255)                                                                                   tGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCaa                                                                    <  8:10493/149‑1 (MQ=255)                                                                                    ggATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAAc                                                                   >  7:197350/1‑149 (MQ=255)                                                                                    ggATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAAc                                                                   >  6:324835/1‑149 (MQ=255)                                                                                      aTGAATGTGTGCACCACGTTGTCCCAGATGTGTTGGTGCTGGCCATCCACGTTCAGGCCCTATGTCGGGCCCAGTGTCTAGTCCGAGGCCGGAGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACgg                                                                 <  1:262591/149‑1 (MQ=255)                                                                                             tGGGCACCACGGTGTCCCAGATGTGGTGGAGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAg                                                          >  3:221339/1‑149 (MQ=255)                                                                                               ggCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGc                                                         <  3:279784/148‑1 (MQ=255)                                                                                                gCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGcc                                                        >  6:259603/1‑148 (MQ=255)                                                                                                    cacGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGCGTCCGAGGCCGGCGAACCTGCATCGCCCAACGAACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGtt                                                   >  6:231786/1‑149 (MQ=255)                                                                                                    cacGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGtt                                                   >  6:244771/1‑149 (MQ=255)                                                                                                    cacGGTGTCCCAGATGTGGTGGTGATGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGtt                                                   >  6:201930/1‑149 (MQ=255)                                                                                                             cctgATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTgcacgcac                                           <  5:137491/145‑1 (MQ=255)                                                                                                                gATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTg                                       >  6:79532/1‑149 (MQ=255)                                                                                                                      gtggtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGCAGCCCAGTTGCACGCACAGTGGCACg                                  >  3:69271/1‑148 (MQ=255)                                                                                                                       tggtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACg                                  >  8:222736/1‑147 (MQ=255)                                                                                                                         gtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAg                                                                        <  1:352743/107‑1 (MQ=255)                                                                                                                         gtgCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAg                                                                        >  2:352743/1‑107 (MQ=255)                                                                                                                               gCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCgg                        <  5:324835/149‑1 (MQ=255)                                                                                                                                ccATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGc                       <  5:79532/149‑1 (MQ=255)                                                                                                                                  aTCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGc                       >  1:160653/1‑147 (MQ=255)                                                                                                                                    ccACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCAAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAg                    >  8:104660/1‑148 (MQ=255)                                                                                                                                     cACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAGGATCGCGGCCAGGa                  <  5:259603/149‑1 (MQ=255)                                                                                                                                     cACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGCTCGCGGCCCGGa                  >  4:83934/1‑149 (MQ=255)                                                                                                                                      aCGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGAt                 >  7:132674/1‑149 (MQ=255)                                                                                                                                       cGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATc                >  3:203010/1‑149 (MQ=255)                                                                                                                                         cctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCgg              >  1:165084/5‑149 (MQ=255)                                                                                                                                          ctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCg                                                                           >  3:80565/4‑87 (MQ=255)                                                                                                                                          ctgggCCCGACCTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCg                                                                           <  4:80565/84‑1 (MQ=255)                                                                                                                                              gCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAAAGGGTCGAAGCCCAGTTGAACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACgg         >  7:341568/1‑149 (MQ=255)                                                                                                                                                ccGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTg       <  5:316329/149‑1 (MQ=255)                                                                                                                                                ccGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTg       <  5:1049/149‑1 (MQ=255)                                                                                                                                                ccGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTg       <  2:148298/149‑1 (MQ=255)                                                                                                                                                    ggTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTTCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGgaga   <  5:201930/149‑1 (MQ=255)                                                                                                                                                     gTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGGAGaa  <  5:363952/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGATGAAGGTGGGCACCACGGTGTCCCAGATGTGGTGGTGCTGGCCATCCACGTTCAGGCCCGAGGTCGGGCCCAGGGTCGAGTCCGAGGCCGGCGAACCTGCATCGCCCAACGCACCGGCAGTACCGATGATGCACACCGTAGCCAACGGGTCGAAGCCCAGTTGCACGCACAGTGGCACGAAGATCGCGGCCAGGATCGGCACGGTGGAGAA  >  NC_002947/4033354‑4033650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑