Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,741,273:1	+C	76.9%	intergenic (+464/+288)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,273	1	.	C	76.9%	132.9 / 39.7	52	intergenic (+464/+288)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (6/6);  new base C (21/19);  total (27/25)
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.89e-01

CAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑T‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCA  >  NC_002947/4741124‑4741408                                                                                                                                                       |                                                                                                                                              cAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgag                                                                                                                                                <  2:1963350/150‑1 (MQ=255) cAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgag                                                                                                                                                <  2:1966072/150‑1 (MQ=255)  aGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagc                                                                                                                                               >  2:786799/1‑149 (MQ=255)  aGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagc                                                                                                                                               <  2:2378549/150‑2 (MQ=255)                                gTCAGGCTGCCATCCGCCTCCTGGACGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CC‑GCGACT‑TGCATTCGTGaatccag                                                                                                            <  2:1196155/150‑6 (MQ=255)                                 tCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGaatccagg                                                                                                           <  1:2199534/150‑7 (MQ=255)                                         ccATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACATTCGTGCATGAG‑CA‑GC‑TGCG‑GCGATT‑TGCATTCGTTAA‑TCCAGGCACAGGcc                                                                                                  >  2:260724/1‑150 (MQ=255)                                                 ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATGCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                          >  1:1451634/1‑150 (MQ=255)                                                  tcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACg                                                                                         >  2:296201/1‑150 (MQ=255)                                                   cctTGATGGTCGTAATGTGCAGCACGGAATCACGAGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGt                                                                                        >  1:332049/1‑150 (MQ=255)                                                    ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc                                                                                       <  1:1961742/150‑1 (MQ=255)                                                         tGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTcctcct                                                                                  >  2:2415303/1‑150 (MQ=255)                                                             cGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAg                                                                              >  1:269999/1‑150 (MQ=255)                                                              gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgc                                                                             <  2:331965/150‑1 (MQ=255)                                                              gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgc                                                                             >  1:331965/1‑150 (MQ=255)                                                                aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgcgc                                                                           >  1:592224/1‑150 (MQ=255)                                                                     gcagcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTT‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                      <  2:269999/150‑1 (MQ=255)                                                                        gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCGAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAcat                                                                   >  1:1094181/1‑148 (MQ=255)                                                                         caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACA‑CGTTGTCTGGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                  <  2:1165220/150‑1 (MQ=255)                                                                          aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                 <  1:786799/150‑1 (MQ=255)                                                                           cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                <  1:2416616/150‑1 (MQ=255)                                                                                  aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                         >  2:120772/1‑150 (MQ=255)                                                                                    gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                       >  2:175499/1‑150 (MQ=255)                                                                                      ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                     <  1:798856/150‑1 (MQ=255)                                                                                         gTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTAACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                  >  1:118490/1‑150 (MQ=255)                                                                                           cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                <  1:1293308/150‑1 (MQ=255)                                                                                             ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                              >  2:2407899/1‑150 (MQ=255)                                                                                             ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                              >  1:1951835/1‑150 (MQ=255)                                                                                               ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                            <  1:2116771/150‑1 (MQ=255)                                                                                                      gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                     <  2:592224/150‑1 (MQ=255)                                                                                                         tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                  >  2:2388218/1‑150 (MQ=255)                                                                                                            tGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCGTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                               <  2:29412/150‑1 (MQ=255)                                                                                                              aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                             <  1:2415303/150‑1 (MQ=255)                                                                                                                cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                           <  2:2247397/150‑1 (MQ=255)                                                                                                                 tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                          >  1:1865906/1‑150 (MQ=255)                                                                                                                  gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                         <  1:766812/150‑1 (MQ=255)                                                                                                                   ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                        <  2:2214123/150‑1 (MQ=255)                                                                                                                     tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                           >  2:1547492/1‑113 (MQ=255)                                                                                                                     tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                           <  1:1547492/113‑1 (MQ=255)                                                                                                                     tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGt                      >  2:1222072/1‑150 (MQ=255)                                                                                                                       taGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTAc                    <  1:260724/150‑1 (MQ=255)                                                                                                                         ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                  >  1:79963/1‑150 (MQ=255)                                                                                                                          gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                        <  2:1794721/95‑1 (MQ=255)                                                                                                                          gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                        >  1:1794721/1‑95 (MQ=255)                                                                                                                          gTCCAGCTCTGGTCACCTTCGTGCTTCTG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                 >  1:665188/1‑150 (MQ=255)                                                                                                                           tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                <  2:1109011/102‑1 (MQ=255)                                                                                                                           tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                >  1:1109011/1‑102 (MQ=255)                                                                                                                             cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                  >  2:1782486/1‑146 (MQ=255)                                                                                                                             cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                  <  1:1782486/146‑1 (MQ=255)                                                                                                                              aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAcc             >  1:275194/1‑150 (MQ=255)                                                                                                                              aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAcc             >  2:306901/1‑150 (MQ=255)                                                                                                                                ctctGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                      >  2:2377254/1‑91 (MQ=255)                                                                                                                                ctctGGTCACCTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                      <  1:2377254/91‑1 (MQ=255)                                                                                                                                ctctGGTCACCTTCGTGCTTGAG‑CA‑GC‑T‑TCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg           <  2:332049/150‑1 (MQ=255)                                                                                                                                      tCACCTTCGTGCTTGAG‑CAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc     >  1:2131376/1‑150 (MQ=255)                                                                                                                                         ccTTCGTGCTTGAGCAG‑CT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCa  <  2:118490/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                              CAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑T‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCA  >  NC_002947/4741124‑4741408

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