Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,909,160	G→A	12.7%	Q45* (CAG→TAG)	PP_1709 ←	fumarylacetoacetate hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,909,160	0	G	A	12.7%	132.1 / 9.3	55	Q45* (CAG→TAG)	PP_1709	fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (5/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TCGCGGGCGATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCT  >  NC_002947/1909015‑1909306                                                                                                                                                   |                                                                                                                                                    tCGCGGGCGATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTgcggc                                                                                                                                                 >  2:800982/1‑150 (MQ=255) tCGCGGGCGATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTgcggc                                                                                                                                                 >  2:676469/1‑150 (MQ=255)         gctGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGgcgccgcgcc                                                                                                                                         >  2:1200653/3‑150 (MQ=255)         gATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGgcgccgcgcc                                                                                                                                         >  2:455395/1‑150 (MQ=255)               ggCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                   <  1:859182/150‑1 (MQ=255)                    cgccTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCg                                                                                                                              <  1:800982/150‑1 (MQ=255)                      ccTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGc                                                                                                                            >  1:655141/1‑150 (MQ=255)                        tttGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                   <  2:48792/141‑1 (MQ=255)                        tttGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                   >  1:48792/1‑141 (MQ=255)                            gCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACggg                                                                                                                      <  2:57364/150‑1 (MQ=255)                                tGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGt                                                                                                                  >  1:777892/1‑150 (MQ=255)                                    cgcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCg                                                                                                              >  2:499183/1‑150 (MQ=255)                                    cgcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCg                                                                                                              <  2:655141/150‑1 (MQ=255)                                     gcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCgg                                                                                                             >  1:849552/1‑150 (MQ=255)                                           aGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAgg                                                                                                       <  1:676469/150‑1 (MQ=255)                                             cTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                     <  1:1098806/150‑1 (MQ=255)                                              tCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCCCGTTGCTTGCCGGTTTCATGAAAAACC‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGccg                                                                                                    >  2:202294/1‑148 (MQ=255)                                                         gCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGcc                                                                                        >  1:1364193/1‑150 (MQ=255)                                                               cTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCgg                                                                                  <  2:945922/150‑1 (MQ=255)                                                                tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAg                                                                                                        <  1:872122/128‑1 (MQ=255)                                                                tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCggg                                                                                 >  1:1038015/1‑150 (MQ=255)                                                                tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAg                                                                                                        >  2:872122/1‑128 (MQ=255)                                                                   tcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAg                                                                              <  2:494965/150‑1 (MQ=255)                                                                    cgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGTATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAgc                                                                             <  1:175642/150‑1 (MQ=255)                                                                                     ggCAGTTCGCCCAGTGCCTCGACCACTTTGCTTGCCGGTTTCATGAAAAACACCAAGC‑CCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACGCCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCgg                                                            >  1:886386/1‑150 (MQ=255)                                                                                      gCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGa                                                           <  2:777892/150‑1 (MQ=255)                                                                                      gCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGCTTAATGAAACCAA‑CGCGCTCCTGCGGCGCCGCGCCGTAGGATTACGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGa                                                           >  1:1096857/1‑150 (MQ=255)                                                                                             gCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAg                                                    <  2:29989/150‑1 (MQ=255)                                                                                               ccAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                     <  1:1094457/100‑1 (MQ=255)                                                                                               ccAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                     >  2:1094457/1‑100 (MQ=255)                                                                                                      cTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGAt                                           >  2:347640/1‑150 (MQ=255)                                                                                                      cTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGAt                                           >  2:1271716/1‑150 (MQ=255)                                                                                                       tCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                          <  2:1096857/150‑1 (MQ=255)                                                                                                       tCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                          <  2:849552/150‑1 (MQ=255)                                                                                                              cGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCcggcg                                   <  1:1164827/150‑1 (MQ=255)                                                                                                                ttgcttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCggg                                 >  1:1215953/1‑150 (MQ=255)                                                                                                                ttgcttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCggg                                 <  1:1271716/150‑1 (MQ=255)                                                                                                                    ttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGGGGCCCATCCGGGATGCCCTCCTGG‑GCAGCTCAGCGGTGCCGCGGATCGGCAGGCTGACGATAGCCGGGGGGGTGa                             >  1:525498/1‑150 (MQ=255)                                                                                                                    ttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCg                                                        <  2:1018280/123‑1 (MQ=255)                                                                                                                    ttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCg                                                        >  1:1018280/1‑123 (MQ=255)                                                                                                                     tgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGaa                            <  1:1200653/150‑1 (MQ=255)                                                                                                                        cGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACAcc                                                                                             >  2:167251/1‑82 (MQ=255)                                                                                                                        cGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACAcc                                                                                             <  1:167251/82‑1 (MQ=255)                                                                                                                            ttCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTaa                     <  2:188869/150‑1 (MQ=255)                                                                                                                                tGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                          <  2:1413540/125‑1 (MQ=255)                                                                                                                                tGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                          >  1:1413540/1‑125 (MQ=255)                                                                                                                                tGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGtga                 >  1:1172274/1‑148 (MQ=255)                                                                                                                                tGAAAAACA‑CAGGCTCCTGCGGCGCCACGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCa                 >  1:831764/1‑150 (MQ=255)                                                                                                                                  aaaaaCA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATg               >  2:232724/1‑150 (MQ=255)                                                                                                                                      aca‑cAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTg           <  1:499183/150‑1 (MQ=255)                                                                                                                                        gccggaatCCTACGGCGCGGCGCCGCAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTgg          >  2:1184896/8‑150 (MQ=255)                                                                                                                                             gCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGAtcct     >  1:297358/1‑150 (MQ=255)                                                                                                                                             aatCCTACGGCGCGGCGCCGCAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGAtcct     >  1:1293549/3‑150 (MQ=255)                                                                                                                                                ccTGCGGCGCCGCGCCGTAGGATTCCGGCACCGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATcctcct  <  2:337599/150‑1 (MQ=255)                                                                                                                                                ccTACGGCGCGGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATcctcct  >  1:360228/1‑150 (MQ=255)                                                                                                                                                  tGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACa                          >  2:278747/1‑124 (MQ=255)                                                                                                                                                  tGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACa                          <  1:278747/124‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TCGCGGGCGATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCT  >  NC_002947/1909015‑1909306

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