Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,545,948	G→A	8.5%	R350H (CGC→CAC)	PP_2238 →	metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,545,948	0	G	A	8.5%	165.7 / 3.4	58	R350H (CGC→CAC)	PP_2238	metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (23/30);  new base A (3/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACTGA  >  NC_002947/2545803‑2546094                                                                                                                                                  |                                                                                                                                                     gCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCa                                                                                                                                                   <  8:155953/149‑1 (MQ=255)     gggCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTc                                                                                                                                               <  5:328339/149‑1 (MQ=255)          cgGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGa                                                                                                                                          <  7:131248/149‑1 (MQ=255)            gctgctGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGc                                                                                                                                        <  3:224505/149‑1 (MQ=255)             ctgctgGTGGACCCGACCCGGGTCACCTGGGGCCTGCTGGACAACCTCGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑a                                                                                                                                     <  4:170607/149‑1 (MQ=255)                     ggACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGgcg                                                                                                                             >  2:103558/1‑149 (MQ=255)                         cggacccgGGTCACCTTCGGCCTGCGGGTCAACCTGGCGGCAGAGGTGGTTGTTGTCGATGGGCTAAACCCCTCGACGCTGAGCAAGTCCTGCAAGGGTGACGACGACGTGGTGCATATTGGCCATGTCATGGAGC‑‑ATGACGGCGCAg                                                                                                                          <  4:115212/146‑1 (MQ=255)                         ccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGc                                                                                                                         >  2:128065/1‑149 (MQ=255)                         ccgacccgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGc                                                                                                                         <  3:170055/149‑1 (MQ=255)                               cgGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTtgtg                                                                                                                   <  2:62706/149‑1 (MQ=255)                                gggTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTAATGGAGC‑‑AGGACGGCGCAGCCAtgtcc                                                                                                                  >  8:176357/1‑147 (MQ=255)                                  gTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGa                                                                                                                <  3:245900/149‑1 (MQ=255)                                       cTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAg                                                                                                                                                  <  1:134545/112‑1 (MQ=255)                                       cTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAg                                                                                                                                                  >  2:134545/1‑112 (MQ=255)                                       cTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGtttt                                                                                                           >  8:33125/1‑149 (MQ=255)                                        tGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGttttt                                                                                                          >  4:29755/1‑149 (MQ=255)                                          cGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTc                                                                                                         <  1:14313/148‑1 (MQ=255)                                          cGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCg                                                                                                        >  7:141670/1‑149 (MQ=255)                                                  tgGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTc                                                                                                >  4:189875/1‑149 (MQ=255)                                                    gACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTc                                                                                                >  5:301745/1‑147 (MQ=255)                                                    gACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCg                                                                                               <  8:141670/148‑1 (MQ=255)                                                         ccTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCg                                                                                          >  1:174429/1‑148 (MQ=255)                                                          cTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGaa                                                                                        >  4:302739/1‑149 (MQ=255)                                                              cggcAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTc                                                                                    <  1:265245/149‑1 (MQ=255)                                                                gcAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTCGAACCAGGCGaaaaactc                                                                                  <  7:33125/149‑7 (MQ=255)                                                                 cAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTCGAACCAGGCGaaaaactcg                                                                                 <  7:176357/149‑8 (MQ=255)                                                                     ggtggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCgg                                                                              >  1:130492/1‑148 (MQ=255)                                                                      gtggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGa                                                                            <  5:58323/149‑1 (MQ=255)                                                                       tggtTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAg                                                                           <  6:239632/149‑1 (MQ=255)                                                                           tCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGcc                                                                                                                        <  7:149023/100‑1 (MQ=255)                                                                           tCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGcc                                                                                                                        >  8:149023/1‑100 (MQ=255)                                                                           tCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGCTGACGACGAACTGGTGCATAATCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGAGAACCTCGGCCGGGAGGtca                                                                       >  5:270008/1‑149 (MQ=255)                                                                               gTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGtcatca                                                                    <  3:189875/148‑1 (MQ=255)                                                                               gTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCAc                                                                   >  4:291047/1‑149 (MQ=255)                                                                                   aGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAg                                                               >  6:277550/1‑149 (MQ=255)                                                                                    ggggCTAAACCCCAACAAGCTGATCAATTCCTTCAATGGGGTAGACGACCTTGTGGATATTCGCCATGTGATTGATCCTAAGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGa                                                                <  1:103558/149‑1 (MQ=255)                                                                                        cTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGAc                                                          <  3:29755/149‑1 (MQ=255)                                                                                              cccACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGa                                                    <  6:301745/149‑1 (MQ=255)                                                                                                 accacGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGa                                                 <  7:78684/149‑1 (MQ=255)                                                                                                      gCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCAAATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAAt                                            >  8:83250/1‑149 (MQ=255)                                                                                                          aGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGa                                         <  7:32504/148‑1 (MQ=255)                                                                                                                tCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCcgc                                  >  8:295214/1‑149 (MQ=255)                                                                                                                  cTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCcgcgc                                <  3:302739/149‑1 (MQ=255)                                                                                                                    gCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGt                              >  7:68077/1‑149 (MQ=255)                                                                                                                           tGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGATGAGCAATTGAGTGCCGCGCGTGCTc                          <  1:286203/146‑1 (MQ=255)                                                                                                                           tGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGATGAGCAATTGAGTGCCGCGCGTGCTc                          >  2:286203/1‑146 (MQ=255)                                                                                                                              cgacgacCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                    <  3:291047/149‑1 (MQ=255)                                                                                                                               catgacCTGGCGAATATGCACCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                    <  5:164367/145‑1 (MQ=255)                                                                                                                               catgacCTGGCGAATATGCACCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGc                    >  6:164367/4‑148 (MQ=255)                                                                                                                                 tgacCTGGCGAATATGCACCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCg                  >  1:173360/2‑148 (MQ=255)                                                                                                                                    ccTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAAtt              <  7:295214/149‑1 (MQ=255)                                                                                                                                    ccTGGTGCATATTCACCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAAtt              >  3:313682/1‑149 (MQ=255)                                                                                                                                    ccTGGTGCATATTCACCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAAtt              <  4:313682/149‑1 (MQ=255)                                                                                                                                     cTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAAttt             >  1:195777/1‑149 (MQ=255)                                                                                                                                        gTGCAAATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCgt          <  7:83250/149‑1 (MQ=255)                                                                                                                                              atTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACt    >  7:39989/1‑149 (MQ=255)                                                                                                                                               ttCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACt    >  1:84837/1‑148 (MQ=255)                                                                                                                                                 cGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACTGa  <  1:128065/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCGGGGGCGCGGCTGCTGGTGGACCCGACCCGGGTCACCTGCGGCCTGCTGGACAACCTGGCGGCAGAGGTGGTTCTTGTCGAGGGGCTAAACCCCACCACGCTGAGCAAGTCCTGCAAGGGTGACGACGACCTGGTGCATATTCGCCAGGTCATGGAGC‑‑AGGACGGCGCAGCCTTGTGCGAGTTTTTCGCCTGGTTCGAGGCGAACCTCGGCCGGGAGGTCATCACCGAGCTGACGGTTGACGAGCAATTGAGTGCCGCGCGTGCTCGACGGCCGAATTTCGTGTCACTGA  >  NC_002947/2545803‑2546094

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