Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,619,466	C→T	8.8%	T345T (ACG→ACA)	PP_4938 ←	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,619,466	0	C	T	8.8%	192.0 / 6.9	68	T345T (ACG→ACA)	PP_4938	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base T (2/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

TTGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAG  >  NC_002947/5619323‑5619609                                                                                                                                                 |                                                                                                                                                ttGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTa                                                                                                                                              >  6:5874/1‑147 (MQ=255)     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)                                cTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGcc                                                                                                                                                 >  7:84029/1‑113 (MQ=255)                                cTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGcc                                                                                                                                                 <  8:84029/113‑1 (MQ=255)                                 tGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTca                                                                                                                                     <  8:173476/124‑2 (MQ=255)                                 tGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTca                                                                                                                                     >  7:173476/1‑123 (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                                                                <  3:284555/149‑1 (MQ=255)                                                                             ccTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGa                                                                <  8:245694/149‑1 (MQ=255)                                                                               tGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGaaa                                                              >  4:166655/1‑149 (MQ=255)                                                                               tGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGaaa                                                              >  1:23238/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                                                >  4:262842/1‑149 (MQ=255)                                                                                             cGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGc                                                >  8:123702/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                                                             >  5:167656/1‑111 (MQ=255)                                                                                                                       gcggcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAac                                                             <  6:167656/111‑1 (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                           <  7:43736/125‑1 (MQ=255)                                                                                                                                        atcgAGCGCTGTACGGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCgg                           >  8:43736/4‑128 (MQ=255)                                                                                                                                            aGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAg  >  4:214075/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                TTGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGG‑CGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGCCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAG  >  NC_002947/5619323‑5619609

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