Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,633,545	T→G	6.9%	V165G (GTC→GGC)	lepB →	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,633,545	0	T	G	6.9%	253.3 / 5.6	83	V165G (GTC→GGC)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base G (2/4);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTG  >  NC_002947/1633398‑1633692                                                                                                                                                      |                                                                                                                                                    aCAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTc                                                                                                                                                    <  3:77456/149‑1 (MQ=255)  cAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTc                                                                                                                                                    >  3:276906/1‑148 (MQ=255)   aagaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGTCCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACAACAACGACAAGCGGCTGTTCGTCaa                                                                                                                                                  <  5:317057/149‑1 (MQ=255)   aagaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCaa                                                                                                                                                  >  7:3409/1‑149 (MQ=255)   aagaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCAGCTACACCACCGACCAGCGGCTGTTCGTCaa                                                                                                                                                  >  5:312096/1‑149 (MQ=255)     gaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCTTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACg                                                                                                                                                <  1:284342/149‑1 (MQ=255)     gaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACg                                                                                                                                                >  6:258366/1‑149 (MQ=255)       agGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGc                                                                                                                                              <  5:295319/149‑1 (MQ=255)           cATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACgggctgg                                                                                                                                          >  7:185764/1‑144 (MQ=255)               gAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGt                                                                                                                                      <  3:170150/149‑1 (MQ=255)                   tCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGcc                                                                                                                                  <  6:356599/149‑1 (MQ=255)                       tGACCCGCAACGCGTGGATGTAATGGTGTTCCGCTACCCGATTGACCCGC‑ATGTCAACTACATCAAGCGTGTGGGT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaac                                                                                                                              <  6:312096/149‑1 (MQ=255)                         aCCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaaca                                                                                                                             <  7:213984/148‑1 (MQ=255)                         aCCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaaca                                                                                                                             <  8:185764/148‑1 (MQ=255)                         aCCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaaca                                                                                                                             <  5:258366/148‑1 (MQ=255)                         aCCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaaca                                                                                                                             <  6:72466/148‑1 (MQ=255)                         aCCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGGGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaacag                                                                                                                            <  3:286431/149‑2 (MQ=255)                           ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                          >  6:360056/1‑149 (MQ=255)                           ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                          >  5:252182/1‑149 (MQ=255)                             gCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtgg                                                                                                                        >  5:106474/1‑149 (MQ=255)                               aaCGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACTTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGACCCACtggtc                                                                                                                      >  8:155518/1‑148 (MQ=255)                               aaCGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                      <  4:168780/149‑1 (MQ=255)                                aCGCGGTGATGTAATGGTGTTCCGCTACACGAGTGACCCTA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCTACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtgg                                                                                                                     <  8:131091/149‑1 (MQ=255)                                     gggATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGc                                                                                                                >  6:123209/1‑149 (MQ=255)                                       gATGTACTGGGGTTCCGCTACCCGAGTGACCCGA‑CTGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGa                                                                                                              <  4:222243/149‑1 (MQ=255)                                       gATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACTACTGGTGGGCAGCGa                                                                                                              <  3:335592/149‑1 (MQ=255)                                       gATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggt                                                                                                                       >  2:69662/1‑140 (MQ=255)                                       gATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggt                                                                                                                       <  1:69662/140‑1 (MQ=255)                                            aaTGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                         >  2:100304/1‑149 (MQ=255)                                             aTGGTGTTCCGCTACCCGAGTGACCCGAGATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                         >  1:23400/1‑149 (MQ=255)                                             aTGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCACCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCggg                                                                                                        >  3:250104/1‑149 (MQ=255)                                              tGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGt                                                                                                       <  6:106474/149‑1 (MQ=255)                                                gtgtTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggt                                                                                                                       <  5:208958/131‑1 (MQ=255)                                                gtgtTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggt                                                                                                                       >  6:208958/1‑131 (MQ=255)                                                  gtTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAccc                                                                                                   >  6:17995/1‑149 (MQ=255)                                                   ttCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTTGGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAccc                                                                                                   >  1:278118/1‑149 (MQ=255)                                                     ccGCTACCCGAGGGACCCGA‑ATGTCAACTACATCACGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                 <  8:78468/148‑1 (MQ=255)                                                       gCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGc                                                                                              >  7:264880/1‑149 (MQ=255)                                                           cccGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGt                                                                                           <  8:112399/148‑1 (MQ=255)                                                              gAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACCCCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTCCCCTGGGCAGTGCt                                                                                        >  8:174450/1‑148 (MQ=255)                                                              gAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTc                                                                                       <  5:17995/149‑1 (MQ=255)                                                               aGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCa                                                                                      <  2:126055/149‑1 (MQ=255)                                                                gTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAg                                                                                     >  7:17971/1‑149 (MQ=255)                                                                   aCCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGGCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCt                                                                                   <  5:360056/148‑1 (MQ=255)                                                                    cccGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTc                                                                                       <  6:165033/143‑1 (MQ=255)                                                                    cccGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTc                                                                                       >  5:165033/1‑143 (MQ=255)                                                                     ccGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAg                                                                                            <  4:265000/137‑1 (MQ=255)                                                                     ccGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAg                                                                                            >  3:265000/1‑137 (MQ=255)                                                                     ccGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCAGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACGGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                >  1:242939/1‑149 (MQ=255)                                                                      cGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                >  2:55021/1‑148 (MQ=255)                                                                       gA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACa                                                                              <  2:278118/149‑1 (MQ=255)                                                                           tGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAgg                                                                           <  8:3409/149‑1 (MQ=255)                                                                             tCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGa                                                                          >  5:226012/1‑148 (MQ=255)                                                                               aaCTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGaaaa                                                                       >  5:57770/1‑149 (MQ=255)                                                                                   aCATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTc                                                                   <  1:354561/149‑1 (MQ=255)                                                                                    cATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCg                                                                  >  3:332976/1‑149 (MQ=255)                                                                                    cATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCg                                                                  <  6:252182/149‑1 (MQ=255)                                                                                    cATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCg                                                                  <  5:162596/149‑1 (MQ=255)                                                                                     aTCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCgg                                                                 <  3:271195/149‑1 (MQ=255)                                                                                      tCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgc                                                                <  7:351040/149‑1 (MQ=255)                                                                                             tgtgTTT‑GGCCTGCCGGTCGACGTTGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCg                                                         <  3:238362/149‑1 (MQ=255)                                                                                                  tt‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAacac                                                    <  5:123209/149‑1 (MQ=255)                                                                                                       ccTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGa                                                >  5:48668/1‑149 (MQ=255)                                                                                                           ccGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATcc                                             <  4:332976/148‑1 (MQ=255)                                                                                                             gggCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                           >  4:220746/1‑148 (MQ=255)                                                                                                              ggCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCaa                                         <  2:23400/149‑1 (MQ=255)                                                                                                                   cGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGaaa                                    >  7:201443/1‑149 (MQ=255)                                                                                                                     tGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                   >  4:1650/1‑148 (MQ=255)                                                                                                                      ggTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                   >  2:79374/1‑147 (MQ=255)                                                                                                                       gTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGa                                 <  6:57770/148‑1 (MQ=255)                                                                                                                          cGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCg                             >  7:285518/1‑149 (MQ=255)                                                                                                                          cGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCg                             >  8:18037/1‑149 (MQ=255)                                                                                                                             tACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACAGGTGGGCAGCGAGCCGGGTACACTGGGCAGTGCTAAGCTGTAAAAAAAAAACCTCGGCGCGGACGAACACCTGATCCGCAAGGAAATGAGCCGTTa                          >  4:48215/1‑149 (MQ=255)                                                                                                                              acacCAACGACAAGCGGCTGTTCGTAAAAGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTa                          >  8:33871/1‑148 (MQ=255)                                                                                                                                   aaCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGTGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCAt                    <  6:347209/149‑1 (MQ=255)                                                                                                                                      gACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATgcc                 <  4:250104/149‑1 (MQ=255)                                                                                                                                           gCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCgg            >  1:19947/1‑149 (MQ=255)                                                                                                                                             ggCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGAc          <  8:17971/149‑1 (MQ=255)                                                                                                                                               cTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACca        <  1:55021/149‑1 (MQ=255)                                                                                                                                               cTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACca        <  6:226012/149‑1 (MQ=255)                                                                                                                                                tGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcag       <  8:201443/149‑1 (MQ=255)                                                                                                                                                  ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcag       >  6:128517/1‑147 (MQ=255)                                                                                                                                                     gTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg  >  5:370352/1‑149 (MQ=255)                                                                                                                                                      tCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg  >  1:85154/1‑148 (MQ=255)                                                                                                                                                      tCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGATCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGACAACCTGATCCGCAAGGAAATGAGCCGTAATCGCATGCCGCCGGACCAGCAGTg  >  6:267676/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGA‑ATGTCAACTACATCAAGCGTGTGGTT‑GGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTG  >  NC_002947/1633398‑1633692

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