Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,707,272	G→C	15.2%	intergenic (+48/+94)	PP_1501 → / ← PP_1502	hypothetical protein/OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,707,272	0	G	C	15.2%	131.8 / 17.6	53	intergenic (+48/+94)	PP_1501/PP_1502	hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base C (4/4);  total (24/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 = 7.20e-01

TTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGA  >  NC_002947/1707126‑1707419                                                                                                                                                    |                                                                                                                                                      ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggcg                                                                                                                                                    <  2:601783/150‑1 (MQ=255) ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggcg                                                                                                                                                    <  1:409107/150‑1 (MQ=255)     acacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                                <  1:1537852/150‑1 (MQ=255)       acCCAGCAGACAGGCCAGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                              >  2:1221938/1‑150 (MQ=255)         ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            >  1:499156/1‑150 (MQ=255)          cagcagAAAGGTCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                           >  1:1560443/1‑150 (MQ=255)          cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                           <  1:744143/150‑1 (MQ=255)            gcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGAGCTGCATGGGCGCGGACACAC‑CGCGGTCGGTGCGGTGGTGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                         <  2:688904/150‑1 (MQ=38)             cagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGAGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                                        >  1:1788106/1‑150 (MQ=255)                 aaGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGCGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                                                    >  1:7252/1‑150 (MQ=255)                          ccGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaa                                                                                                                           <  1:672127/150‑1 (MQ=255)                            ggCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGACACAC‑CGCGGCAGGTGCGCTAAAACTACTGCGGTGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaaca                                                                                                                         <  2:1005031/150‑1 (MQ=37)                                   tggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCgg                                                                                                                  <  1:1223345/150‑1 (MQ=255)                                    ggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                                 >  1:1872119/1‑150 (MQ=255)                                     gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCg                                                                                                                >  2:1110675/1‑150 (MQ=255)                                        gACAACCGCAAGCGCGCGGTGCGGTGCATGGGCGAGGAAACAC‑GTCGTCATGTGCAGTGATGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAAc                                                                                                             <  2:591629/150‑1 (MQ=38)                                           aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaaa                                                                                                          >  2:345175/1‑150 (MQ=255)                                           aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACaaa                                                                                                          >  2:537801/1‑150 (MQ=255)                                             ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                                        <  1:537801/150‑1 (MQ=255)                                              cGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaa                                                                                                       >  2:327166/1‑150 (MQ=255)                                                 aaGCGCGCGGTGCGCTGCATGGGCGAGGAAACACACGC‑GAAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCCCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAg                                                                                                    <  1:576114/150‑1 (MQ=38)                                                   gcgcgcGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            >  2:1376721/1‑108 (MQ=25)                                                   gcgcgcGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            <  1:1376721/108‑1 (MQ=25)                                                                      ggCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCt                                                                             >  1:404419/1‑150 (MQ=255)                                                                       gCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCtt                                                                            <  1:12223/150‑1 (MQ=255)                                                                        cGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGGAAAGCCCCGCGATCGA‑CGC‑GGGGCttt                                                                           <  2:1560841/150‑2 (MQ=38)                                                                           ggAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTcc                                                                        <  2:1881016/150‑1 (MQ=255)                                                                              aaCAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                     <  1:1183859/150‑1 (MQ=255)                                                                              aaCAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                     <  1:1110675/150‑1 (MQ=255)                                                                               acac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTc                                                                    <  2:404419/150‑1 (MQ=255)                                                                                cac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCa                                                                   >  1:1302791/1‑150 (MQ=255)                                                                                 ac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑gggg                                                                               >  1:1855173/1‑137 (MQ=255)                                                                                 ac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑gggg                                                                               <  2:1855173/137‑1 (MQ=255)                                                                                 ac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAg                                                                  >  1:997932/1‑150 (MQ=255)                                                                                 ac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAg                                                                  >  1:23407/1‑150 (MQ=255)                                                                                 ac‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAg                                                                  >  2:1517035/1‑150 (MQ=255)                                                                                    cgcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGACGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAGGCCCCGCG‑TCGATCGCGGGGGCTTGTCCTTTCAGc                                                                 >  1:167399/1‑150 (MQ=255)                                                                                     gcgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTa                                                               <  1:1221938/150‑1 (MQ=255)                                                                                                 gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCATCTAGCCAGCCTTAgc                                                   >  2:1844312/1‑150 (MQ=255)                                                                                                   gAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgcgc                                                 <  2:1788106/150‑1 (MQ=255)                                                                                                                  caTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAcc                                  <  2:997932/150‑1 (MQ=255)                                                                                                                    tttttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGa                                <  1:1813087/150‑1 (MQ=255)                                                                                                                        tGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCTATCTCACAGGCGCCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑GCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgt                            <  1:580437/150‑1 (MQ=255)                                                                                                                         ggCGCCCGCCATATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtg                           <  1:1516010/150‑1 (MQ=255)                                                                                                                                   agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCgg                 <  2:485099/150‑1 (MQ=255)                                                                                                                                   agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCgg                 <  2:485108/150‑1 (MQ=255)                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg              >  1:388667/1‑150 (MQ=39)                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATcgc‑g                                                                                  <  1:764594/82‑1 (MQ=255)                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATcgc‑g                                                                                  >  2:764594/1‑82 (MQ=255)                                                                                                                                       cGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                     <  2:1240473/94‑1 (MQ=38)                                                                                                                                       cGAGCGCCGCCCGCGCGGCGCTCGATCTCGCAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCttt                                                                     >  1:1240473/1‑94 (MQ=38)                                                                                                                                         aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt           <  1:345175/150‑1 (MQ=255)                                                                                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCGGGCTCACCGGGCGTGGTGCCACCGGCTGGTTGTCg       >  2:1535347/1‑150 (MQ=255)                                                                                                                                               gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt     <  1:1494178/150‑1 (MQ=255)                                                                                                                                               gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt     >  2:583766/1‑150 (MQ=255)                                                                                                                                                  cgggcggCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGga  <  2:499156/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACAC‑CGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGC‑GGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGA  >  NC_002947/1707126‑1707419

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