Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,616,385	A→T	5.4%	E919V (GAG→GTG)	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,385	0	A	T	5.4%	248.9 / 3.1	84	E919V (GAG→GTG)	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (2/3);  total (42/42)
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

CTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTACACCTACGACGC  >  NC_002947/4616243‑4616529                                                                                                                                               |                                                                                                                                                 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                                                                                                                                        <  8:299787/129‑1 (MQ=38)                         gCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGcc                                                                                                                                        >  7:299787/1‑129 (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:231600/149‑1 (MQ=255)                                   tGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGGCGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGTGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGaa                                                                                                          <  8:237667/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                                                                                     >  2:110778/1‑145 (MQ=37)                                                            ccATGTGCCCTCGCCTGACGAGCGGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGc                                                                                     <  1:110778/145‑1 (MQ=37)                                                              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)                                                                               gAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATctgct                                                              >  2:56156/1‑149 (MQ=37)                                                                                 gcagcCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATgg                                                                                                     >  4:135087/1‑108 (MQ=40)                                                                                    gccgccCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCa                                                         >  2:291380/1‑149 (MQ=37)                                                                                        ccccTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCAc                                                     <  2:340667/149‑1 (MQ=37)                                                                                                cGACAACGCCCCGCCGACACGGACGACTTCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTAACGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTca                                                >  6:378197/1‑146 (MQ=255)                                                                                                 gACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTAc                                                                                        >  4:372828/1‑105 (MQ=35)                                                                                                 gACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTAc                                                                                        <  3:372828/105‑1 (MQ=35)                                                                                                 gACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGa                                            <  1:56156/149‑1 (MQ=37)                                                                                                 gACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTGCACAGa                                            >  4:251916/1‑149 (MQ=37)                                                                                                   cAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACa                                          <  4:230400/149‑1 (MQ=35)                                                                                                         cccGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGCCGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAg                                    >  5:401314/1‑149 (MQ=25)                                                                                                         cccGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAg                                    <  7:241224/149‑1 (MQ=35)                                                                                                          ccgccgACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAgg                                   <  3:284219/149‑1 (MQ=37)                                                                                                           cgccgACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGga                                  <  7:319611/149‑1 (MQ=38)                                                                                                                acacGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGCCAATAAAGGAGaaa                              >  3:20846/1‑148 (MQ=255)                                                                                                                  acGGACGACTACCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAAc                           <  4:373557/149‑1 (MQ=255)                                                                                                                    ggACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCCTCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAg                         <  5:397650/149‑1 (MQ=255)                                                                                                                    ggACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATGACCTTCACAGACAATAAAGGAGAAAAACAg                         >  1:305111/1‑149 (MQ=255)                                                                                                                      acgacTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCa                       <  5:43154/149‑1 (MQ=39)                                                                                                                           tGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTc                   >  6:114832/1‑148 (MQ=37)                                                                                                                           tGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTc                   >  6:68093/1‑148 (MQ=37)                                                                                                                           tGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCg                  >  6:109278/1‑149 (MQ=37)                                                                                                                             ccATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGAc                <  1:43788/149‑1 (MQ=37)                                                                                                                              cATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTc                   <  4:114148/145‑1 (MQ=37)                                                                                                                              cATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTc                   >  3:114148/1‑145 (MQ=37)                                                                                                                              cATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTAAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACCACGCGGCCGACAATCTCCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGAc                >  1:305950/1‑148 (MQ=255)                                                                                                                               aTCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTa              >  7:359276/1‑149 (MQ=37)                                                                                                                               aTCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTa              >  3:343389/1‑149 (MQ=37)                                                                                                                                tCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTa              >  6:111334/1‑148 (MQ=37)                                                                                                                                  cGCCATGACTTCTTTCTGTACGCCTTCGGCCTGCTTGTCAGCAAGCGCAACGAAGATGTCACCACCGAATAACGCTACTACGCGGCCGACAATCTGCTGTCCATCACCTTAACATACAATAAAGGAGAAAAACAGCAACTCGACTacac           <  8:359276/149‑1 (MQ=11)                                                                                                                                  cGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTacac           <  3:147883/149‑1 (MQ=37)                                                                                                                                           ttCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAg                                    <  6:244285/115‑1 (MQ=34)                                                                                                                                           ttCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAg                                    >  5:244285/1‑115 (MQ=34)                                                                                                                                           ttCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTACACCTACGACGc  >  5:115891/1‑149 (MQ=37)                                                                                                                                               |                                                                                                                                                 CTGGATGGCAGTGGCCGCCTCTATGCCTACAACGTGCTGGATGAGGTCATTCGTCTTACCCATGTGCCCTCGCCTGACGAGCAGCCGCCCCTGTCCGACAACGCCCCGCCGACACGGACGACTGCCATCCGCCATGACTTCGAGCGGGACGCCATCGGCCGGCTTGTCAGCAAGCGCACCGAAGATGGCACCACCGACTACCGCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAATAAAGGAGAAAAACAGCAACTCGACTACACCTACGACGC  >  NC_002947/4616243‑4616529

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