Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,935,747	C→G	7.1%	R17G (CGC→GGC)	PP_1736 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,935,747	0	C	G	7.1%	219.6 / 6.3	70	R17G (CGC→GGC)	PP_1736	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base G (3/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATGA  >  NC_002947/1935599‑1935892                                                                                                                                                       |                                                                                                                                                  ccTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGccc                                                                                                                                                   <  7:99550/149‑1 (MQ=255)  cTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCcg                                                                                                                                                  >  7:230375/1‑149 (MQ=255)       ggCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGcgccg                                                                                                                                             <  2:270782/149‑1 (MQ=255)       ggCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGcgcc                                                                                                                                              <  8:232264/148‑1 (MQ=255)       ggCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGcgcc                                                                                                                                              <  2:208399/148‑1 (MQ=255)          gcACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCt                                                                                                                                          <  5:7333/149‑1 (MQ=255)          gcACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCt                                                                                                                                          <  5:2690/149‑1 (MQ=255)           cACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATAGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCGGGCgccgccgc                                                                                                                                         <  4:323680/149‑3 (MQ=255)            aCCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTAc                                                                                                                                        >  4:199937/1‑149 (MQ=255)             ccgccgGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGCTGCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTAc                                                                                                                                        >  1:263899/1‑149 (MQ=255)                   gTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCg                                                                                                                                 <  4:216260/149‑1 (MQ=255)                    ttCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAATACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCg                                                                                                                                 >  1:93883/1‑149 (MQ=255)                      cAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGTACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGgt                                                                                                                               <  4:224799/148‑1 (MQ=255)                       aaGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGgtgt                                                                                                                             >  3:182670/1‑149 (MQ=255)                          tCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCt                                                                                                                          >  1:184481/1‑149 (MQ=255)                             cgGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGc                                                                                                                       <  3:199937/149‑1 (MQ=255)                             cgGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGc                                                                                                                       <  2:263899/149‑1 (MQ=255)                              ggCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCg                                                                                                                      <  3:128569/149‑1 (MQ=255)                                    aGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCgcgc                                                                                                                                               <  2:335463/118‑1 (MQ=255)                                      ggTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGc                                                                                                              <  8:95236/149‑1 (MQ=255)                                          cgccgcGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGa                                                                                                           <  6:47115/148‑1 (MQ=255)                                            ccgcGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGtt                                                                                                        >  6:17029/1‑149 (MQ=255)                                             cgcgGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATc                                                                                                                >  4:22143/1‑140 (MQ=255)                                             cgcgGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATc                                                                                                                <  3:22143/140‑1 (MQ=255)                                                      aTGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGTCGGCGCCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCccgcc                                                                                              >  5:269273/1‑149 (MQ=255)                                                                gaagC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGAc                                                                                     >  6:9778/1‑148 (MQ=255)                                                                  agC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAAc                                                                                  >  2:219846/1‑149 (MQ=255)                                                                    c‑gcTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAAccc                                                                                <  7:180779/149‑1 (MQ=255)                                                                    c‑gcTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAAccc                                                                                >  7:111222/1‑149 (MQ=255)                                                                           aa‑aCACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTcc                                                                          <  2:184481/149‑1 (MQ=255)                                                                            a‑aCACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCg                                                                         >  3:283911/1‑149 (MQ=255)                                                                              acacATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCgg                                                                        <  4:283911/149‑1 (MQ=255)                                                                              acacATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCgg                                                                        <  1:177155/149‑1 (MQ=255)                                                                                acaTCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTc                                                                      <  8:230375/149‑1 (MQ=255)                                                                                     tACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCAtcgt                                                                 >  1:116352/1‑149 (MQ=255)                                                                                        cccAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCgg                                                              <  5:17029/149‑1 (MQ=255)                                                                                         ccAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGc                                                             <  2:93883/149‑1 (MQ=255)                                                                                          cAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCa                                                            <  1:219846/149‑1 (MQ=255)                                                                                          cAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCa                                                            >  3:222043/1‑149 (MQ=255)                                                                                            ggaaggaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCAcc                                                          <  4:222043/149‑1 (MQ=255)                                                                                             gaaggaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCACCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCt                                                         >  1:315021/1‑149 (MQ=255)                                                                                               aggaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTcc                                                       >  1:62290/1‑149 (MQ=255)                                                                                                 gaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATc                                                                                                                <  2:283354/90‑1 (MQ=255)                                                                                                 gaCGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATc                                                                                                                >  1:283354/1‑90 (MQ=255)                                                                                                    gCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCgg                                                  <  2:116352/149‑1 (MQ=255)                                                                                                     ccATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCggg                                                 <  8:111222/149‑1 (MQ=255)                                                                                                      cATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGGCAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCgggg                                                >  1:30792/1‑149 (MQ=255)                                                                                                            accaccCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATc                                           >  3:229181/1‑148 (MQ=255)                                                                                                                cccccGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCgggg                                                >  5:287249/1‑139 (MQ=255)                                                                                                                cccccGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCgggg                                                <  6:287249/139‑1 (MQ=255)                                                                                                                           cGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTgg                           >  8:184922/1‑149 (MQ=255)                                                                                                                              tctcATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGc                          >  6:852/1‑147 (MQ=255)                                                                                                                                  aTCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAgtgg                    >  2:294374/1‑149 (MQ=255)                                                                                                                                    ccTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAgtggtg                  >  8:248411/1‑149 (MQ=255)                                                                                                                                     cTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGa                                                                                      >  4:185354/1‑80 (MQ=255)                                                                                                                                     cTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGa                                                                                      <  3:185354/80‑1 (MQ=255)                                                                                                                                      tGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcg                >  6:32167/1‑149 (MQ=255)                                                                                                                                      tGTCTGGCGGCGGCGCCAGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCACTggtgcg                >  6:153061/1‑149 (MQ=255)                                                                                                                                       gTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCa                        <  3:6678/140‑1 (MQ=255)                                                                                                                                       gTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCa                        >  4:6678/1‑140 (MQ=255)                                                                                                                                        tCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcgg               <  7:248411/148‑1 (MQ=255)                                                                                                                                         tggtaggcggcgCGGGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcggtg             >  2:16846/6‑149 (MQ=255)                                                                                                                                          tGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcggtg             >  5:242198/1‑148 (MQ=255)                                                                                                                                           ggcggcggcGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCg                                    >  7:206332/1‑124 (MQ=255)                                                                                                                                           ggcggcggcGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCg                                    <  8:206332/124‑1 (MQ=255)                                                                                                                                            gcggcggcgCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCAc          >  8:39133/1‑149 (MQ=255)                                                                                                                                              ggcggcgCGGGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcggtg             <  6:79129/144‑1 (MQ=255)                                                                                                                                              ggcggcgCGGGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTggtgcggtg             >  5:79129/1‑144 (MQ=255)                                                                                                                                              ggcggcgCGGGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACtt        >  2:224415/1‑149 (MQ=255)                                                                                                                                              ggcggcgCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACtt        <  7:184922/149‑1 (MQ=255)                                                                                                                                                  gcgcCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAAt    >  6:271955/1‑149 (MQ=255)                                                                                                                                                    gcCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATGa  <  7:39133/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CCTGTGGGCGCACCGCCGGTTCAAGTCGCGGCCCGAGGGTTCGCCGCGGCTGTATGACAAGAAGAAGC‑GCTGAAA‑ACACATCTACCCCAGGAAGGACGCCATGGCACCACCCCCGGTAACCGGTCTCATCCTGTCTGGCGGCGGCGCCCGCGCCGCCTACCAGGTCGGTGTGCTGGCGGGCATCGCCGAGTTGCTGCCCGCCGGCGCGGACAACCCTTTTCCGGTCATCGTCGGCACCTCCGCCGGGGCCATCAACGCCGTCACCCTGGCCAGTGGTGCGGTGCACTTCAATGA  >  NC_002947/1935599‑1935892

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