Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,551,221	G→C	6.5%	G125R (GGC→CGC)	PP_3135 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,551,221	0	G	C	6.5%	178.8 / 2.1	57	G125R (GGC→CGC)	PP_3135	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base C (2/2);  total (28/29)
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

GCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGACTGGTT  >  NC_002947/3551083‑3551366                                                                                                                                               |                                                                                                                                                     gcgcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGc                                                                                                                                             <  4:33412/148‑1 (MQ=255) gcgcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCCATCGCGAccctcc                                                                                                                                            <  8:21568/149‑5 (MQ=255)  cgcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAg                                                                                                                                           <  8:212788/149‑1 (MQ=255)   gcACCTCAACCGAAATTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCCATCGCGAccctccag                                                                                                                                          <  8:308613/149‑7 (MQ=255)        tCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGt                                                                                                                                     >  8:110291/1‑149 (MQ=255)          aaGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCGGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGtt                                                                                                                                    <  6:193442/148‑1 (MQ=255)              gCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGc                                                                                                                               <  7:18762/149‑1 (MQ=255)                aaGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTa                                                                                                                             <  6:264191/149‑1 (MQ=255)                            ctgctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACAcc                                                                                                                 <  2:64748/149‑1 (MQ=255)                              gctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGACAACCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGAGCCTGGAGGTACGCGATGGCGACCTCCAGCCAAGTTGACGCTAATTCCCCCCACCg                                                                                                                >  1:284839/1‑149 (MQ=255)                               ctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                              >  7:149423/1‑149 (MQ=255)                               ctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                              >  5:353201/1‑149 (MQ=255)                               ctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                              >  4:17991/1‑149 (MQ=255)                                   aTACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAAt                                                                                                          <  6:353201/149‑1 (MQ=255)                                       cGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATgtgt                                                                                                      >  3:339976/1‑149 (MQ=255)                                       cGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATgtgt                                                                                                      >  8:13393/1‑149 (MQ=255)                                         aTGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTtc                                                                                                    <  4:283838/149‑1 (MQ=255)                                             cTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGGAGTACATGGAC‑GCCCACCCGGTGTGCGGCGTACT‑‑GGGGGTGCGCGTGGAGGTGCGCGATGGCGACCTGCCGCCCAGTTGTGGGTAGTTCCCACCACCGGTCATTGTGttgggc                                                                                                 >  1:231181/1‑145 (MQ=255)                                              ttCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCg                                                                                                                                <  6:46469/116‑1 (MQ=255)                                              ttCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCg                                                                                                                                >  5:46469/1‑116 (MQ=255)                                               tCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCggg                                                                                              >  6:249661/1‑149 (MQ=255)                                                  ttGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGt                                                                                           >  1:181221/1‑149 (MQ=255)                                                       gccgACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCgg                                                                                      >  4:319220/1‑149 (MQ=255)                                                          gACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACTGGGGGGGTGCGCCTGGAGGGTCGCGATGGCGACATGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCTTCGGGCGTACCggg                                                                                     >  1:285238/1‑149 (MQ=255)                                                          gACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGtt                                                                                   <  3:17991/149‑1 (MQ=255)                                                           aCTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTg                                                                                  <  5:249661/149‑1 (MQ=255)                                                               gCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTgggg                                                                               <  1:311958/148‑1 (MQ=255)                                                                cTGCAAAAAACGCTC‑GACTACAGGGAC‑GCCCCCCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTCGAGCGTCGCGATCGCGCCCTGCAGCCCAGGGGTCCCTACTGCCACCCACCGCGCAATGGGATCGTCTGGCCTACCCGGTGGGGGCg                                                                             >  2:216012/1‑149 (MQ=255)                                                                                gACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGt                                                              >  3:64796/1‑149 (MQ=255)                                                                                     cATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑Gaa                                                          <  7:110291/147‑1 (MQ=255)                                                                                     cATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAg                                                         <  2:181221/148‑1 (MQ=255)                                                                                      aTGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGAGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑t                                                     >  8:277918/1‑149 (MQ=255)                                                                                             gCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTcga                                               <  2:285238/149‑1 (MQ=255)                                                                                             gCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTcga                                               <  4:11516/149‑1 (MQ=255)                                                                                              cccACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTcgac                                              <  3:319220/149‑1 (MQ=255)                                                                                              cccACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTcgac                                              >  6:152349/1‑149 (MQ=255)                                                                                               ccACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGATTCTTCCCAGGGGT‑GAAG‑A‑TGGTc                                                 <  6:156862/145‑1 (MQ=255)                                                                                               ccACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGATTCTTCCCAGGGGT‑GAAG‑A‑TGGTc                                                 >  5:156862/1‑145 (MQ=255)                                                                                                  cccGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACa                                          <  5:152349/149‑1 (MQ=255)                                                                                                             gcgAACT‑‑GGGGCTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTGGTCGCTACTTCCCCACACCGCTCAATCTGATCGTCGGGCCGACCGGGTTGGGGCCGGTCTTCCCAGGGGG‑GAAG‑A‑TGGTCGACGACATGACGTgggcc                               >  2:183666/1‑147 (MQ=255)                                                                                                               gTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGttctt                                                                       >  6:34909/1‑110 (MQ=255)                                                                                                               gTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGttctt                                                                       <  5:34909/110‑1 (MQ=255)                                                                                                                       ggggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGcc                         <  1:18661/147‑1 (MQ=255)                                                                                                                       ggggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCt                        <  4:339976/148‑1 (MQ=255)                                                                                                                        gggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCGCACACCGCTCAATGGGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCGCAGGGGT‑GAAG‑A‑TGGTCGACGGCATGAGCTGGGCCCACGCGTCg                      >  6:264305/1‑149 (MQ=255)                                                                                                                        gggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGAcc                              >  5:174974/1‑141 (MQ=255)                                                                                                                        gggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGAcc                              <  6:174974/141‑1 (MQ=255)                                                                                                                              gcCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTgcgc                >  5:221875/1‑149 (MQ=255)                                                                                                                                  ggAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAgt            >  2:170503/1‑149 (MQ=255)                                                                                                                                   gcaggTCGCCATCGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAAGTGGGACCACGCCTCGGTGGGGGAgtg           >  4:193378/4‑149 (MQ=255)                                                                                                                                   gAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTGGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGGCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAgt            <  7:277918/148‑1 (MQ=255)                                                                                                                                      ggTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGGGAGTAGTGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAgtg           >  1:119956/1‑149 (MQ=255)                                                                                                                                      ggTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGa        <  7:13393/149‑1 (MQ=255)                                                                                                                                       gTCGCGATGGCGCCCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTGCGTCGGGCGTACCGGGTTGGGGGGGTTCTTCGCAGGGGT‑GTAG‑A‑TGGGCGACGACATGAGCTGGGCCCTCGCGTCGGTGGGCGAGTGTGAc       >  8:130041/1‑149 (MQ=255)                                                                                                                                            gATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGACTGGtt  <  1:141161/149‑1 (MQ=255)                                                                                                                                              tGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTgc                  <  7:35954/131‑1 (MQ=255)                                                                                                                                              tGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTgc                  >  8:35954/1‑131 (MQ=255)                                                                                                                                               |                                                                                                                                                     GCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTC‑GACTACATGGAC‑GCCCACCCGGAGTGCGGCGTACT‑‑GGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGT‑GAAG‑A‑TGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGACTGGTT  >  NC_002947/3551083‑3551366

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