Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,633,545	T→G	7.1%	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	7.1%	259.9 / 9.0	84	V165G (GTC→GGC)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base G (3/3);  total (44/40)
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.34e-01

GACAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTG  >  NC_002947/1633397‑1633692                                                                                                                                                      |                                                                                                                                                    gACAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGt                                                                                                                                                     >  7:325505/1‑149 (MQ=255) gACAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGt                                                                                                                                                     <  3:39007/149‑1 (MQ=255)    aagaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCaa                                                                                                                                                  >  3:193426/1‑149 (MQ=255)    aagaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCaa                                                                                                                                                  >  4:85041/1‑149 (MQ=255)    aagaagGTCATCGAGGTCGGTGACCCACAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCaa                                                                                                                                                  >  1:199136/1‑149 (MQ=255)      gaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACg                                                                                                                                                <  8:24478/149‑1 (MQ=255)      gaagGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACg                                                                                                                                                <  4:294800/149‑1 (MQ=255)          gTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCACCGACACGCGGCTGTTCGTCACCGGc                                                                                                                                              >  4:226996/1‑147 (MQ=255)          gTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGc                                                                                                                                              >  2:302541/1‑147 (MQ=255)            cATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACgggctgg                                                                                                                                          >  8:342420/1‑144 (MQ=255)             aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGcc                                                                                                                                         >  2:56841/1‑149 (MQ=255)             aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGcc                                                                                                                                         >  1:140669/1‑149 (MQ=255)                gAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGt                                                                                                                                      <  3:141452/149‑1 (MQ=255)                            ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                          >  3:306416/1‑149 (MQ=255)                            ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                          >  3:101586/1‑149 (MQ=255)                            ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                          >  1:245381/1‑149 (MQ=255)                            ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACt                                                                                                                          >  1:240220/1‑149 (MQ=255)                              gCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtgg                                                                                                                        >  1:22583/1‑149 (MQ=255)                              gCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaacagccgc                                                                                                                        <  7:342420/149‑6 (MQ=39)                               cAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggt                                                                                                                       <  4:101586/149‑1 (MQ=255)                                    cgGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCa                                                                                                                  >  8:126435/1‑149 (MQ=255)                                     ggggATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGTCAGCCCGTTGCCGAACAACTGGTGGGCAg                                                                                                                 >  4:233025/1‑149 (MQ=255)                                        gATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGa                                                                                                              <  3:85041/149‑1 (MQ=255)                                        gATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGa                                                                                                              <  2:240220/149‑1 (MQ=255)                                        gATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGGGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGa                                                                                                              <  2:245381/149‑1 (MQ=255)                                          tGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGc                                                                                                            <  4:193426/149‑1 (MQ=255)                                          tGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGAGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGc                                                                                                            <  1:302541/149‑1 (MQ=255)                                            tAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGcc                                                                                                           <  2:344048/148‑1 (MQ=255)                                             aaTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGcc                                                                                                           <  1:121738/147‑1 (MQ=255)                                             aaTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGcc                                                                                                           >  2:121738/1‑147 (MQ=255)                                             aaTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                         >  8:14529/1‑149 (MQ=255)                                              aTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCggg                                                                                                        >  7:201805/1‑149 (MQ=255)                                              aTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCggg                                                                                                        >  5:52990/1‑149 (MQ=255)                                              aTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                         <  1:56841/148‑1 (MQ=255)                                                    ttCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCt                                                                                                  <  4:306416/149‑1 (MQ=255)                                                        gCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGc                                                                                              <  2:140669/149‑1 (MQ=255)                                                         cTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                             <  2:68927/149‑1 (MQ=255)                                                         cTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaac                                                                                                                              <  7:321278/116‑1 (MQ=39)                                                         cTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaac                                                                                                                              >  8:321278/1‑116 (MQ=39)                                                          tACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                               >  7:95413/1‑146 (MQ=255)                                                            cccGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                             <  8:95413/146‑1 (MQ=255)                                                            cccGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTg                                                                                          <  7:126435/149‑1 (MQ=255)                                                            cccGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTg                                                                                          <  8:325505/149‑1 (MQ=255)                                                               gAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                         <  7:313382/131‑1 (MQ=255)                                                               gAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                         >  8:313382/1‑131 (MQ=255)                                                                aGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCa                                                                                      <  2:22583/149‑1 (MQ=255)                                                                 gTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAg                                                                                     >  4:267220/1‑149 (MQ=255)                                                                    aCCCTAATGTCACCTAGCATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCt                                                                                   <  1:127895/149‑1 (MQ=255)                                                                     cccGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTg                                                                                  >  5:66712/1‑148 (MQ=255)                                                                      ccGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGt                                                                                 >  3:279318/1‑148 (MQ=255)                                                                      ccGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                <  3:299619/149‑1 (MQ=255)                                                                           tGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAgg                                                                           <  6:52990/149‑1 (MQ=255)                                                                                   a‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTc                                                                   >  1:168559/1‑149 (MQ=255)                                                                                     cATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCg                                                                  <  4:279318/149‑1 (MQ=255)                                                                                      aTCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCgg                                                                 <  3:267220/149‑1 (MQ=255)                                                                                       tCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgc                                                                >  3:277644/1‑149 (MQ=255)                                                                                       tCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaac                                                                                                                              <  2:368695/87‑1 (MQ=37)                                                                                       tCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaac                                                                                                                              >  1:368695/1‑87 (MQ=37)                                                                                        cAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGt                                                                                           >  8:29293/1‑121 (MQ=255)                                                                                        cAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGt                                                                                           <  7:29293/121‑1 (MQ=255)                                                                                        cAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgcg                                                               >  7:227331/1‑149 (MQ=255)                                                                                             gtgtGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGcc                                                          >  8:344047/1‑149 (MQ=255)                                                                                                 ggTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCGGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAac                                                      <  8:30583/149‑1 (MQ=255)                                                                                                  gTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGcc                                                                                                           >  3:267448/1‑95 (MQ=255)                                                                                                  gTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGcc                                                                                                           <  4:267448/95‑1 (MQ=255)                                                                                                   ttGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAacac                                                    <  7:344047/149‑1 (MQ=255)                                                                                                   ttGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAaca                                                     >  4:59438/1‑148 (MQ=255)                                                                                                     ggCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACAcc                                                   <  3:49930/148‑1 (MQ=255)                                                                                                        cTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                               >  7:107060/1‑101 (MQ=255)                                                                                                                 gACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAgg                                       >  3:149082/1‑148 (MQ=255)                                                                                                                  aCGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTCCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGa                                      >  1:300806/1‑148 (MQ=255)                                                                                                                        tCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGc                               >  3:198571/1‑149 (MQ=255)                                                                                                                          cGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATCAGCCg                             >  7:61786/1‑149 (MQ=255)                                                                                                                             tACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACAGGTTGGCAGAGATAAAGGTAACATTGGAAGAGCTCAGATGTAAAAAGAAAAGCTCGTCGCGGCCGAACACCTGATCCTCAAGGAAATGAGACGTTa                          >  2:308811/1‑149 (MQ=255)                                                                                                                                 ccAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGc                      <  2:95963/149‑1 (MQ=255)                                                                                                                                  cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                   >  3:310491/1‑135 (MQ=255)                                                                                                                                  cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAAt                                   <  4:310491/135‑1 (MQ=255)                                                                                                                                  cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCa                     >  1:154821/1‑149 (MQ=255)                                                                                                                                  cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCa                     >  7:196092/1‑149 (MQ=255)                                                                                                                                  cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGATATCGCa                     >  4:352126/1‑149 (MQ=255)                                                                                                                                   aaCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCTCAAGGAAATGAGCCGTTATCGCAt                    <  8:227331/149‑1 (MQ=255)                                                                                                                                      gACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATTCGCAAGGAAATGAGCCGTTATCGCATgcc                 >  1:53223/1‑149 (MQ=255)                                                                                                                                        cAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATgccgc               <  4:149082/149‑1 (MQ=255)                                                                                                                                           gCGGCTGTTCGTCAACGGTCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCgg            <  3:233025/149‑1 (MQ=255)                                                                                                                                                  ttCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagca     <  4:123348/149‑1 (MQ=255)                                                                                                                                                     gTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg  <  8:201805/149‑1 (MQ=255)                                                                                                                                                      tCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg  >  4:17346/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GACAAGAAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTA‑CATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTG  >  NC_002947/1633397‑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: ‑