Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,616,319	A→G	5.7%	D897G (GAC→GGC)	PP_4085 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,616,319	0	A	G	5.7%	199.6 / 2.1	70	D897G (GAC→GGC)	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

TTACCAACGAGAATGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTG  >  NC_002947/4616172‑4616467                                                                                                                                                    |                                                                                                                                                     ttACCAACGAGAATGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGAc                                                                                                                                                     <  6:279538/149‑1 (MQ=33)      aaCGAGAATGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTTCAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGag                                                                                                                                                  >  1:99483/1‑147 (MQ=17)        cGAGAATGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGagcagc                                                                                                                                              >  5:267354/1‑149 (MQ=32)         gagaATGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAgcc                                                                                                                                             <  6:49412/149‑1 (MQ=32)            aaTGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAgccgcc                                                                                                                                          >  7:211464/1‑149 (MQ=32)            aaTGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAgccgcc                                                                                                                                          >  4:45025/1‑149 (MQ=32)            aaTGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAgccgcc                                                                                                                                          >  2:80995/1‑149 (MQ=32)            aaTGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAgccgcc                                                                                                                                          >  4:286613/1‑149 (MQ=32)                      cATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGa                                                                                                                                 <  4:104983/148‑1 (MQ=32)                           cGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGc                                                                                                                           >  5:382422/1‑149 (MQ=37)                                    tGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGa                                                                                                                   >  7:282745/1‑148 (MQ=37)                                      ggACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGacac                                                                                                                <  1:349262/149‑1 (MQ=37)                                       gACGAACTGGACCGGCTGGTGGCACAAAAAGAGCTGGATGGCAGTGGCCGCCTCTATGCCGACAAGGTGGTGGATGAGGTCAGTCGTGTTGCGCATGTGCCCTCGCCGGACGAGCAGCCGCGCCTGTCGGACAACGCCCCGGCGACACg                                                                                                               >  8:241224/1‑149 (MQ=11)                                          gAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGgac                                                                                                            <  2:133357/149‑1 (MQ=37)                                            aCTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGgacga                                                                                                          >  3:245607/1‑149 (MQ=38)                                                gACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTg                                                                                                       <  6:267354/148‑1 (MQ=39)                                                gACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTg                                                                                                       <  2:60605/148‑1 (MQ=39)                                                     gCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATcc                                                                                                 >  5:152164/1‑149 (MQ=38)                                                              tcaAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACt                                                                                        >  3:230400/2‑149 (MQ=37)                                                                aaaaaGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTc                                                                                      >  2:237562/1‑149 (MQ=37)                                                                        cTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTACCATCCGCCATGACTTCGAGCGGGa                                                                              >  3:373557/1‑149 (MQ=37)                                                                            aTGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGcc                                                                          <  1:80995/149‑1 (MQ=37)                                                                              ggCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGAAACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCAt                                                                        >  3:261596/1‑149 (MQ=37)                                                                               gCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGc                                                                                  >  8:90972/1‑138 (MQ=37)                                                                               gCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGc                                                                                  <  7:90972/138‑1 (MQ=37)                                                                                  gTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATcggc                                                                    >  4:284219/1‑149 (MQ=37)                                                                                   tGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATcggcc                                                                   >  1:156954/1‑149 (MQ=37)                                                                                    ctccgccTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATcggccg                                                                  <  7:325433/147‑1 (MQ=37)                                                                                      ccgccTCTATGCCTACAACGTGCGGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATcggccggc                                                                <  3:286613/149‑1 (MQ=37)                                                                                       cgccTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATcggccg                                                                  >  8:325433/1‑146 (MQ=37)                                                                                        gccTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCtt                                                              <  6:382422/149‑1 (MQ=37)                                                                                          ctctATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGt                                                            <  4:245607/149‑1 (MQ=37)                                                                                               tGCCTACAACTTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCa                                                          <  8:211464/146‑1 (MQ=37)                                                                                               tGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCa                                                          <  3:330041/146‑1 (MQ=37)                                                                                               tGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCa                                                       <  6:152164/149‑1 (MQ=37)                                                                                                gCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGcc                                                                          >  7:299787/1‑129 (MQ=38)                                                                                                gCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGcc                                                                          <  8:299787/129‑1 (MQ=38)                                                                                                 ccTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCaa                                                      <  1:237562/148‑1 (MQ=37)                                                                                                  cttcAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAgc                                                    <  2:99483/146‑1 (MQ=37)                                                                                                  cTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCTGGCTTGt                                                            <  5:130819/141‑1 (MQ=37)                                                                                                  cTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCGGGCTTGt                                                            >  6:130819/1‑141 (MQ=37)                                                                                                     cAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCa                                                 >  4:147883/1‑149 (MQ=37)                                                                                                        cGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGGCGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGTGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCAcc                                               >  7:237667/1‑148 (MQ=255)                                                                                                        cGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGGCGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGTGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCAcc                                               >  7:231600/1‑148 (MQ=255)                                                                                                          tGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGGCGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGTGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGaa                                            <  8:237667/149‑1 (MQ=255)                                                                                                          tGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGGCGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGTGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGaa                                            <  8:231600/149‑1 (MQ=255)                                                                                                                aTGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGc                                      <  3:45025/149‑1 (MQ=37)                                                                                                                aTGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGc                                      <  5:8466/149‑1 (MQ=37)                                                                                                                     gTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCAACCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGcaccac                                 >  5:24009/1‑149 (MQ=37)                                                                                                                              cTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCg                        >  6:95772/1‑149 (MQ=37)                                                                                                                                 aCCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGca                                     >  8:364818/1‑133 (MQ=37)                                                                                                                                 aCCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGc                       <  8:282745/147‑1 (MQ=37)                                                                                                                                 aCCCATGTGCCCTCGCCTGACGAGCAGCCGCCACTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGca                                     <  7:364818/133‑1 (MQ=21)                                                                                                                                  cccATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCGTCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTac                    <  4:12526/149‑1 (MQ=37)                                                                                                                                   ccATGTGCCCTCGCCTGACGAGCGGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGc                       <  1:110778/145‑1 (MQ=37)                                                                                                                                   ccATGTGCCCTCGCCTGACGAGCGGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGc                       >  2:110778/1‑145 (MQ=37)                                                                                                                                     aTGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACtt                                                                                       <  7:80544/79‑1 (MQ=255)                                                                                                                                     aTGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACtt                                                                                       >  8:80544/1‑79 (MQ=255)                                                                                                                                     aTGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGc                                                        <  4:41529/110‑1 (MQ=35)                                                                                                                                     aTGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGc                                                        >  3:41529/1‑110 (MQ=35)                                                                                                                                      tgtgCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTacgacg                >  1:340667/1‑149 (MQ=37)                                                                                                                                        tgCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGAcgcg              >  6:43154/1‑149 (MQ=37)                                                                                                                                          cccTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGCCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGc            >  4:81952/1‑149 (MQ=37)                                                                                                                                           ccTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGcc           >  8:319611/1‑149 (MQ=37)                                                                                                                                            cTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCg          <  5:95772/149‑1 (MQ=37)                                                                                                                                               gCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGAc        >  2:90471/1‑148 (MQ=37)                                                                                                                                               gCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCAACCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGa         <  6:24009/147‑1 (MQ=37)                                                                                                                                                 cTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAAt     >  6:397650/1‑149 (MQ=37)                                                                                                                                                  tGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGCCGCGGCCGACAATc    >  8:323187/1‑149 (MQ=37)                                                                                                                                                   gACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTAc                          <  3:98437/126‑1 (MQ=38)                                                                                                                                                    aCGAGCAGCCGCGCCGGGCCGACAACGCCCCGCCGTCACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATctg  <  4:261596/149‑1 (MQ=255)                                                                                                                                                    aCGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATctg  <  3:81952/149‑1 (MQ=37)                                                                                                                                                    |                                                                                                                                                     TTACCAACGAGAATGACGAGTCATACCGCTTCGGCTGGGACGAACTGGACCGGCTGGTGGCACAAAAAGACCTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTG  >  NC_002947/4616172‑4616467

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