Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,989,571	G→C	7.3%	G274A (GGT→GCT)	PP_1780 →	mannosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,989,571	0	G	C	7.3%	207.0 / 3.9	65	G274A (GGT→GCT)	PP_1780	mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base C (2/3);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

TTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCTGATGAGCGTATGG  >  NC_002947/1989423‑1989717                                                                                                                                                     |                                                                                                                                                     ttCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTggg                                                                                                                                                      >  6:70996/1‑149 (MQ=255)    gCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGCAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGt                                                                                                                                                     >  6:38649/1‑147 (MQ=255)     ccAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAg                                                                                                                                                  >  4:320909/1‑149 (MQ=255)     ccAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCAAGCAACTTCAAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATACGGCATTGTAGCAAGTTAACCTGGAGCTGGCTCAg                                                                                                                                                  >  6:11253/1‑149 (MQ=255)           ccTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTGACCCAGCTCCAg                                                                                                                                            <  7:267210/149‑1 (MQ=255)            cTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTGACCCAGCTCCAgg                                                                                                                                           <  5:133389/149‑1 (MQ=255)            cTTCGATGCGGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGTGCTGACCCAGCTCCAgg                                                                                                                                           <  5:366297/149‑1 (MQ=255)             tgcGATGATGAGTATGCGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGACGAAAGCTGCCATTAGGCATTGTCGCAAGTTAACCTGGACCTGGGTCAGCACCAgg                                                                                                                                           >  6:210825/3‑148 (MQ=255)                  tGCTGAGTATGGGGGGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGccc                                                                                                                                     <  7:49866/149‑1 (MQ=255)                      gAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTtgt                                                                                                                                 >  5:262432/1‑149 (MQ=255)                       aGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGCAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTtgt                                                                                                                                 <  5:38649/148‑1 (MQ=255)                       aGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTtgt                                                                                                                                 <  5:70763/148‑1 (MQ=255)                        gTAGGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTgg                                                                                                                               <  2:122346/149‑1 (MQ=255)                         tATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGa                                                                                                                              >  5:249972/1‑149 (MQ=255)                         tATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGa                                                                                                                              >  6:221944/1‑149 (MQ=255)                          aTGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGC‑AGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGAtt                                                                                                                            >  8:202517/1‑149 (MQ=255)                             gggTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGAttt                                                                                                                           >  6:349834/1‑148 (MQ=255)                              ggTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGa                                                                                                                         <  6:249972/149‑1 (MQ=255)                              ggTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGa                                                                                                                         <  3:320909/149‑1 (MQ=255)                               gTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAg                                                                                                                        >  2:248254/1‑149 (MQ=255)                                    ccTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGTCCTTGTGGATTTGAG‑CCAg                                                                                                                   >  6:118169/1‑148 (MQ=255)                                        gAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGaa                                                                                                                <  5:221944/147‑1 (MQ=255)                                          ggACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTc                                                                                                             <  5:77363/148‑1 (MQ=255)                                           gACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAg                                                                                                                                            <  3:371140/117‑1 (MQ=255)                                           gACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAg                                                                                                                                            >  4:371140/1‑117 (MQ=255)                                                   tGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATc                                                                                                   >  4:100018/1‑149 (MQ=255)                                                     aGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATcgc                                                                                                 >  1:89417/1‑149 (MQ=255)                                                      gcgcTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTAGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTAGTGGATTTAAC‑CCAGGAAATCGATCGtactagct                                                                                                >  8:121139/1‑142 (MQ=255)                                                           aTGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGa                                                                                                           <  7:11053/133‑1 (MQ=255)                                                           aTGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGa                                                                                                           >  8:11053/1‑133 (MQ=255)                                                            tGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTAGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGc                                                                                                         <  7:25475/134‑1 (MQ=255)                                                            tGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTAGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGc                                                                                                         >  8:25475/1‑134 (MQ=255)                                                              aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACt                                                                                         <  5:70996/148‑1 (MQ=255)                                                              aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                        <  6:262432/149‑1 (MQ=255)                                                              aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                        >  6:174829/1‑149 (MQ=255)                                                              aagaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                        >  4:185727/1‑149 (MQ=255)                                                               agaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTgg                                                                                       <  3:369579/149‑1 (MQ=255)                                                               agaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                        >  2:197409/1‑148 (MQ=255)                                                                gaaAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTg                                                                                        >  2:293382/1‑147 (MQ=255)                                                                  aaGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGgaga                                                                                    <  5:349834/149‑1 (MQ=255)                                                                     gctggctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGcc                                                                                 <  6:9056/149‑1 (MQ=255)                                                                     gctggctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGcc                                                                                 <  6:14287/149‑1 (MQ=255)                                                                     gctggctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑‑CAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCt                                                                                >  1:203251/1‑149 (MQ=255)                                                                         gctAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTcaa                                                                             <  7:301726/149‑1 (MQ=255)                                                                               aaCTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAgt                                                                       <  3:100018/149‑1 (MQ=255)                                                                                agttCCAGCAAGTGGGCGACAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGGTGGGTCCGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAgtg                                                                      <  5:210825/147‑1 (MQ=255)                                                                                 cTTCCATCAAGTTGTCGAAAGCGTCGAATGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCTGGCAGGCGGTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGATCCTCAAGCAAgtgg                                                                     <  5:11253/149‑1 (MQ=255)                                                                                    ccAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAgtggt                                                                    <  5:120504/147‑1 (MQ=255)                                                                                          aGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAGCCCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGc                                                             <  1:293382/149‑1 (MQ=255)                                                                                               gCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAac                                                       <  2:163297/149‑1 (MQ=255)                                                                                                 gAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAacca                                                     <  1:242050/149‑1 (MQ=255)                                                                                                 gAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAacca                                                     <  3:233850/149‑1 (MQ=255)                                                                                                   aaGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAaccacc                                                   <  1:248254/149‑1 (MQ=255)                                                                                                        tCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTg                                              >  4:365645/1‑149 (MQ=255)                                                                                                              gCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGAt                                        <  7:216810/149‑1 (MQ=255)                                                                                                               cTGACATTCTGCATTCTCGAAAGTTAACCTTGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATa                                       <  7:121139/149‑1 (MQ=255)                                                                                                                tGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATAc                                      >  7:313026/1‑149 (MQ=255)                                                                                                                   cATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑‑CAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑g                                 <  2:203251/149‑1 (MQ=255)                                                                                                                    aTTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑g                                 <  5:174829/149‑1 (MQ=255)                                                                                                                         gCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACg                            >  1:169567/1‑149 (MQ=255)                                                                                                                         gCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACg                            >  7:21103/1‑149 (MQ=255)                                                                                                                                 gCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGctc                    >  2:323737/1‑149 (MQ=255)                                                                                                                                     gTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGc                >  2:48481/1‑149 (MQ=255)                                                                                                                                        tgccTGGTGCTGACCCAGCTCCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCtga             >  7:199847/3‑149 (MQ=255)                                                                                                                                           cTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CAAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCCCCATCTGTCCGATACCCA‑GCAACGGATGGCTCTTGCtgat            >  6:253865/1‑147 (MQ=255)                                                                                                                                              gAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCACGCAACGGCTGGCTCTTGCTGATGAGc        <  1:197409/149‑1 (MQ=255)                                                                                                                                                    ggTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCTGATGAGCGTATgg  <  3:185727/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TTCGCCAGCGCCTTCGATGCTGAGTATGGGGTGACCCTTGAGGACCTGGCTGAGCGCTATGAAGAAAGGCTGGCTAGCAACTTCCAGCAAGTGGGCGAAAGCGTCGAAAGCTGCCATTCGGCATTGTCGCAAGTTAACCTGGAGCTGGGTCAGCACCAGGCAGCCCTTGTGGATTTGAG‑CCAGGAAGTCGAGCGTCATCGCGCCGAACTGGAGAGCCTCAAGCAAGTGGTGAACGGCCTGAACCACCATCTGTCCGATACCCA‑GCAACGGCTGGCTCTTGCTGATGAGCGTATGG  >  NC_002947/1989423‑1989717

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