Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,273,064	G→A	9.5%	intergenic (+40/+21)	PP_3743 → / ← glcC	hypothetical protein/transcriptional regulator GlcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,273,064	0	G	A	9.5%	174.9 / 7.5	63	intergenic (+40/+21)	PP_3743/glcC	hypothetical protein/transcriptional regulator GlcC
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base A (3/3);  total (34/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 = 8.07e-01

GTCGGCGCTTGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCTGCTCGGGCAG  >  NC_002947/4272919‑4273210                                                                                                                                                  |                                                                                                                                                     gTCGGCGCTTGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTATCGGCCCTCTCGCGTCGAACCGCCGCGAAAgg                                                                                                                                                   <  8:314394/149‑1 (MQ=38) gTCGGCGCTTGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTATCGGCCCTCTCGCGTCGAACCGCCGCGAAAgg                                                                                                                                                   <  4:159314/149‑1 (MQ=38)    ggCGCTTGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGCCCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGcc                                                                                                                                                >  1:80982/1‑149 (MQ=255)       gcTTGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCGTTCGCGGCGGTGCGCCGCGAGCGGGcc                                                                                                                                                >  6:199672/1‑146 (MQ=255)         ttGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCgaga                                                                                                                                                     >  6:309551/1‑139 (MQ=255)         ttGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCgaga                                                                                                                                                     <  5:309551/139‑1 (MQ=255)            aaCAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACaa                                                                                                                                         <  2:249961/148‑1 (MQ=255)             aCAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAgg                                                                                                                                       <  8:68575/149‑1 (MQ=255)                   ccGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCtt                                                                                                                                 >  6:156480/1‑149 (MQ=255)                      cAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCAc                                                                                                                              <  1:339422/149‑1 (MQ=255)                          gACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCAccccc                                                                                                                          >  8:42916/1‑149 (MQ=255)                          gACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCAccccc                                                                                                                          <  5:156480/149‑1 (MQ=255)                                  tGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGccatc                                                                                                                   >  1:345570/1‑148 (MQ=255)                                     ggtGCTGGCCAACGACGAAGAGCTGGCCCGGGGCATCCTGTCGCAAGATGACCCTGTGGATGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATa                                                                                                               <  8:277598/147‑1 (MQ=255)                                         cAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCa                                                                                                           >  3:82432/1‑149 (MQ=255)                                           ggCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCa                                                                                                           <  7:184137/147‑1 (MQ=255)                                           ggCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCa                                                                                                           >  8:184137/1‑147 (MQ=255)                                              cAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCa                                                                                                                               >  6:348374/1‑124 (MQ=255)                                              cAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTc                                                                                                    >  5:47759/1‑149 (MQ=255)                                               aaCGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCg                                                                                                   <  4:82432/149‑1 (MQ=255)                                                  gacgaAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGccc                                                                                                <  8:237706/149‑1 (MQ=255)                                                   acgaAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCt                                                                                               <  6:106448/149‑1 (MQ=255)                                                   acgaAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCt                                                                                               <  4:188344/149‑1 (MQ=255)                                                     gaAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGa                                                                                             >  7:208061/1‑149 (MQ=255)                                                        gagCTGGACAGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAACCCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCa                                                                                          <  6:271044/149‑1 (MQ=255)                                                              gCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGCCGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGgcgc                                                                                     <  5:199672/148‑1 (MQ=255)                                                              gCCCGGGCCATCCTGTCGCAAGATGACCCTGGGGCTGACCATTTATTAAACCTGTTCCGGCCCTTTCGCGGGGGTTCTACGCTAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGgcgc                                                                                     <  2:313063/148‑1 (MQ=255)                                                                ccGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCgctgc                                                                                  >  8:205043/1‑149 (MQ=255)                                                                 cGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCgctgct                                                                                 >  5:353588/1‑149 (MQ=255)                                                                       aTCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCt                                                                           >  2:102677/1‑149 (MQ=255)                                                                               gCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGAtc                                                                   >  2:247153/1‑149 (MQ=255)                                                                                 aaGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGAtctc                                                                 >  5:338403/1‑149 (MQ=255)                                                                                               gCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGcc                                                                                                 <  6:288247/103‑1 (MQ=255)                                                                                               gCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGcc                                                                                                 >  5:288247/1‑103 (MQ=255)                                                                                               gCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTcgc                                                   >  7:318450/1‑149 (MQ=255)                                                                                                     cATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGc                                             >  3:50214/1‑149 (MQ=255)                                                                                                      aTTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTAGACGCGAGAGGGCCGATACAAGGCAGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCTCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTAGCTTAGGcc                                            >  2:262756/1‑149 (MQ=255)                                                                                                       ttGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTcg                                                    <  1:97696/140‑1 (MQ=255)                                                                                                       ttGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTcg                                                    >  2:97696/1‑140 (MQ=255)                                                                                                         gAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGt                                          <  1:102677/148‑1 (MQ=255)                                                                                                          aGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGCCCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGtt                                         >  5:80120/1‑148 (MQ=255)                                                                                                          aGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAACCTGTCGCGTATGCCGTTg                                        >  5:331237/1‑149 (MQ=255)                                                                                                          aGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGtt                                         <  6:338403/148‑1 (MQ=255)                                                                                                           gTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGa                                       <  8:182579/149‑1 (MQ=255)                                                                                                            tAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGa                                       >  2:47527/1‑148 (MQ=255)                                                                                                                cttGTATCGGCCCTCTCGCGTCGAACCGCCGCGAAAGGGCCGGTACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGc                                  >  3:43595/3‑149 (MQ=37)                                                                                                                ccTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGctc                                                                                >  2:237593/1‑103 (MQ=255)                                                                                                                ccTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGctc                                                                                <  1:237593/103‑1 (MQ=255)                                                                                                                 ttGTATCGGCCCTCTCGCGTCGAACCGCCGCGAAAGGGCCGGTACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTcgc                                                   >  5:315534/2‑131 (MQ=25)                                                                                                                 ttGTATCGGCCCTCTCGCGTCGAACCGCCGCGAAAGGGCCGGTACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTcgc                                                   <  6:315534/130‑1 (MQ=25)                                                                                                                 cTGTACCGGCCCTTTCGCGGCGGTTCGCCGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCt                                 >  2:236160/1‑149 (MQ=255)                                                                                                                    tatcGGCCCTCTCGCGTCGAACCGCCGCGAAAGGGCCGGTACAGGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCAc                              >  3:22650/4‑149 (MQ=37)                                                                                                                           ccTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGcgccgc                       >  8:222532/1‑149 (MQ=255)                                                                                                                           ccTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGAATG‑‑GTCGCGCTGACGAGGCGCTGCTCTTCCTGCACGATCTCCCGGAAACTGTCGCGTATGCCGTTGTGATGCTGACGcgccgc                       >  8:222535/1‑149 (MQ=255)                                                                                                                             tttCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATGTCGGCGCCATGACCAGGCGCTGCTCTTCCTGATCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGcgccgc                       >  1:280120/1‑149 (MQ=255)                                                                                                                               tCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGccgcccg                    <  6:47759/148‑1 (MQ=255)                                                                                                                                 gcggcgGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGtt                                         >  7:182287/1‑125 (MQ=255)                                                                                                                                 gcggcgGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGtt                                         <  8:182287/125‑1 (MQ=255)                                                                                                                                   ggcggTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGtt                                         <  3:329379/123‑1 (MQ=255)                                                                                                                                   ggcggTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGtt                                         >  4:329379/1‑123 (MQ=255)                                                                                                                                      ggTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCt            <  1:247153/149‑1 (MQ=255)                                                                                                                                           gACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTc                                                                             >  3:338854/1‑79 (MQ=255)                                                                                                                                           gACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTc                                                                             <  4:338854/79‑1 (MQ=255)                                                                                                                                           gACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCTGCTCg       <  6:353588/149‑1 (MQ=255)                                                                                                                                            aCGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCTGCTCgg      <  7:42916/149‑1 (MQ=255)                                                                                                                                                gagaGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGcgccgc                       >  7:317977/1‑128 (MQ=255)                                                                                                                                                gagaGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGcgccgc                       <  8:317977/128‑1 (MQ=255)                                                                                                                                                gagaGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCTGCTCGGGCAg  <  1:47527/149‑1 (MQ=255)                                                                                                                                                gagaGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCTGCTCGGGCAg  <  1:236160/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTCGGCGCTTGAACAGGCCCGCAAGGACCATCCTGGGCTGCAGGCCAACGACGAAGAGCTGGCCCGGGCCATCCTGTCGCAAGATGACCCTGTGGCTGACCATTGAGTAAACCTGTACCGGCCCTTTCGCGGCGGTTCGACGCGAGAGGGCCGATACAAGGCCGCTTCACCCCCAGCCATCCATACGCATG‑‑GTCGCCCTGACCAGGCGCTGCTCTTCCTGCTCGATCTCCCGCAAACTGTCGCGTATGCCGTTGATATGCTCACGCGCCGCCCGCTGCGCCTGCTCGGGCAG  >  NC_002947/4272919‑4273210

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