Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,009,722	T→G	5.9%	intergenic (+82/‑90)	PP_4412 → / → PP_4413	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,009,722	0	T	G	5.9%	193.6 / 3.7	68	intergenic (+82/‑90)	PP_4412/PP_4413	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (2/2);  total (34/34)
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.33e-01

TGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTTCTCTC  >  NC_002947/5009583‑5009866                                                                                                                                               |                                                                                                                                                 tGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCt                                                                                                                                         >  6:197139/1‑149 (MQ=255)  gggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCt                                                                                                                                         >  7:164534/1‑148 (MQ=255)   ggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTa                                                                                                                                       <  1:325844/149‑1 (MQ=255)   ggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGACGC‑C‑GCCATACGGCGAGCGGCTTggttca                                                                                                                                         <  1:47744/147‑5 (MQ=38)     cagcatgagGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGGGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCTGGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTAc                                                                                                                                      <  1:92235/149‑1 (MQ=255)     cagcagcagGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACg                                                                                                                                     <  2:109661/149‑1 (MQ=255)            agGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGAGTCCTGGGCGCGCCCTTCGTGCC‑GGCTGTCGTGAACCAAGCATCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCtac                                                                                                                               >  7:352752/1‑146 (MQ=255)                ttCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATcc                                                                                                                          >  1:25576/1‑149 (MQ=255)                 tCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTCGGGCCTAAGGCG‑GCGTGATTGCCCTTACCGGCttgccta                                                                                                                           >  1:125600/1‑143 (MQ=255)                    gACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTcgc                                                                                                                      <  6:80935/149‑1 (MQ=255)                      cGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGt                                                                                                                    <  8:164534/149‑1 (MQ=255)                       gATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGGGGTCGTGAACCAAGCCCCAC‑GCCGTATGGCC‑GCGTCTTGGCCCTTACGGGCTTCCAACCTCGCGtc                                                                                                                   >  4:192151/1‑149 (MQ=255)                                cctcctCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTc                                                                                                          >  7:223186/1‑149 (MQ=255)                                 ctcctcCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCACGCGTCTCCGCCTTcc                                                                                                         >  3:240639/1‑149 (MQ=255)                                     tccAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGc                                                                                                     >  4:298839/1‑149 (MQ=255)                                        aaaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgc                                                                                                  <  4:240639/149‑1 (MQ=255)                                          aCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAACAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCGC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgcgc                                                                                                <  5:197701/149‑1 (MQ=255)                                            ggggAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgcgcgc                                                                                              >  5:225086/1‑149 (MQ=255)                                            ggggAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTAGGGCC‑GGCGGTCGTGAACCAATCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTAAGGGCTTGCATACTAGCCGCGCCGCCTTCCTGGCTtcgcgc                                                                                              >  3:214860/1‑149 (MQ=255)                                                  ctctTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCt                                                                                        >  5:353859/1‑149 (MQ=255)                                                     ttctAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGGCTCCGCCTTCCTGGCTGCGCGCTCCGCGTgg                                                                                     >  4:43005/1‑148 (MQ=255)                                                               gCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAAc                                                                           >  6:380728/1‑149 (MQ=255)                                                                  aGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGt                                                                         >  2:122885/1‑148 (MQ=255)                                                                      ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCaaa                                                                     <  3:64587/148‑1 (MQ=255)                                                                      ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCaaa                                                                     <  3:192151/148‑1 (MQ=255)                                                                      ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAAt                                                                    >  7:13230/1‑149 (MQ=255)                                                                       tgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGGTCTGCCCAACGTCAAAt                                                                    >  1:349009/1‑148 (MQ=255)                                                                       tgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAAt                                                                    >  6:350412/1‑148 (MQ=255)                                                                         cAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGTTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  <  1:297911/148‑1 (MQ=255)                                                                         cAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  <  1:170567/148‑1 (MQ=255)                                                                         cAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAg                                                                 <  8:223186/149‑1 (MQ=255)                                                                          aaaaGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGc                                                                >  6:169458/1‑149 (MQ=255)                                                                            atgCGTCCTGGGCGGGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGTGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGt                                                              <  8:352752/147‑1 (MQ=255)                                                                            aaGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCCTCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  >  8:329286/1‑145 (MQ=255)                                                                            aaGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  <  7:329286/145‑1 (MQ=255)                                                                            aaGCGTCCTGGGCGCGCCCATCGGGCC‑GGCTGTGGTGAAACAACCCGCAC‑GCCGTATTGGG‑GCGTCGTTGCCCGTTAGGGCTTGCCTACTAGCGGGTCCGCCGTGCTGGCATCGCGCGCCGCGTGCAAGGGCCCACGTCAAATGAGCgg                                                              >  4:370784/1‑148 (MQ=255)                                                                             aGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGtc                                                             <  5:380728/149‑1 (MQ=255)                                                                              gCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGGCTCCGCCTGCCTGGCTGCGCGCTCCGCTTGCTCTGCCCACCGTCAAATGCGCGtct                                                            >  6:192398/1‑149 (MQ=255)                                                                                   cTGGGCGCGCGCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTGTGCCCAACGTCAAATGAGCGtctttctt                                                       <  3:266226/149‑1 (MQ=255)                                                                                   cTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGtctttctt                                                       <  7:123625/149‑1 (MQ=255)                                                                                      ggCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTg                                                      <  4:139119/147‑1 (MQ=255)                                                                                       gcgcgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGc                                                    >  8:133350/1‑148 (MQ=255)                                                                                       gcgcgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCa                                                   >  8:153118/1‑149 (MQ=255)                                                                                         gcgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATc                                                 <  3:43005/149‑1 (MQ=255)                                                                                          cgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCt                                                <  7:133350/149‑1 (MQ=255)                                                                                           gcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTg                                               >  4:195759/1‑149 (MQ=255)                                                                                              cTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCAt                                            <  3:298839/149‑1 (MQ=255)                                                                                                 cGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATg                                           <  2:25576/147‑1 (MQ=255)                                                                                                          cTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGtt                                 <  4:27605/149‑1 (MQ=255)                                                                                                            gTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTaa                               >  7:330225/1‑149 (MQ=255)                                                                                                            gTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTaa                               >  7:322814/1‑149 (MQ=255)                                                                                                             tCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAAc                              <  5:197139/149‑1 (MQ=255)                                                                                                                tGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCa                           >  7:270477/1‑149 (MQ=255)                                                                                                                    ccAAGACGC‑C‑GCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGTCCAACGTCAAATGAGCGTCTTTATTGGCATCTGCATGCCACCCTGTTAAATCActt                        >  2:336138/1‑148 (MQ=255)                                                                                                                    ccAAGACGC‑C‑GCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtc                       >  4:266954/1‑149 (MQ=255)                                                                                                                     cAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtct                      <  8:330225/149‑1 (MQ=255)                                                                                                                      aaGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtctt                     >  7:172753/1‑149 (MQ=255)                                                                                                                      aaGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtctt                     >  2:58445/1‑149 (MQ=255)                                                                                                                       aGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCttt                    <  6:225086/149‑1 (MQ=255)                                                                                                                        gCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTg                   <  1:58445/149‑1 (MQ=255)                                                                                                                        agacgcc‑GCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCttt                    <  3:266954/143‑1 (MQ=255)                                                                                                                                gCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTa            <  7:153118/149‑1 (MQ=255)                                                                                                                                gCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTa            <  3:370784/149‑1 (MQ=255)                                                                                                                                   gTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGcc         >  2:69072/1‑149 (MQ=255)                                                                                                                                        gcg‑tcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtct     <  8:322814/148‑1 (MQ=255)                                                                                                                                        gcg‑gcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtct     <  1:69072/148‑1 (MQ=255)                                                                                                                                         cg‑gcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtctc    >  4:35082/1‑148 (MQ=255)                                                                                                                                          g‑gcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtctctc  <  6:353859/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 TGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTTCTCTC  >  NC_002947/5009583‑5009866

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