Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,060,099	Δ1 bp	12.3%	coding (242/891 nt)	iciA →	chromosome initiation inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,060,099	0	A	.	12.3%	162.6 / 18.9	57	coding (242/891 nt)	iciA	chromosome initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base . (4/3);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

TTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAA  >  NC_002947/1059955‑1060247                                                                                                                                                 |                                                                                                                                                        tttCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTAc                                                                                                                                                    <  1:707657/150‑1 (MQ=255)    cccAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑gg                                                                                                                                                <  1:161239/150‑1 (MQ=255)     ccAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑Ggc                                                                                                                                               >  1:568567/1‑150 (MQ=255)                  ttGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCC‑GGTACCTGgcgc                                                                                                                                             >  2:660850/1‑139 (MQ=255)                  ttGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCC‑GGTACCTGgcgc                                                                                                                                             <  1:660850/139‑1 (MQ=255)                    gCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAggg                                                                                                                                >  2:449235/1‑150 (MQ=255)                    gCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAggg                                                                                                                                <  2:479821/150‑1 (MQ=255)                      ttGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACTGTGGGCTGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCa                                                                                                                              >  1:146227/1‑150 (MQ=255)                      ttGAGGCGCGGGACGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGAGCCAGGTACC‑GGCGCTGGACGAAGAGGGCa                                                                                                                              >  1:669340/1‑150 (MQ=255)                           gcgcGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGACGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑Ggcgc                                                                                                                                             <  1:749795/130‑1 (MQ=255)                           gcgcGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑Ggcgc                                                                                                                                             >  2:749795/1‑130 (MQ=255)                             gcgGGTCGGCCAGACGGTGCCGGTACGGGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGa                                                                                                                       >  2:595761/1‑150 (MQ=255)                                 gTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGCGGGCATGCCGGAgcgc                                                                                                                   >  2:275472/1‑150 (MQ=255)                                            ggTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGc                                                                                                        <  2:1099319/150‑1 (MQ=255)                                             gTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGcc                                                                                                       >  2:892132/1‑150 (MQ=255)                                             gTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGcc                                                                                                       >  2:816443/1‑150 (MQ=255)                                              tGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGcc‑c                                                                                                     <  2:125819/150‑1 (MQ=255)                                                     aCGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGc                                                                                              <  1:533656/150‑1 (MQ=255)                                                       gTGCCTCGCCGCCCGGCGCGACCGAAGGCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCg                                                                                            >  2:40158/1‑150 (MQ=255)                                                        tGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGa                                                                                           >  2:200066/1‑150 (MQ=255)                                                              gccgccCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCt                                                                                    >  2:747240/1‑150 (MQ=255)                                                               ccgccCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTg                                                                                   <  1:459858/150‑1 (MQ=255)                                                                cgccCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTgg                                                                                  <  1:557922/150‑1 (MQ=255)                                                                   ccGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGa                                                                                                                                         <  2:1032730/94‑1 (MQ=255)                                                                   ccGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGa                                                                                                                                         >  1:1032730/1‑94 (MQ=255)                                                                             cGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGAt                                                                                          <  1:613480/130‑1 (MQ=255)                                                                             cGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGAt                                                                                          >  2:613480/1‑130 (MQ=255)                                                                              gAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGcc                                                                    >  2:255538/1‑150 (MQ=255)                                                                               aaGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCg                                                                   <  1:255538/150‑1 (MQ=255)                                                                                   cGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGgcgc                                                               <  1:200066/150‑1 (MQ=255)                                                                                    ggCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGAAGGAGAGGGCATGCCGGAGCGCCTCAGGATCGTC‑CTCAACGGCGATAG‑CCTGGCCACCTGGTGGGCCGGcgcg                                                              >  1:870583/1‑150 (MQ=255)                                                                                        gccAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGttg                                                          <  2:568567/150‑3 (MQ=255)                                                                                        gccAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTgg                                                          <  1:595761/150‑1 (MQ=255)                                                                                           aGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCa                                                       >  2:532139/1‑150 (MQ=255)                                                                                            gCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCaa                                                      >  1:1057705/1‑150 (MQ=255)                                                                                              tACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGc                                                                     >  1:355799/1‑133 (MQ=255)                                                                                              tACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGc                                                                     <  2:355799/133‑1 (MQ=255)                                                                                                     ccATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTgcgc                                             >  1:369350/1‑150 (MQ=255)                                                                                                      cATGTGCAGCAGTTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGca                                            >  1:341983/1‑150 (MQ=255)                                                                                                              gcagGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAACAGGCGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGaa                                    <  1:449235/150‑1 (MQ=255)                                                                                                                 ggTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGt                                 <  2:369350/150‑1 (MQ=255)                                                                                                                   tGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGACTTCTTCTGCGCACAACAGAACGTGc                               <  1:766666/150‑1 (MQ=255)                                                                                                                        cTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTgacc                          <  1:892132/150‑1 (MQ=255)                                                                                                                                  gcgGCCGGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGACTGAGCATCGCC‑CTAAACGAAGATAG‑CCTGGCCACCTAGTAGGTCGTCGCAGTTGGGAAGTGCAGCGGACAGCAGAACGGTCTGAGCGACCtggtgg                >  1:777168/1‑150 (MQ=255)                                                                                                                                    gACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCCGTGGGCAACTTCTGCGCACAACACAACGTGCTGACCGACCTGGTGGtt              <  2:669340/150‑1 (MQ=255)                                                                                                                                    gACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGCAGGAGGTGTG‑‑TCGCCTCACAATTCGGGTTGTTCTGTCCACCGGTTGGTCCGGCGCGGTGGGCAACTTCTGCTCACAACAGAACGTGCTGACCGACCTGGTGGtt              <  1:275472/150‑1 (MQ=255)                                                                                                                                      ccTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGa            >  1:834582/1‑150 (MQ=255)                                                                                                                                        tccAGCGCC‑GGTACCTGGCGCTGCACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAgg          >  2:343037/3‑150 (MQ=255)                                                                                                                                        tccAGCGCC‑GGTACCTGGCGCTGCACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAgg          <  1:343037/148‑1 (MQ=255)                                                                                                                                         gCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCt                              <  1:41111/129‑1 (MQ=255)                                                                                                                                         gCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCt                              >  2:41111/1‑129 (MQ=255)                                                                                                                                         ccAGCGCC‑GGTACCTGGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAGGa         >  2:1048770/2‑150 (MQ=255)                                                                                                                                          cAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAGGAc        <  1:40158/150‑1 (MQ=255)                                                                                                                                           aGCGCC‑GGTACCTGGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACGACAGAACGTGCTGACCGACCtggtgg                >  2:645498/1‑141 (MQ=255)                                                                                                                                           aGCGCC‑GGTACCTGGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCtggtgg                <  1:645498/141‑1 (MQ=255)                                                                                                                                             cgcCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAGGACCAg     <  1:730198/150‑1 (MQ=255)                                                                                                                                                cAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGaa  <  1:816443/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                        TTTCCCAGCGCATCAAGTTGCTTGAGGCGCGGGTCGGCCAGCCGGTGCTGGTACGTGCCTCGCCGCCCGGCGCGACCGAAGTCGGCCGCCAGCTACTCAACCATGTGCAGCAGGTGCGCCTGCTCGAACGCGACCTGCAGCGCCAGGTACC‑GGCGCTGGACGAAGAGGGCATGCCGGAGCGCCTGCGCATCGCC‑CTCAACGCCGATAG‑CCTGGCCACCTGGTGGGCCGGCGCGGTGGGCAACTTCTGCGCACAACAGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAA  >  NC_002947/1059955‑1060247

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