Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,941,126	Δ1 bp	8.1%	coding (118/531 nt)	PP_3475 ←	type II secretion system pseudopilin PulG‑like signal peptide protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,941,126	0	T	.	8.1%	173.4 / 10.8	62	coding (118/531 nt)	PP_3475	type II secretion system pseudopilin PulG‑like signal peptide protein
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base . (3/2);  total (33/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 = 2.13e-01

TGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGT  >  NC_002947/3940978‑3941273                                                                                                                                                      |                                                                                                                                                    tGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATCCTCTGTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGGG‑GCGCTc                                                                                                                                                    <  1:308239/150‑1 (MQ=255) tGCGGGTAGCGGTCGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTc                                                                                                                                                    >  2:208932/1‑150 (MQ=255)        aGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGa                                                                                                                                             <  1:105697/150‑1 (MQ=255)           ggTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCg                                                                                                                                          >  2:555400/1‑150 (MQ=255)                tCCTGTAGCAGGTCGGCCAGATCCTGGGGGTGCTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCg                                                                                                                                     >  2:745796/1‑150 (MQ=255)                 ccTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCgg                                                                                                                                    >  1:395970/1‑150 (MQ=255)                  cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCggg                                                                                                                                   >  1:382439/1‑149 (MQ=255)                  cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGt                                                                                                                                   >  2:306854/1‑150 (MQ=255)                  cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGt                                                                                                                                   <  1:97521/150‑1 (MQ=255)                   tGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTg                                                                                                                                  <  1:208932/150‑1 (MQ=255)                    gTAGCAGGTCGGCCAGATGCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTgg                                                                                                                                 >  1:37843/1‑150 (MQ=255)                         aGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGgccgcc                                                                                                                            >  1:114465/1‑150 (MQ=255)                               gCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGcc                                                                                                                      >  2:161657/1‑150 (MQ=255)                                ccAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCa                                                                                                                     >  2:427286/1‑150 (MQ=255)                                  aGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCCAGGCTTGGGGATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGCTGGCCGCCCGGGacaca                                                                                                                   >  2:171537/1‑149 (MQ=255)                                  aGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACg                                                                                                                   >  2:39256/1‑150 (MQ=255)                                   gATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGC‑CGCCCAg                                                                                                                                              >  1:832496/1‑122 (MQ=255)                                   gATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGC‑CGCCCAg                                                                                                                                              <  2:832496/122‑1 (MQ=255)                                      ccTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACgctgc                                                                                                               <  1:427286/150‑1 (MQ=255)                                       cTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAg                                                                                                                          <  2:726630/138‑1 (MQ=255)                                       cTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAg                                                                                                                          >  1:726630/1‑138 (MQ=255)                                       cTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACgctgct                                                                                                              >  2:225435/1‑150 (MQ=255)                                           gggTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTTGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGAcc                                                                                                          <  1:977853/150‑1 (MQ=255)                                             gTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGAACGACTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGc                                                                                                        >  2:588093/1‑150 (MQ=255)                                              cccTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCa                                                                                                       >  2:671936/3‑150 (MQ=255)                                                cTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAAt                                                                                                     <  2:395970/150‑1 (MQ=255)                                                    gCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagc                                                                                                 <  1:654659/150‑1 (MQ=255)                                                      cAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                               <  1:555400/150‑1 (MQ=255)                                                            tGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCAc                                                                                         >  2:151113/1‑150 (MQ=255)                                                                  ggCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAg                                                                                   <  2:359592/150‑1 (MQ=255)                                                                   gcgcGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAgg                                                                                  <  1:164592/150‑1 (MQ=255)                                                                             agCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCAAGGTATGTAAAgc                                                                        >  1:768114/1‑150 (MQ=255)                                                                                    cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTAt                                                                               <  2:765906/136‑1 (MQ=255)                                                                                    cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTAt                                                                               >  1:765906/1‑136 (MQ=255)                                                                                     aaGGCTTGTGCATCCTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTg                                                                <  1:500423/150‑1 (MQ=255)                                                                                      aGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTgg                                                               <  1:225435/150‑1 (MQ=255)                                                                                            gggCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                         <  1:671936/150‑1 (MQ=255)                                                                                              gCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCa                                                       >  1:584111/1‑150 (MQ=255)                                                                                                    tGGCTGTCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGc                                                 >  2:481504/1‑150 (MQ=255)                                                                                                    tGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTTGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAg                                                                         >  1:630147/1‑126 (MQ=255)                                                                                                    tGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAg                                                                         <  2:630147/126‑1 (MQ=255)                                                                                                        tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                             <  1:306854/150‑1 (MQ=255)                                                                                                        tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                             >  2:818179/1‑150 (MQ=255)                                                                                                            aaCCCACAGCAACTGGCGTCCGCGGCCGCGCTGGGCG‑GCCCCCGCCCAGATCGGGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCGGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGc                                         <  1:1088860/149‑1 (MQ=255)                                                                                                             aCCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                        >  1:498958/1‑150 (MQ=255)                                                                                                                 acaGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGAAGTTCGGCGGACCATTCAGCGAAGAGGCTGCCa                                    >  2:529291/1‑150 (MQ=255)                                                                                                                       aCTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCAGTCAGCGAAGAGGCTGCCATCAcgc                              <  1:151113/150‑1 (MQ=255)                                                                                                                          ggCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGcc                           <  1:529291/150‑1 (MQ=255)                                                                                                                          ggCGTGCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGGCGCGATAGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGcc                           <  1:243503/150‑1 (MQ=255)                                                                                                                             gTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACgct                                                                                                                 >  2:54934/1‑61 (MQ=255)                                                                                                                             gTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACgct                                                                                                                 <  1:54934/61‑1 (MQ=255)                                                                                                                             gTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGccctc                        >  2:96722/1‑150 (MQ=255)                                                                                                                                    gTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTgg                 <  1:393562/150‑1 (MQ=255)                                                                                                                                    gTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGAAGCAGCACCCCCAGGTATGTAAAGCCACCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTgg                 <  2:37843/150‑1 (MQ=255)                                                                                                                                       gcgcTGGGCG‑GCACTCGCCCAGGGCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCAc              <  1:324711/150‑1 (MQ=255)                                                                                                                                       gatcTGGGCGAGCGC‑CGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCa                                    >  2:555507/4‑128 (MQ=255)                                                                                                                                       gatcTGGGCGAGCGC‑CGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCa                                    <  1:555507/125‑1 (MQ=255)                                                                                                                                         tcTGGGCGAGCGC‑CGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGCTGGGGCTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCg            >  1:618694/2‑150 (MQ=255)                                                                                                                                         gcTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGt                                                                             >  2:601490/1‑85 (MQ=255)                                                                                                                                         gcTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGt                                                                             <  1:601490/85‑1 (MQ=255)                                                                                                                                         gcTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCg            >  2:455189/1‑150 (MQ=255)                                                                                                                                                    gcTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGt  >  1:659216/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCG‑GCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGT  >  NC_002947/3940978‑3941273

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