Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,872,138	A→T	10.4%	H81Q (CAT→CAA)	PP_3420 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,872,138	0	A	T	10.4%	181.6 / 10.5	67	H81Q (CAT→CAA)	PP_3420	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base T (4/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATCAC  >  NC_002947/3871993‑3872281                                                                                                                                                  |                                                                                                                                                ccAGCAAGGTGTCCGGCAGGTGCTCGACATTATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGACGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGc                                                                                                                                              <  7:17479/149‑1 (MQ=255)  cAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGcc                                                                                                                                             <  7:267641/149‑1 (MQ=255)      aaGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCATTCGCAATTGCCGcct                                                                                                                                         >  4:175730/1‑147 (MQ=255)       aGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCATTCGCAATTGCCGcctg                                                                                                                                        <  3:175730/149‑4 (MQ=255)                        tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGc                                                                                                                             >  4:174187/1‑143 (MQ=255)                        tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGc                                                                                                                             <  3:174187/143‑1 (MQ=255)                        tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                       >  8:195911/1‑149 (MQ=255)                        tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                       >  4:39432/1‑149 (MQ=255)                          gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                       <  1:183283/147‑1 (MQ=255)                          gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                       <  8:32910/147‑1 (MQ=255)                          gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                       <  3:39432/147‑1 (MQ=255)                          gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCt                                                                                                                     <  7:272434/149‑1 (MQ=255)                             aTGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCtgaa                                                                                                                  <  1:147619/149‑1 (MQ=255)                             aTGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCtgaa                                                                                                                  <  1:190455/149‑1 (MQ=255)                              tGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCtgaat                                                                                                                 <  8:86609/149‑1 (MQ=255)                                          tCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTg                                                                                                           <  1:75591/143‑1 (MQ=255)                                          tCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTg                                                                                                           >  2:75591/1‑143 (MQ=255)                                            gATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGt                                                                                                     <  4:303980/147‑1 (MQ=255)                                            gATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGt                                                                                                     >  4:124108/1‑147 (MQ=255)                                                cGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCAt                                                                                               >  1:80867/1‑149 (MQ=255)                                                 gggTCGATTCGCCCACCGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCCCCGCTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGCCCATTCCCTGAATGAACTGCGCGGTAAGCATg                                                                                              >  4:165941/1‑149 (MQ=255)                                                           ggCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGACAGGCGGCAATTGCGAATGCCGCACTTCGGACTGGGACATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGctct                                                                                    >  8:177139/1‑149 (MQ=255)                                                            gCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCAtt                                                                                                                         >  5:28298/1‑111 (MQ=255)                                                            gCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCAtt                                                                                                                         <  6:28298/111‑1 (MQ=255)                                                             ccAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTcc                                                                                  <  2:245240/149‑1 (MQ=255)                                                                  gCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCGGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCt                                                                             <  7:195911/149‑1 (MQ=255)                                                                          gACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGca                                                                     >  2:132843/1‑149 (MQ=255)                                                                             cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCgcgc                                                                      >  3:353205/1‑145 (MQ=255)                                                                             cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCgcgc                                                                      <  4:353205/145‑1 (MQ=255)                                                                             cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcccca                                                                  >  6:341414/1‑149 (MQ=255)                                                                             cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                  >  7:98083/1‑149 (MQ=255)                                                                             cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                  >  8:190083/1‑149 (MQ=255)                                                                             cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                  >  2:377430/1‑149 (MQ=255)                                                                                cAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGGGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACg                                                                >  1:133416/1‑148 (MQ=255)                                                                                  aGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAg                                                             <  8:47239/149‑1 (MQ=255)                                                                                  aGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAg                                                             <  7:177139/149‑1 (MQ=255)                                                                                   gTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacacagtat                                                             >  3:271481/1‑141 (MQ=255)                                                                                        gcAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCTCGCTGTCGCGCACCACGAATCGCGTc                                                       >  8:173358/1‑149 (MQ=255)                                                                                             tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCGAGGCTGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATTAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                  <  6:237534/149‑1 (MQ=255)                                                                                             tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                  <  1:132843/149‑1 (MQ=255)                                                                                             tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                  <  1:188373/149‑1 (MQ=255)                                                                                             tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                  <  7:190083/149‑1 (MQ=255)                                                                                                    gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAgcg                                                           >  4:208623/1‑133 (MQ=255)                                                                                                    gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAgcg                                                           <  3:208623/133‑1 (MQ=255)                                                                                                    gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAAcc                                            >  3:237006/1‑148 (MQ=255)                                                                                                            aGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCAcc                                   >  5:329130/1‑149 (MQ=255)                                                                                                            aGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCAcc                                   >  3:110921/1‑149 (MQ=255)                                                                                                                aCTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAg                                         <  3:33983/139‑1 (MQ=255)                                                                                                                aCTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAg                                         >  4:33983/1‑139 (MQ=255)                                                                                                                     cccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCAcc                            >  1:25985/1‑147 (MQ=255)                                                                                                                      ccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGt                                                        <  6:86536/118‑1 (MQ=255)                                                                                                                      ccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGt                                                        >  5:86536/1‑118 (MQ=255)                                                                                                                           gtggtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTaa                     <  4:345452/148‑1 (MQ=255)                                                                                                                            tggtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATg                   <  4:237006/149‑1 (MQ=255)                                                                                                                              gtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGAt                 <  4:110921/149‑1 (MQ=255)                                                                                                                               cgaagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGAt                 >  4:82833/7‑148 (MQ=255)                                                                                                                                 ccAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa              <  2:25985/149‑1 (MQ=255)                                                                                                                                 aagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGa                  <  1:294575/141‑1 (MQ=255)                                                                                                                                 aagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGGCGAAACGTAACCCAGCCCCCCTGCCACCGCGGTAATGa                  >  2:294575/5‑145 (MQ=255)                                                                                                                                   gtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAgg                        <  3:193776/135‑1 (MQ=255)                                                                                                                                   gtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAgg                        >  4:193776/3‑137 (MQ=255)                                                                                                                                   aGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa              >  8:241557/1‑147 (MQ=255)                                                                                                                                   aGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa              >  5:92026/1‑147 (MQ=255)                                                                                                                                     gcgggAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGAGCGACGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAgg          >  7:158551/1‑149 (MQ=255)                                                                                                                                          aaTTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGAATGGCATCGCGCTCGACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGAt     >  5:259243/1‑149 (MQ=255)                                                                                                                                           aTTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATc    <  3:124108/149‑1 (MQ=255)                                                                                                                                              gCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATCAc  <  6:92026/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                CCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATCAC  >  NC_002947/3871993‑3872281

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