Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	590,703:1	+A	6.7%	coding (827/885 nt)	PP_0500 →	dTDP‑4‑rhamnose reductase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	590,703	1	.	A	6.7%	237.9 / 11.4	60	coding (827/885 nt)	PP_0500	dTDP‑4‑rhamnose reductase‑like protein
Reads supporting (aligned to +/- strand):  ref base . (26/30);  new base A (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GATCCTCGCCGAGGCCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATT  >  NC_002947/590557‑590849                                                                                                                                                    |                                                                                                                                                    gATCCTCGCCGAGGCCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCc‑gc                                                                                                                                                   >  6:349495/1‑149 (MQ=255)   tCCTCGCCGAGGCCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGt                                                                                                                                                 <  2:234909/149‑1 (MQ=255)               ccggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCt                                                                                                                                     <  4:136588/149‑1 (MQ=255)                cggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGgcacgcg                                                                                                                                     >  3:98639/1‑148 (MQ=255)                cggccAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTg                                                                                                                                    <  3:339974/149‑1 (MQ=255)                      aTACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGc                                                                                                                              <  6:330022/149‑1 (MQ=255)                       tACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGcc                                                                                                                             <  3:77470/149‑1 (MQ=255)                        aCCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCt                                                                                                                                             <  7:102655/132‑1 (MQ=255)                        aCCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCt                                                                                                                                             >  8:102655/1‑132 (MQ=255)                         ccGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCCAGCGCGCCAGGCACGC‑GGCTTgata                                                                                                                            <  4:98639/148‑4 (MQ=255)                         ccGCAAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                           >  5:184638/1‑149 (MQ=255)                          cGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAcc                                                                                                                          >  7:31847/1‑149 (MQ=255)                          cGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAcc                                                                                                                          >  7:37188/1‑149 (MQ=255)                             cAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACct                                                                                                                         <  8:332309/147‑1 (MQ=255)                              aGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTa                                                                                                                      >  4:139601/1‑149 (MQ=255)                               gCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                     <  3:139601/149‑1 (MQ=255)                               gCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTAc                                                                                                                     <  2:316253/149‑1 (MQ=255)                                    ccGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGAc                                                                                                                >  1:170527/1‑149 (MQ=255)                                    ccGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGAc                                                                                                                >  8:306515/1‑149 (MQ=255)                                     cGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCa                                                                                                                            >  6:211079/1‑136 (MQ=255)                                     cGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCa                                                                                                                            <  5:211079/136‑1 (MQ=255)                                      gTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTgg                                                                                                                  >  4:170492/1‑145 (MQ=255)                                      gTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTgg                                                                                                                  <  3:170492/145‑1 (MQ=255)                                       ttCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCgg                                                                                                             >  4:190644/1‑149 (MQ=255)                                            gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTa                                                                                                       <  5:349495/149‑1 (MQ=255)                                            gCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTa                                                                                                       <  6:124670/149‑1 (MQ=255)                                             cccccACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCAc                                                                                                                           >  6:186191/1‑129 (MQ=255)                                                ccACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGc                                                                                                   <  2:170527/149‑1 (MQ=255)                                                 cACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGc                                                                                                   >  5:224074/1‑148 (MQ=255)                                                 cACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGc                                                                                                   >  7:236551/1‑148 (MQ=255)                                                   cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCAt                                                                                                <  5:129357/149‑1 (MQ=255)                                                   cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCAt                                                                                                <  3:190644/149‑1 (MQ=255)                                                   cgccgcAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCAt                                                                                                <  8:236551/149‑1 (MQ=255)                                                           gCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACg                                                                                        >  3:69146/1‑149 (MQ=255)                                                            cccATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGc                                                                                       <  8:31847/149‑1 (MQ=255)                                                            cccATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGc                                                                                       <  8:37188/149‑1 (MQ=255)                                                                          caggATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCa                                                                         >  2:266550/3‑149 (MQ=255)                                                                               tGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACcggc                                                                    <  6:224074/149‑1 (MQ=255)                                                                                              acaGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                     >  5:23392/1‑149 (MQ=255)                                                                                               caGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTTCCTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                     <  1:266550/149‑1 (MQ=255)                                                                                                   aCGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTc                                                <  4:69146/149‑1 (MQ=255)                                                                                                    cgcgGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTcc                                               >  6:320932/1‑149 (MQ=255)                                                                                                            tGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACcg                                                                      <  6:169651/118‑1 (MQ=255)                                                                                                            tGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACcg                                                                      >  5:169651/1‑118 (MQ=255)                                                                                                                  gCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                     >  3:49166/1‑129 (MQ=255)                                                                                                                  gCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAtt                                                     <  4:49166/129‑1 (MQ=255)                                                                                                                       aaaaTCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCgttgt                            <  5:320932/149‑1 (MQ=255)                                                                                                                           tCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGa                        >  4:284159/1‑149 (MQ=255)                                                                                                                              ttCATACCTTCGGTATCAAGCC‑GCGTGCCGGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCAt                                                      <  6:248709/116‑1 (MQ=255)                                                                                                                                cATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                    <  6:184638/148‑1 (MQ=255)                                                                                                                                       tCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCt                                          >  3:39861/1‑119 (MQ=255)                                                                                                                                       tCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCt                                          <  4:39861/119‑1 (MQ=255)                                                                                                                                       tCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGc            >  6:363587/1‑149 (MQ=255)                                                                                                                                         ggTATCTAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTG‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTCGACAAAGGCTACGcc           <  6:192725/148‑1 (MQ=255)                                                                                                                                           tATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTa        >  1:272088/1‑149 (MQ=255)                                                                                                                                           tATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTa        >  5:29205/1‑149 (MQ=255)                                                                                                                                             tCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTCATCACCGCCGGTGCCGGCTTCATTGTCTCCCACCTGTGCCATGCGTTGTTGGAc                       >  4:50798/1‑132 (MQ=255)                                                                                                                                             tCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCAACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGAc                       <  3:50798/132‑1 (MQ=255)                                                                                                                                             gcAAGCCAGCGCGCCAGGCACGC‑GGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                    >  5:237632/2‑135 (MQ=255)                                                                                                                                             gcAAGCCAGCGCGCCAGGCACGC‑GGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACaaa                    <  6:237632/134‑1 (MQ=255)                                                                                                                                                 gCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCAtt  >  7:171571/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GATCCTCGCCGAGGCCGGCCAATACCGCCAGCTGGCCGTTCAGGCCCCCACGCCGCAGGCCCATGCCGCGCGCCCGGATGCCAGCGAAGAGCCACAGCACGCGGTACTGGCCTGCAAGAAAATCCTTCATACCTTCGGTATCAAGCC‑GCGTGCCTGGCGCGCTGGCTTGCCACCTCTACTGGACCGGTT‑CTACCGCCATGGCTGACGCCCCCATCCTGATCACCGGCGGTGCCGGCTTCATTGGCTCCCACCTGTGCGATGCGTTGTTGGACAAAGGCTACGCCGTACGCATT  >  NC_002947/590557‑590849

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