Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,981,954	T→G	5.1%	S867A (TCA→GCA)	PP_5242 →	GAF domain‑containing sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,981,954	0	T	G	5.1%	276.3 / 2.1	88	S867A (TCA→GCA)	PP_5242	GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (46/37);  new base G (3/2);  total (49/39)
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

TGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTGCACGGCGTTGCCCCTGG  >  NC_002947/5981806‑5982102                                                                                                                                                       |                                                                                                                                                     tGACTATCAGGCGACCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGt                                                                                                                                                      >  4:238354/1‑149 (MQ=255)     tATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACt                                                                                                                                                  >  3:381821/1‑149 (MQ=255)     tATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACt                                                                                                                                                  >  2:332053/1‑149 (MQ=255)              gcCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCAc                                                                                                                                         <  1:401801/149‑1 (MQ=255)                cATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCAcc                                                                                                                                        >  8:290705/1‑148 (MQ=255)                 aTCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGa                                                                                                                                      >  4:340819/1‑149 (MQ=255)                  tCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTACCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCa                                                                                                                                                    >  6:17931/1‑134 (MQ=255)                  tCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTACCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCa                                                                                                                                                    <  5:17931/134‑1 (MQ=255)                    gACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACt                                                                                                                                    <  7:332304/148‑1 (MQ=255)                    gACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGC‑GGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCg                                                                                                                                  <  5:63232/149‑1 (MQ=255)                       aGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGTCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCgga                                                                                                                                <  1:322833/149‑1 (MQ=255)                       aGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCgga                                                                                                                                <  3:238354/149‑1 (MQ=255)                            ccTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCAc                                                                                                                                                   <  8:261534/125‑1 (MQ=255)                            ccTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCAc                                                                                                                                                   >  7:261534/1‑125 (MQ=255)                            ccTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCggagg                                                                                                                              >  7:72793/1‑146 (MQ=255)                            ccTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGaa                                                                                                                            <  8:72793/148‑1 (MQ=255)                            ccTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAg                                                                                                                           <  7:290705/149‑1 (MQ=255)                                gccggccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGa                                                                                                                                      <  2:29994/134‑1 (MQ=255)                                gccggccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGa                                                                                                                                      >  1:29994/1‑134 (MQ=255)                                   ggccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGt                                                                                                                    >  2:351731/1‑149 (MQ=255)                                   ggccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑GTGCCTG‑GAGTCCCTGCCGCGCACCGCCTCGGAGGAAGTGACGGt                                                                                                                    >  2:387260/1‑149 (MQ=255)                                     ccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGtt                                                                                                                   >  6:42947/1‑148 (MQ=255)                                     ccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGtt                                                                                                                   <  5:42947/148‑1 (MQ=255)                                     ccTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTg                                                                                                                  <  5:111260/149‑1 (MQ=255)                                      cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGTCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGc                                                                                                                 <  5:266272/149‑1 (MQ=255)                                      cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTg                                                                                                                  >  7:386209/1‑148 (MQ=255)                                      cTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGc                                                                                                                 >  4:412386/1‑149 (MQ=255)                                           cTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCa                                                                                                            <  2:375806/149‑1 (MQ=255)                                            ttACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAg                                                                                                           >  5:366736/1‑149 (MQ=255)                                                 gggTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGc                                                                                                                                                >  6:306253/1‑107 (MQ=255)                                                 gggTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGc                                                                                                                                                <  5:306253/107‑1 (MQ=255)                                                  ggTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCctgct                                                                                                     <  4:381821/149‑1 (MQ=255)                                                  ggTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCctgct                                                                                                     >  1:409575/1‑149 (MQ=255)                                                   gTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTGTGAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCctgct                                                                                                     <  1:332053/149‑1 (MQ=255)                                                        aGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTTATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGc                                                                                                  >  4:368594/1‑146 (MQ=255)                                                        aGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTTATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGc                                                                                                  <  3:368594/146‑1 (MQ=255)                                                        aGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCAc                                                                                               >  2:351089/1‑149 (MQ=255)                                                             ggTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGt                                                                                          >  1:299005/1‑149 (MQ=255)                                                              gTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGtc                                                                                         <  2:44680/149‑1 (MQ=255)                                                                          gACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCcttctt                                                                             >  2:2496/1‑149 (MQ=255)                                                                          gACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCcttctt                                                                             <  2:299005/149‑1 (MQ=255)                                                                                cagcagGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGaa                                                                       <  1:351089/149‑1 (MQ=255)                                                                                cagcagGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGaa                                                                       <  8:187915/149‑1 (MQ=255)                                                                                     ggTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCAGTTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGgac                                                                   >  1:114221/1‑149 (MQ=255)                                                                                      gTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGgacga                                                                 >  3:350388/1‑149 (MQ=255)                                                                                                 tCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGc                                                                                                          <  3:154719/97‑1 (MQ=255)                                                                                                 tCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGc                                                                                                          >  4:154719/1‑97 (MQ=255)                                                                                                   aTGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGtt                                                                                                                   <  8:291221/86‑1 (MQ=255)                                                                                                   aTGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGtt                                                                                                                   >  7:291221/1‑86 (MQ=255)                                                                                                     gCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCa                                                  >  3:308381/1‑149 (MQ=255)                                                                                                       ccATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCctgct                                                                                                     >  4:45288/1‑96 (MQ=255)                                                                                                       ccATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCctgct                                                                                                     <  3:45288/96‑1 (MQ=255)                                                                                                       ccATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCttcttc                                                                            >  8:298102/1‑121 (MQ=255)                                                                                                       ccATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCttcttc                                                                            <  7:298102/121‑1 (MQ=255)                                                                                                          tgtgCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTgg                                             <  6:366736/149‑1 (MQ=255)                                                                                                            tgCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGcc                                           <  3:340819/149‑1 (MQ=255)                                                                                                             gCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGcc                                           >  6:112699/1‑148 (MQ=255)                                                                                                              cAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCGACTCTGTCATCATCGTTGTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCaa                                         <  8:386209/149‑1 (MQ=255)                                                                                                              cAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCGCCTGGcc                                           >  6:158164/1‑147 (MQ=255)                                                                                                                 gCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCt                                      >  1:169454/1‑149 (MQ=255)                                                                                                                   gTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCa                                    >  3:264365/1‑149 (MQ=255)                                                                                                                      gATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTg                                                    >  4:259088/1‑130 (MQ=255)                                                                                                                      gATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTg                                                    <  3:259088/130‑1 (MQ=255)                                                                                                                       aTCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATgg                                 <  6:416451/148‑1 (MQ=255)                                                                                                                       aTCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGa                                <  3:412386/149‑1 (MQ=255)                                                                                                                        tCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAg                               >  3:67775/1‑149 (MQ=255)                                                                                                                        tCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGCAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCCTCATCCTTCTTCGTGAAGGACGACATCAGATTCATGCACCTGGCCAAGCTCCTGGAg                               >  2:279556/1‑149 (MQ=255)                                                                                                                          gAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTC‑GAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGtt                             <  4:208364/148‑1 (MQ=255)                                                                                                                           aaCAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTg                            >  7:335195/1‑149 (MQ=255)                                                                                                                           aaCAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTg                            >  4:53483/1‑149 (MQ=255)                                                                                                                            tgaGAATGACCTGCT‑TTGCCTG‑GATGCACTGCCGCGCACCGACTCGGATGAAGTGACTGTTGGCATCAGCCTGCTGGCCACTCTGTGATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGt                           <  1:149682/147‑1 (MQ=255)                                                                                                                              aGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGgacg                                                                  <  8:30415/108‑1 (MQ=255)                                                                                                                              aGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGgacg                                                                  >  7:30415/1‑108 (MQ=255)                                                                                                                              aGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTc                          >  2:74718/1‑148 (MQ=255)                                                                                                                              aGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTc                          <  1:74718/148‑1 (MQ=255)                                                                                                                              aGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCg                         >  1:23541/1‑149 (MQ=255)                                                                                                                                aaGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGc                        >  3:424516/1‑148 (MQ=255)                                                                                                                                aaGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGc                        >  2:395163/1‑148 (MQ=255)                                                                                                                                   gcgcGGCA‑GTGACTC‑CAGGCAATGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACt                    >  5:143161/4‑149 (MQ=255)                                                                                                                                   gACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACt                    >  5:126445/1‑149 (MQ=255)                                                                                                                                   gACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACt                    <  2:409575/149‑1 (MQ=255)                                                                                                                                     gcGGCA‑GTGACTC‑CAGGCAATGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGc                        >  3:330081/2‑143 (MQ=255)                                                                                                                                     gcGGCA‑GTGACTC‑CAGGCAATGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGc                        <  4:330081/142‑1 (MQ=255)                                                                                                                                     gcGGCA‑GTGACTC‑CAGGCAATGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTg                   >  5:275843/2‑148 (MQ=255)                                                                                                                                     ccGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACgg                                                                                                                     <  4:299444/50‑1 (MQ=255)                                                                                                                                     ccGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACgg                                                                                                                     >  3:299444/1‑50 (MQ=255)                                                                                                                                      cGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTGCa                 >  4:350520/1‑149 (MQ=255)                                                                                                                                            ttGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTGCACGGCg            >  3:94900/1‑149 (MQ=255)                                                                                                                                                  g‑gAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTGCACGGCGTTGccc      >  2:165468/1‑149 (MQ=255)                                                                                                                                                       tCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTGCACGGCGTTGCCCCTgg  <  2:114221/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     TGACTATCAGGCGGCCATCGACAGCGCCCTGGCCGGCCTGGACTTACTGGGTATCAGCCCGGTGCGGGGCCTCGACTGGCAGCAGGTCGAGGCCAGTCATGCCCATGTGCAGGCGTTGATCGAACAGAAGGACCGGCA‑TTGCCTG‑GAGTCACTGCCGCGCACCGACTCGGAGGAAGTGACGGTTGCCATCAGCCTGCTGGCCACTCTGTCATCATCCTTCTTCGTGAAGGACGACATCAGCTTCCTGCACCTGGCCAAGCTCATGGAGTTGTCGCTACTGCACGGCGTTGCCCCTGG  >  NC_002947/5981806‑5982102

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