Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,511,522	G→A	8.5%	pseudogene (87/1851 nt)	PP_1325 ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,511,522	0	G	A	8.5%	152.0 / 2.6	56	pseudogene (87/1851 nt)	PP_1325	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (2/3);  total (28/28)
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.97e-01

GAATGCGTGCGGCGCGTGGGAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCAATTGTACCCAAGCACTGGCCGGG  >  NC_002947/1511374‑1511666                                                                                                                                                     |                                                                                                                                                   aaaggcgtgcgGCGCGTGGGGCGTCCTTCTGCTCGCTGTTCCGGTCTGCTGCGCTCATGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGCGGTGCGTGGTAGTTCGc                                                                                                                                                   <  2:836001/146‑1 (MQ=255)    tgcgtgcgGCGCGTGGGAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCa                                                                                                                                                <  2:1256517/150‑1 (MQ=255)    tgcgtgcgGCGCGTGGGAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCa                                                                                                                                                <  2:1256133/150‑1 (MQ=255)                 tGGGAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAggg                                                                                                                                   >  2:1526993/1‑150 (MQ=255)                   ggAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCg                                                                                                                                 <  2:1267480/150‑1 (MQ=255)                    gAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGa                                                                                                                                >  2:772638/1‑150 (MQ=255)                          attCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGgc                                                                                                                          >  1:957121/2‑150 (MQ=255)                                   gTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGAATTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGAGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CCGGcgac                                                                                                                >  2:317449/1‑147 (MQ=255)                                       gCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCACCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCcagc                                                                                                            >  2:1332794/1‑150 (MQ=255)                                       gCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCACCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCcagc                                                                                                            >  2:1332695/1‑150 (MQ=255)                                        ccAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGa                                                                                                                                >  2:1384505/1‑130 (MQ=255)                                        ccAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGa                                                                                                                                <  1:1384505/130‑1 (MQ=255)                                           ggTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAggg                                                                                                        <  2:696448/150‑1 (MQ=255)                                            gTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGggc                                                                                                       >  2:376091/1‑150 (MQ=255)                                              cggcggCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCACCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGggcgg                                                                                                     <  1:1332695/150‑1 (MQ=255)                                              cggcggCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCACCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGggcgg                                                                                                     <  1:1332794/150‑1 (MQ=255)                                                  ggcgCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAg                                                                                                 <  2:495982/150‑1 (MQ=255)                                                  ggcgCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAg                                                                                                 <  1:1167173/150‑1 (MQ=255)                                                   gcgcTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAgg                                                                                                <  2:1269854/150‑1 (MQ=255)                                                       tCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCaga                                                                                            >  2:1349100/1‑150 (MQ=255)                                                             gcagcaGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGAAGGTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCtg                                                                                      <  1:1251274/150‑1 (MQ=255)                                                              cagcagCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCtgt                                                                                     <  1:317449/150‑1 (MQ=255)                                                               agcagcGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCtgtg                                                                                    >  2:79796/1‑150 (MQ=255)                                                                gcagcGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAAC‑GCCGTCAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGa                                                                                   >  1:836525/1‑150 (MQ=255)                                                                  agcGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGag‑c                                                                                <  2:786207/150‑1 (MQ=255)                                                                      cggcgTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGAAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGcc                                                                            <  1:878577/150‑1 (MQ=255)                                                                       ggcgTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCg                                                                           >  2:1433244/1‑150 (MQ=255)                                                                            cTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCT‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGc                                                                      >  1:1439725/1‑150 (MQ=255)                                                                                     aCTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCa                                                             >  2:1534814/1‑150 (MQ=255)                                                                                       ttGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATg                                                           >  1:948916/1‑150 (MQ=255)                                                                                       ttGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATg                                                           >  1:1312568/1‑150 (MQ=255)                                                                                        tGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGt                                                          <  2:957121/150‑1 (MQ=255)                                                                                        tGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGt                                                          <  2:1312568/150‑1 (MQ=255)                                                                                             gAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAgg                                                     >  2:1442373/1‑150 (MQ=255)                                                                                               aGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTc                                                   <  1:446545/150‑1 (MQ=255)                                                                                                  tgCCGTTTCGAGCGGTTGCTCGGTGTTGGCGTTCGTCGTTCGCGGTGTAGCAGCGATGCGGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTg                                                <  2:836525/149‑1 (MQ=255)                                                                                                   gCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGt                                               >  1:1289491/1‑150 (MQ=255)                                                                                                    ccTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAg                                                                                                                                     <  1:883580/65‑1 (MQ=255)                                                                                                    ccTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAg                                                                                                                                     >  2:883580/1‑65 (MQ=255)                                                                                                    ccTTTTCGAGCAGTTGCTCGATGCTGGAAGCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGtt                                              >  1:1533324/1‑150 (MQ=255)                                                                                                     cTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAg                                                                                                 <  2:1504906/99‑1 (MQ=255)                                                                                                     cTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAg                                                                                                 >  1:1504906/1‑99 (MQ=255)                                                                                                           gAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCACATGATCATGTAAAGGTCCTGTTACTCGAt                                       <  1:1534814/150‑1 (MQ=255)                                                                                                              cAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCaa                                    >  1:1317273/1‑150 (MQ=255)                                                                                                                  tGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTg                                <  1:502519/150‑1 (MQ=255)                                                                                                                       gATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCgg                           >  2:1559965/1‑150 (MQ=255)                                                                                                                         tGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCa                         >  2:20221/1‑150 (MQ=255)                                                                                                                          gCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCaa                        <  2:1533324/150‑1 (MQ=255)                                                                                                                          gCTGGCGTCCGGGGTGCGTGGTAGTCCGCCCAGAATGGTCGGTTGCTAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCaa                        <  1:772638/150‑1 (MQ=255)                                                                                                                             ggCGTCCGGGGAAAGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCAATTg                     >  2:919546/1‑150 (MQ=255)                                                                                                                                          cGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAg                                                      >  1:1009137/1‑104 (MQ=255)                                                                                                                                          cGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAg                                                      <  2:1009137/104‑1 (MQ=255)                                                                                                                                          cGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCAATTGTACCCAAGCACTg        <  2:1289491/150‑1 (MQ=255)                                                                                                                                            tGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCACGTTGGCCGGCAATTGTACCCAAGCACTGGc      >  2:443167/1‑150 (MQ=255)                                                                                                                                              gTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCAATTGTACCCAAGCACTGGCCg    >  2:461009/1‑150 (MQ=255)                                                                                                                                                aGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGAAGAAGGCTACGAGTCAAAC‑AAGCCAAGCGATCACTTAAAGGTTCTGTTACTCGATCAAGTTGGCCGGCAATTGTACCCAAGCACTGGCCggg  <  1:1433244/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GAATGCGTGCGGCGCGTGGGAAGTCCTTCTGCTTGTAGGCCAGGTCGGCGGCGCTCAGGCGCAGCAGCGCGGCGTCTTCCGCAGACTTGCTGGAAGCGGCCTTTTCGAGCAGTTGCTCGATGCTGGCGTCCGGGGTGCGTGGTAGTTCGCCCAGGCTGGATGAGGGCGAGCTGGCG‑CAGGCTGCCAGCAGGGCGGCGAGGCAGAGGGCTGTGAG‑CAGCCGCAGGCAAGCGATCATGTAAAGGTCCTGTTACTCGATCAAGTTGGCCGGCAATTGTACCCAAGCACTGGCCGGG  >  NC_002947/1511374‑1511666

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