Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,744,600	A→C	7.5%	S76A (TCC→GCC)	PP_4198 ←	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,744,600	0	A	C	7.5%	195.8 / 4.1	63	S76A (TCC→GCC)	PP_4198	OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base C (3/2);  total (35/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 = 1.00e+00

GCCAGTTTGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTG  >  NC_002947/4744452‑4744748                                                                                                                                                      |                                                                                                                                                      gCCAGTTTGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGa                                                                                                                                                       <  3:369172/149‑1 (MQ=255)  ccAGTTTGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGaa                                                                                                                                                      <  6:316050/149‑1 (MQ=255)         ggCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGc                                                                                                                                               >  2:176130/1‑149 (MQ=255)          gCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGCACCGGTTcc                                                                                                                                              <  4:145629/149‑1 (MQ=255)               aCGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGAt                                                                                                                                           <  8:356324/147‑1 (MQ=255)               aCGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATg                                                                                                                                          <  2:1993/148‑1 (MQ=255)                 gCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCGGAGGTAAGCCTAATGCCG‑CACATCGGCAGCGGTGGCGAGGtcg                                                                                                                                       >  4:79612/1‑147 (MQ=255)                    aGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTg                                                                                                                                    <  1:206851/149‑1 (MQ=255)                     ggCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTg                                                                                                                                    >  2:71668/1‑148 (MQ=255)                       ccAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                  >  2:49249/1‑148 (MQ=255)                        cAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                  >  7:205294/1‑147 (MQ=255)                         aGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAg                                                                                                                               <  4:173435/149‑1 (MQ=255)                          gCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAgg                                                                                                                              <  2:306351/149‑1 (MQ=255)                              cAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGgcg                                                                                                                            <  6:4457/147‑1 (MQ=255)                                 gcgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑gg                                                                                                                      <  7:14380/149‑1 (MQ=255)                                    ctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑Ggcgc                                                                                                                   <  2:120297/149‑1 (MQ=255)                                       gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGa                                                                                                                 <  8:205294/148‑1 (MQ=255)                                             cAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGc                                                                                                           <  1:176130/148‑1 (MQ=255)                                               ggTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGAc                                                                                                        >  3:297615/1‑149 (MQ=255)                                                 tCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACt                                                                                                       >  4:196840/1‑148 (MQ=255)                                                    aGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCAcc                                                                                                   <  4:15590/149‑1 (MQ=255)                                                            cTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGaa                                                                                           >  7:293977/1‑149 (MQ=255)                                                               gcgcTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAAtca                                                                                        >  2:187064/1‑149 (MQ=255)                                                                 gcTCCTGGTTCAGCGCCAGCTTGGCATTCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCTGGCGCGAGCCAGCGACTCACCCGCCCGAAtcac                                                                                       <  1:49249/149‑1 (MQ=255)                                                                          tCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACgcggc                                                                             >  8:10239/1‑149 (MQ=255)                                                                             gcgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGgcgcg                                                                          >  4:73197/1‑149 (MQ=255)                                                                              cgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTGGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCa                                                                                                     <  2:157203/121‑1 (MQ=255)                                                                              cgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTGGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCa                                                                                                     >  1:157203/1‑121 (MQ=255)                                                                                 cAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTgc                                                                      >  8:244468/1‑149 (MQ=255)                                                                                  aGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTgcg                                                                     <  3:196840/149‑1 (MQ=255)                                                                                  aGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTgcg                                                                     <  8:293977/149‑1 (MQ=255)                                                                                          aTGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGCATCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACAt                                                             >  7:104393/1‑149 (MQ=255)                                                                                                gggggCAATCTCGCTGTATCGCTGGGTGATGGAATCATAATGCCGCCTAAAGGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCg                                                           <  1:121623/144‑1 (MQ=255)                                                                                                 cgGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAat                                                       >  7:141230/1‑148 (MQ=255)                                                                                                    gCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCAGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCTCAGCACATCGGAatcatc                                                   >  1:103275/1‑149 (MQ=255)                                                                                                      aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGaa                                                        <  8:258092/142‑1 (MQ=255)                                                                                                      aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGaa                                                        >  7:258092/1‑142 (MQ=255)                                                                                                       aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCttt                                                <  7:10239/149‑1 (MQ=255)                                                                                                       aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCttt                                                <  4:297615/149‑1 (MQ=255)                                                                                                       aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCt                                                  >  1:131696/1‑147 (MQ=255)                                                                                                         ctcGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGa                                              >  1:185574/1‑149 (MQ=255)                                                                                                              tGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGAccccc                                         >  5:155759/1‑149 (MQ=255)                                                                                                                tAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTg                                       >  1:204472/1‑149 (MQ=255)                                                                                                                  gcgcTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTgg                                      <  1:187064/148‑1 (MQ=255)                                                                                                                  gcgcTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTgg                                      >  5:110095/1‑148 (MQ=255)                                                                                                                     ctggctgAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAgg                                  <  5:315352/149‑1 (MQ=255)                                                                                                                          tgAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCt                                                  <  3:194386/128‑1 (MQ=255)                                                                                                                          tgAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCt                                                  >  4:194386/1‑128 (MQ=255)                                                                                                                               tAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCTGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGc                        <  1:71668/149‑1 (MQ=255)                                                                                                                               tAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGc                        >  5:235463/1‑149 (MQ=255)                                                                                                                                aaGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGcc                       >  7:9301/1‑149 (MQ=255)                                                                                                                                  gCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCt                      >  5:311139/1‑148 (MQ=255)                                                                                                                                  gCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCt                      <  6:311139/148‑1 (MQ=255)                                                                                                                                   cATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCCCCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTc                     >  7:170838/1‑148 (MQ=255)                                                                                                                                      aaTGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATc                  >  4:219685/1‑148 (MQ=255)                                                                                                                                      aaTGCCG‑CACATCGGAACCGGTGCCGATGAAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATc                  >  5:205332/1‑148 (MQ=255)                                                                                                                                           cG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcgc            <  5:20063/149‑1 (MQ=255)                                                                                                                                            gttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTg                          <  5:241147/132‑1 (MQ=255)                                                                                                                                            gttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTg                          >  6:241147/4‑135 (MQ=255)                                                                                                                                            gttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcgc            >  4:266772/4‑149 (MQ=255)                                                                                                                                            g‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGcgcg           <  2:131696/149‑1 (MQ=255)                                                                                                                                              cacaTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCgg          >  5:222339/1‑149 (MQ=255)                                                                                                                                                 aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCAc       >  5:60972/1‑149 (MQ=255)                                                                                                                                                 aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAAGCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGCCCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCAc       >  8:288464/1‑149 (MQ=255)                                                                                                                                                    ggAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgc    >  1:221014/1‑149 (MQ=255)                                                                                                                                                      aaCCGGTTCCGATGTAACTGGACAGGCGCGC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCTCGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGTCCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg  <  8:104393/149‑1 (MQ=255)                                                                                                                                                      aaCCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg  <  8:267357/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCCAGTTTGGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCG‑CACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTC‑GGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTG  >  NC_002947/4744452‑4744748

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