Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,619,471	A→G	8.9%	C344R (TGT→CGT)	PP_4938 ←	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,619,471	0	A	G	8.9%	183.5 / 7.4	67	C344R (TGT→CGT)	PP_4938	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/4);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAGGC  >  NC_002947/5619327‑5619617                                                                                                                                                  |                                                                                                                                                   cgccTTCGAGGCAGCCTTGCTCACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGCTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGCTGCACTGAGCGCCGTACcgcg                                                                                                                                                 >  4:38047/1‑148 (MQ=255)   ccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTc                                                                                                                                              <  8:2481/149‑1 (MQ=255)   ccTTCGAGCCAGCCTTGCTGACACGCTGAAGGTCAGATGCGCACCAGGCTCTGCGACCCGATGCAGAAATCCTCCCCGGCCTGCTCAGGTGCCAGTGACAGG‑CGCTCCTTAGAGGCGAGTGTCGCTGCACTGAGCTCCCTACAGCtctc                                                                                                                                              >  4:122210/1‑149 (MQ=255)     ttCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGa                                                                                                                                            >  7:123215/1‑149 (MQ=255)       cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                           <  5:17539/148‑4 (MQ=255)        gAGGCATCCTTGCTGACACGCTGAACGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGTCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                           <  1:75531/148‑4 (MQ=255)                cTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGcc                                                                                                                                 >  4:75991/1‑149 (MQ=255)                   gCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGccc                                                                                                                                <  6:241625/147‑1 (MQ=255)                      gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGgcg                                                                                                                           >  5:153963/1‑149 (MQ=255)                       acacGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGgc                                                                                                                            <  5:5874/147‑1 (MQ=255)                             tGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTca                                                                                                                                             >  7:173476/1‑123 (MQ=255)                             tGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTca                                                                                                                                             <  8:173476/124‑2 (MQ=255)                              gAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTACAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCTCTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTtcg                                                                                                                     <  1:217639/147‑1 (MQ=255)                              gAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTACAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTtcg                                                                                                                     >  2:217639/1‑147 (MQ=255)                               aaGGTCAGCTGCGCAACAGGCTCTGCGAGTCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTtcgtcg                                                                                                                  <  3:253999/149‑1 (MQ=255)                               aaGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTtcgtcg                                                                                                                  <  1:123599/149‑1 (MQ=255)                                     aGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAAc                                                                                                             <  3:181277/148‑1 (MQ=255)                                         gcgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCgg                                                                                                          >  2:241883/1‑147 (MQ=255)                                          cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                       <  3:75991/149‑1 (MQ=255)                                             aaCAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTa                                                                                                    <  5:241696/149‑1 (MQ=255)                                                        cGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCgtgt                                                                                         >  7:101644/1‑149 (MQ=255)                                                        cGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCg                                                                                            <  5:48625/146‑1 (MQ=255)                                                           gCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgc                                                                                      <  1:241883/149‑1 (MQ=255)                                                            ccGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgcg                                                                                     <  6:73263/149‑1 (MQ=255)                                                            ccGCTGCAGAACGCCTGCCAGGACTGCTCCGGTGCCAGTGCCTGG‑CGATGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTGTCGTCGGCACCCGGCAGGTAGCCACCCGTGTcgcg                                                                                     >  5:88894/1‑149 (MQ=255)                                                                  cAGAACGCCTGCCAGGCCTGCTCAGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACACGTGTCGCGCAACTg                                                                               >  8:236511/1‑149 (MQ=255)                                                                   aGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTgg                                                                              >  1:263643/1‑148 (MQ=255)                                                                   aGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGc                                                                              >  7:245694/1‑149 (MQ=255)                                                                        gCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgc                                                                                      <  4:150743/136‑1 (MQ=255)                                                                        gCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgc                                                                                      >  3:150743/1‑136 (MQ=255)                                                                         ccTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGa                                                                        <  6:153963/149‑1 (MQ=255)                                                                         ccTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGa                                                                        <  8:245694/149‑1 (MQ=255)                                                                         ccTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGa                                                                        <  3:284555/149‑1 (MQ=255)                                                                           tGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGaaa                                                                      >  1:23238/1‑149 (MQ=255)                                                                           tGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGaaa                                                                      >  4:166655/1‑149 (MQ=255)                                                                            gCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGaaa                                                                      >  6:67576/1‑148 (MQ=255)                                                                             ccAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGc                                                                                                         >  5:94535/1‑112 (MQ=255)                                                                             ccAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGc                                                                                                         <  6:94535/112‑1 (MQ=255)                                                                                 gCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTc                                                                <  8:123215/149‑1 (MQ=255)                                                                                    tGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAg                                                             >  7:113567/1‑149 (MQ=255)                                                                                         cGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGc                                                        >  8:123702/1‑149 (MQ=255)                                                                                         cGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGc                                                        >  4:262842/1‑149 (MQ=255)                                                                                           gTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTg                                                      <  7:120479/149‑1 (MQ=255)                                                                                           gTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTg                                                      >  6:201725/1‑149 (MQ=255)                                                                                                aGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGCCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCacac                                                 >  3:70070/1‑149 (MQ=255)                                                                                                     cTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCg                                            <  3:166655/149‑1 (MQ=255)                                                                                                     cTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCg                                            <  8:101644/149‑1 (MQ=255)                                                                                                     cTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCg                                            >  1:145925/1‑149 (MQ=255)                                                                                                            ctgctTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACAt                                      >  1:141615/1‑149 (MQ=255)                                                                                                              gctTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgcg                                                                                     >  6:11610/1‑100 (MQ=255)                                                                                                              gctTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgcg                                                                                     <  5:11610/100‑1 (MQ=255)                                                                                                               ctTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCgg                                   <  7:236511/149‑1 (MQ=255)                                                                                                                ttAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGc                                  >  3:187514/1‑149 (MQ=255)                                                                                                                   gcggcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAac                                                                     <  6:167656/111‑1 (MQ=255)                                                                                                                   gcggcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAac                                                                     >  5:167656/1‑111 (MQ=255)                                                                                                                        gATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCaa                                                                                  >  3:233300/1‑93 (MQ=255)                                                                                                                        gATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCaa                                                                                  <  4:233300/93‑1 (MQ=255)                                                                                                                        gATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCg                          >  1:194321/1‑149 (MQ=255)                                                                                                                             tCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCctct                     <  8:113567/149‑1 (MQ=255)                                                                                                                             tCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCctct                     >  3:9601/1‑149 (MQ=255)                                                                                                                               gATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTg                   <  4:9601/149‑1 (MQ=255)                                                                                                                                aTGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGc                  <  3:262842/149‑1 (MQ=255)                                                                                                                                 tGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACAt                                      <  3:262404/128‑1 (MQ=255)                                                                                                                                 tGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAAAACTTCCAGGTCGCTGCACACCACCGGTACAt                                      >  4:262404/1‑128 (MQ=255)                                                                                                                                    atcgAGCGCTGTACGGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCgg                                   >  8:43736/4‑128 (MQ=255)                                                                                                                                    atcgAGCGCTGTACGGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCgg                                   <  7:43736/125‑1 (MQ=255)                                                                                                                                        aGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAg          >  4:214075/1‑149 (MQ=255)                                                                                                                                               taaaGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTAGGCCACCGGCAGGTAGCCACCCGGGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAgg   <  6:88894/146‑1 (MQ=255)                                                                                                                                                aCAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAGGc  >  6:242306/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCAGGC  >  NC_002947/5619327‑5619617

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