Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,767,819	A→G	9.1%	intergenic (+33/+36)	fpvA → / ← PP_4218	TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,767,819	0	A	G	9.1%	183.6 / 6.7	66	intergenic (+33/+36)	fpvA/PP_4218	TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (2/4);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAG  >  NC_002947/4767674‑4767965                                                                                                                                                   |                                                                                                                                                   acGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGt                                                                                                                                                 <  1:175526/149‑1 (MQ=255)     aTTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCgg                                                                                                                                             >  5:369805/1‑149 (MQ=255)      ttCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGt                                                                                                                                            >  6:145277/1‑149 (MQ=255)        cGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtt                                                                                                                                           <  5:145277/148‑1 (MQ=255)         gACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGc                                                                                                                                               >  2:89506/1‑143 (MQ=255)         gACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGc                                                                                                                                               <  1:89506/143‑1 (MQ=255)           cAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACTCAAAAAACCGCACCTGGGTGCGGttt                                                                                                                                          <  2:357800/146‑1 (MQ=255)           cAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGtttttt                                                                                                                                       <  6:240813/149‑1 (MQ=255)           cAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGttt                                                                                                                                          >  1:357800/1‑146 (MQ=255)             agaagTACTACCCGGGCATTGCCAGCAACTACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCGGGGTGCGGttttttgc                                                                                                                                     <  1:56382/149‑3 (MQ=255)             agaagTACTACCATGGCATTGCCAGCAACAACTCCTTGAGCCCGT‑ACGACGTGTATGTTGACCCACGCAACTTCACCATCACCTCCAAGTACAGCTTCTGATCTCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                     <  7:283965/149‑1 (MQ=255)               aagTACTACCAGGTCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGttttttgcgt                                                                                                                                   <  1:177590/149‑5 (MQ=255)                 gTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGa                                                                                                                                 >  7:218406/1‑149 (MQ=255)                 gTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGa                                                                                                                                 >  8:326054/1‑149 (MQ=255)                  tactacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGAt                                                                                                                                <  7:326054/149‑1 (MQ=255)                  tactacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGAt                                                                                                                                >  2:343212/1‑149 (MQ=255)                  tactacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGAt                                                                                                                                <  1:120957/149‑1 (MQ=255)                   actacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGAt                                                                                                                                <  2:108913/148‑1 (MQ=255)                   actacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATg                                                                                                                               >  2:179547/1‑149 (MQ=255)                    ctacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGc                                                                                                                              <  1:343212/149‑1 (MQ=255)                      acCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTGACGACGTTTATGGTGACCAACGAAACTTCACCATCACCGCCAAGTACAGGTTGTGATCGCATGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGca                                                                                                                             >  1:67918/1‑148 (MQ=255)                       ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                     >  7:271381/1‑139 (MQ=255)                       ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                     <  8:271381/139‑1 (MQ=255)                       ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCt                                                                                                                            >  7:278008/1‑148 (MQ=255)                        cAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTg                                                                                                                                                >  1:366624/1‑127 (MQ=255)                        cAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTg                                                                                                                                                <  2:366624/127‑1 (MQ=255)                         aGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGa                                                                                                                          >  5:31758/1‑148 (MQ=255)                              aTTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCt                                                                                                                    >  5:197172/1‑149 (MQ=255)                               ttGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                                   <  8:218406/149‑1 (MQ=255)                                 gCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTcc                                                                                                                  <  7:69357/148‑1 (MQ=255)                                 gCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTcc                                                                                                                  <  1:179547/148‑1 (MQ=255)                                  ccAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCt                                                                                                                >  2:325389/1‑149 (MQ=255)                                     gCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTa                                                                                                              >  7:19348/1‑148 (MQ=255)                                     gCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTAc                                                                                                             <  7:294736/149‑1 (MQ=255)                                            cTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgct                                                                                                      >  1:361666/1‑149 (MQ=255)                                               cTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgct                                                                                                   >  7:147240/1‑149 (MQ=255)                                               cTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgc                                                                                                    <  8:19348/148‑1 (MQ=255)                                               cTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctg                                                                                                     <  2:361666/147‑1 (MQ=255)                                                  gAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                                  <  4:121532/147‑1 (MQ=255)                                                                gtgtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTAc                                                                                                             <  8:207611/123‑1 (MQ=255)                                                                gtgtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTAc                                                                                                             >  7:207611/1‑123 (MQ=255)                                                                        tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGa                                                                                                          >  5:222668/1‑118 (MQ=255)                                                                        tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGa                                                                                                          <  6:222668/118‑1 (MQ=255)                                                                                gCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAgg                                                                                                                      >  4:303035/1‑98 (MQ=255)                                                                                gCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAgg                                                                                                                      <  3:303035/98‑1 (MQ=255)                                                                                     ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttt                                                                                                                                         <  1:189386/74‑1 (MQ=255)                                                                                     ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttt                                                                                                                                         >  2:189386/1‑74 (MQ=255)                                                                                     ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGc                                                              >  8:163115/1‑149 (MQ=255)                                                                                         caaacACCGCCAAGTACAGCTTCTGATCGACTGACGTTGAACGCAAAAAACCGCACCCACGTGCGGATTTTTCATCGATGCCAGAAAGCCTCCCTTACAGATGCTGCTGCAGTCCCTCGACCACCATCGCCAGGCCTTCAACGGactcc                                                          >  2:190991/5‑149 (MQ=255)                                                                                          cATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcct                                                         >  5:129888/1‑149 (MQ=255)                                                                                                 gCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgatgat                                                  <  8:147240/149‑1 (MQ=255)                                                                                                  ccAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                                 <  6:129888/149‑1 (MQ=255)                                                                                                  ccAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                                 <  8:278008/149‑1 (MQ=255)                                                                                                       tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAcc                                            >  2:369739/1‑149 (MQ=255)                                                                                                         cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAccc                                           >  5:236699/1‑148 (MQ=255)                                                                                                         cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCa                                          <  2:67918/149‑1 (MQ=255)                                                                                                                   tCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATg                                  >  1:56621/1‑147 (MQ=255)                                                                                                                        tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                           >  1:20173/1‑149 (MQ=255)                                                                                                                                   tgaaaaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTccc                >  5:299811/3‑149 (MQ=255)                                                                                                                                      aaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg             <  6:236699/149‑1 (MQ=255)                                                                                                                                        aaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGAt           <  7:163115/149‑1 (MQ=255)                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCa                                                                   >  7:349434/1‑84 (MQ=255)                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCa                                                                   <  8:349434/84‑1 (MQ=255)                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATggg                        >  5:37547/1‑127 (MQ=255)                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATggg                        <  6:37547/127‑1 (MQ=255)                                                                                                                                                 ccAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAg  >  3:405925/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGT‑ACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAG  >  NC_002947/4767674‑4767965

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