Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,788,869	G→C	6.1%	V1672V (GTC→GTG)	pvdI ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,788,869	0	G	C	6.1%	239.4 / 5.1	76	V1672V (GTC→GTG)	pvdI	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (36/35);  new base C (3/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TCAGTGCGCCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAGGATCGCCA  >  NC_002947/4788721‑4789017                                                                                                                                                       |                                                                                                                                                     tCAGTGCGCCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑GCC‑TGCGTGa                                                                                                                                                     >  2:876691/1‑150 (MQ=255) tCAGTGCGCCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGa                                                                                                                                                     <  2:780773/150‑1 (MQ=255)    gTGCGCCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACgc                                                                                                                                                  >  1:653255/1‑150 (MQ=255)      gcgcCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACgcag                                                                                                                                                <  1:52657/150‑1 (MQ=255)       cgcCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTCACgcagg                                                                                                                                               >  2:361307/1‑150 (MQ=255)         cccTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACgcaggca                                                                                                                                             >  1:1084092/1‑150 (MQ=255)          ccTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTCACgcaggtgt                                                                                                                                            <  2:543987/150‑4 (MQ=255)                       cacaCCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGcaggcagg                                                                                                                                           <  1:80250/138‑1 (MQ=255)                       cacaCCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGcaggcagg                                                                                                                                           >  2:80250/1‑138 (MQ=255)                       cacaCCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAgctg                                                                                                                               >  2:632123/1‑150 (MQ=255)                              ttGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGctcactc                                                                                                                        >  2:255871/1‑150 (MQ=255)                                  gccgccCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACATGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGc                                                                                                                    >  2:73790/1‑150 (MQ=255)                                     gccCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAgctgct                                                                                                                             <  2:360760/138‑1 (MQ=255)                                     gccCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAgctgct                                                                                                                             >  1:360760/1‑138 (MQ=255)                                     gccCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTgg                                                                                                                 >  2:372476/1‑150 (MQ=255)                                          gTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTgg                                                                                                            <  2:123548/150‑1 (MQ=255)                                          gTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTgg                                                                                                            <  2:243327/150‑1 (MQ=255)                                              aGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGctcactca                                                                                                                       >  2:418690/1‑135 (MQ=255)                                              aGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGctcactca                                                                                                                       <  1:418690/135‑1 (MQ=255)                                                cccGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTgg                                                                                                            <  2:606794/144‑1 (MQ=255)                                                cccGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTgg                                                                                                            >  1:606794/1‑144 (MQ=255)                                                cccGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTcagc                                                                                                      >  2:4851/1‑150 (MQ=255)                                                 ccGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACATGTGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTcagca                                                                                                     <  1:73790/150‑1 (MQ=255)                                                  cGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTcagcag                                                                                                    >  1:159435/1‑150 (MQ=255)                                                   gAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCagcagc                                                                                                   <  1:546730/150‑1 (MQ=255)                                                      gtgtACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAag                                                                                                <  1:90560/150‑1 (MQ=255)                                                      gtgtACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCATGGTTCAGCAGCAag                                                                                                <  2:653255/150‑1 (MQ=255)                                                           cATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAaggccag                                                                                           >  1:596900/1‑150 (MQ=255)                                                           cATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAaggccag                                                                                           >  1:134619/1‑150 (MQ=255)                                                            aTGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACgcag                                                                                                                                                <  1:496910/96‑1 (MQ=255)                                                            aTGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACgcag                                                                                                                                                >  2:496910/1‑96 (MQ=255)                                                            aTGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAaggccagg                                                                                          >  2:763843/1‑150 (MQ=255)                                                            aTGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGAAGCGGCGCTGCCTTACATGA‑‑CAGACAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAaggccagg                                                                                          >  1:1177819/1‑150 (MQ=255)                                                             tGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTcagca                                                                                                     >  2:381593/1‑138 (MQ=255)                                                             tGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTcagca                                                                                                     <  1:381593/138‑1 (MQ=255)                                                                    ggCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTc                                                                                  >  1:287190/1‑150 (MQ=255)                                                                    ggCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTc                                                                                  <  2:287190/150‑1 (MQ=255)                                                                        aGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCg                                                                              >  2:507033/1‑150 (MQ=255)                                                                            tGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGGTCa                                                                          >  1:531606/1‑150 (MQ=255)                                                                              gCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATg                                                                        <  2:1084092/150‑1 (MQ=255)                                                                                 gggTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGAGCGTCGAGGAAC‑ACA‑TGCGTGACGCAGGCAAGCACCGCAAGCGACGCGCGCAGTCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAg                                                                     >  2:1041709/1‑150 (MQ=255)                                                                                   gTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAgctgc                                                                                                                              >  2:748511/1‑91 (MQ=255)                                                                                   gTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAgctgc                                                                                                                              <  1:748511/91‑1 (MQ=255)                                                                                   gTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATgcacgc                                                                   <  2:134619/150‑6 (MQ=255)                                                                                    tCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTCACgcaggtgt                                                                                                                                            <  1:870799/76‑4 (MQ=255)                                                                                    tCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTCACgcaggtgt                                                                                                                                            >  2:870799/1‑73 (MQ=255)                                                                                      gATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAg                                                                >  1:1050982/1‑150 (MQ=255)                                                                                            gACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTc                                                                               <  2:323863/129‑1 (MQ=255)                                                                                            gACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTc                                                                               >  1:323863/1‑129 (MQ=255)                                                                                               cTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGc                                                       >  1:514863/1‑150 (MQ=255)                                                                                                  cGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCAGCCgg                                                    <  1:268968/150‑1 (MQ=255)                                                                                                   gggTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCTATGCCGCTGTCCTCGATCATGTAGGCCa                                                                 <  1:220988/136‑1 (MQ=255)                                                                                                   gggTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCa                                                                 >  2:220988/1‑136 (MQ=255)                                                                                                          aaGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAg                                            <  1:372476/150‑1 (MQ=255)                                                                                                           aGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAgg                                           <  2:514863/150‑1 (MQ=255)                                                                                                            gCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAggg                                          <  1:507033/150‑1 (MQ=255)                                                                                                              gTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTc                                        >  1:842740/1‑150 (MQ=255)                                                                                                                  ccGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGAt                                                                            >  1:206263/1‑110 (MQ=255)                                                                                                                  ccGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGAt                                                                            <  2:206263/110‑1 (MQ=255)                                                                                                                    gTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCgg                                  <  2:764846/150‑1 (MQ=255)                                                                                                                       ttCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGcac                               >  1:680258/1‑150 (MQ=255)                                                                                                                           agcagcGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAg                           <  2:159435/150‑1 (MQ=255)                                                                                                                              agcGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTc                                                          >  2:465577/1‑116 (MQ=255)                                                                                                                              agcGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTc                                                          <  1:465577/116‑1 (MQ=255)                                                                                                                               gcgcGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATc                                                                           <  2:1098645/98‑1 (MQ=255)                                                                                                                               gcgcGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATc                                                                           >  1:1098645/1‑98 (MQ=255)                                                                                                                                  cgTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATACAGTACGCCGGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAgcgc                    <  1:876691/150‑1 (MQ=38)                                                                                                                                   gTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGcc                   <  2:842740/150‑1 (MQ=255)                                                                                                                                        ggCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCt              >  2:949287/1‑150 (MQ=255)                                                                                                                                           ac‑acC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGGTCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGa           <  2:531606/150‑1 (MQ=255)                                                                                                                                            c‑acC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAg          <  1:632123/150‑1 (MQ=255)                                                                                                                                            c‑acC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGGCCAGCGGCACATTGGCAGTGCCCACCGTGAg          <  2:920257/150‑1 (MQ=255)                                                                                                                                             tgcc‑TGCGTCACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAg          >  2:790578/3‑150 (MQ=255)                                                                                                                                                c‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAGGAt       >  1:199223/1‑150 (MQ=255)                                                                                                                                                     gTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAGGATCGcc   >  1:890232/1‑150 (MQ=255)                                                                                                                                                      tGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAGGATCGCCa  <  2:890232/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     TCAGTGCGCCCTGGGTGATCGCCACACCCTTGGGCCGCCCGGTGGAGCCCGAGGTGTACATGACGTAGGCCAGCGTGGCAGGGTCGATGGCGACCTGCGGGTTGGAAGCGTCACCGTCTTCGAGCAGCGCGTCGAGGCAC‑ACC‑TGCGTGACGCAGGCAGGCACCGCCAGCTGCTCACTCAGGCTGGCCTGGGTCAGCAGCAAGGCCAGGCCGCTGTCCTCGATCATGTAGGCCAGGCGTTCGGCCGGGTAGTCAGGGTCCAGCGGCACATAGGCAGCGCCCACCTTGAGGATCGCCA  >  NC_002947/4788721‑4789017

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑