Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 713,088 Δ1 bp 8.5% coding (77/531 nt) lspA → lipoprotein signal peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947713,0880T. 8.5% 209.6 / 14.5 71coding (77/531 nt)lspAlipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base T (27/36);  new base . (3/3);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GCGCCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGG  >  NC_002947/712940‑713236
                                                                                                                                                     |                                                                                                                                                     
gcgcCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggt                                                                                                                                                       >  4:57317/1‑149 (MQ=255)
 cgcCGACGTCGGCAGCCATCCGGAGCACCCGGATATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtc                                                                                                                                                      <  8:90374/149‑1 (MQ=255)
  gcCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcc                                                                                                                                                     <  4:63658/149‑1 (MQ=255)
   ccGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcct                                                                                                                                                    >  4:279643/1‑149 (MQ=255)
    cGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctg                                                                                                                                                   <  4:290825/149‑1 (MQ=255)
       cGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctggtc                                                                                                                                                <  5:72220/149‑1 (MQ=255)
           ggCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctgggcctggtcgat                                                                                                                                             >  8:14206/1‑145 (MQ=255)
           ggCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATc                                                                                                                                             <  2:378152/148‑1 (MQ=255)
                   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcc                                                                                                                                   >  6:292270/1‑149 (MQ=255)
                     cGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGGCCTGGTCGATGACCAggacca                                                                                                                                  <  7:14206/149‑5 (MQ=255)
                     cGGAGCACCCGGAAATCTGCAGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCGATGACCAggacca                                                                                                                                  <  2:77053/149‑5 (MQ=255)
                        aGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCaa                                                                                                                              <  2:135661/149‑1 (MQ=255)
                         gCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAg                                                                                                                             <  5:292270/149‑1 (MQ=255)
                         gCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAg                                                                                                                             <  3:141240/149‑1 (MQ=255)
                            cccGGAAATCTGGCGCCGTTGCGTCGACTACCTGTCCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                          <  6:168349/149‑1 (MQ=255)
                            cccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGTTGTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                          <  1:328761/149‑1 (MQ=255)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctg                                                                                                                                                   >  1:171406/1‑124 (MQ=255)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctg                                                                                                                                                   <  2:171406/124‑1 (MQ=255)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtc                                                                                                                                                      >  6:183824/1‑121 (MQ=255)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtc                                                                                                                                                      <  5:183824/121‑1 (MQ=255)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                         <  1:64005/149‑3 (MQ=38)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                          >  7:83455/1‑148 (MQ=255)
                             ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTg                                                                                                                          >  2:41708/1‑148 (MQ=255)
                               ggAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAc                                                                                                                       >  2:253426/1‑149 (MQ=255)
                               ggAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAc                                                                                                                       <  5:113764/149‑1 (MQ=255)
                                       gCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaacaa                                                                                                               >  1:77600/1‑149 (MQ=255)
                                         ggCCGTTGCGTCGACAACCTGTCCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACg                                                                                                             <  8:382797/149‑1 (MQ=255)
                                           ccGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGcc                                                                                                           <  1:253426/149‑1 (MQ=255)
                                                  gTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCGATGACCAggac                                                                                                                                    <  3:408267/118‑3 (MQ=255)
                                                  gTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCGATGACCAggac                                                                                                                                    >  4:408267/1‑116 (MQ=255)
                                                  gTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCa                                                                                                    >  3:22746/1‑149 (MQ=255)
                                                   tCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCAt                                                                                                   >  3:30817/1‑149 (MQ=255)
                                                     gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACAATg                                                                                                  >  6:126312/1‑148 (MQ=255)
                                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                               >  3:41312/1‑148 (MQ=255)
                                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                               >  2:134218/1‑148 (MQ=255)
                                                        aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                               >  7:327281/1‑148 (MQ=255)
                                                          ccTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcag                                                                                            >  2:176876/1‑149 (MQ=255)
                                                           cTGACCGGCTCGGGTGCTGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagc                                                                                           <  7:305319/149‑1 (MQ=255)
                                                           cTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAcc                                                                                                     >  2:235180/1‑139 (MQ=255)
                                                           cTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAcc                                                                                                     <  1:235180/139‑1 (MQ=255)
                                                           cTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagc                                                                                           <  1:134218/149‑1 (MQ=255)
                                                               ccGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGAt                                                                                       <  3:57317/149‑1 (MQ=255)
                                                                 ggCTCGGGTGAGGTGCGCCACTTTGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAg                                                                                                                             >  7:260589/1‑109 (MQ=255)
                                                                 ggCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAg                                                                                                                             <  8:260589/109‑1 (MQ=255)
                                                                   cTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                   <  6:382469/149‑1 (MQ=255)
                                                                   cTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                   >  2:260429/1‑149 (MQ=255)
                                                                    tCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTgg                                                                                  <  1:176876/149‑1 (MQ=255)
                                                                     cGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGt                                                                                 <  4:30817/149‑1 (MQ=255)
                                                                           aGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                            <  2:77600/148‑1 (MQ=255)
                                                                           aGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                            <  8:327281/148‑1 (MQ=255)
                                                                                    aCTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                         <  5:170681/94‑3 (MQ=18)
                                                                                    aCTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                         >  6:170681/1‑92 (MQ=18)
                                                                                            tAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCt                                                                          >  3:20459/1‑133 (MQ=255)
                                                                                            tAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCt                                                                          <  4:20459/133‑1 (MQ=255)
                                                                                             aaCCATGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAc                                                         >  3:62200/1‑149 (MQ=255)
                                                                                              aCCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAcc                                                        >  3:42268/1‑149 (MQ=255)
                                                                                              aCCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAcc                                                        <  1:41708/149‑1 (MQ=255)
                                                                                              aCCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAcc                                                        >  4:20939/1‑149 (MQ=255)
                                                                                                ccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCAt                                                                                                                                              >  6:8511/1‑62 (MQ=255)
                                                                                                ccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCAt                                                                                                                                              <  5:8511/62‑1 (MQ=255)
                                                                                                 cTGCAGGGGGGCGCTTCGGGCGCCGTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTg                                                     <  4:62200/149‑1 (MQ=255)
                                                                                                 cTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCTTGACTACTTCAGCTGGACCCTg                                                     <  7:62464/149‑1 (MQ=255)
                                                                                                 cTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTg                                                     >  8:91933/1‑149 (MQ=255)
                                                                                                    cAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCt                                                  <  3:279643/149‑1 (MQ=255)
                                                                                                         gggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACtt                                                                                                                     >  4:380822/1‑77 (MQ=255)
                                                                                                         gggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACtt                                                                                                                     <  3:380822/77‑1 (MQ=255)
                                                                                                         gggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTacaa                                             <  7:91933/149‑1 (MQ=255)
                                                                                                            cgcTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAacac                                          <  3:20939/149‑1 (MQ=255)
                                                                                                             gcTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACAcc                                         <  8:83455/149‑1 (MQ=255)
                                                                                                               ttCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGGTTACAACACCg                                        >  1:256806/1‑149 (MQ=255)
                                                                                                                 cGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGt                                      <  5:126312/148‑1 (MQ=255)
                                                                                                                   ggCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGcc                                   <  2:256806/149‑1 (MQ=255)
                                                                                                                    gcgcCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcg                                  >  3:125681/1‑149 (MQ=255)
                                                                                                                           gCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTc                             <  4:22746/147‑1 (MQ=255)
                                                                                                                                      gctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGcc                                   >  4:318918/1‑130 (MQ=255)
                                                                                                                                      gctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGcc                                   <  3:318918/130‑1 (MQ=255)
                                                                                                                                      gcggAGCTCGCTGGGCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAAGAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCt                 <  6:126696/145‑1 (MQ=255)
                                                                                                                                       ctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAc                                                                                                      <  2:345108/62‑1 (MQ=255)
                                                                                                                                         caGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACg             >  4:167003/2‑149 (MQ=38)
                                                                                                                                          aGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGtt                              <  5:15327/131‑1 (MQ=37)
                                                                                                                                          aGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGtt                              >  6:15327/1‑131 (MQ=37)
                                                                                                                                                gCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTAACAACACCGGTGCCGCGTTCAGCTTCCTGGATGACGGAGCgg       >  8:182962/1‑149 (MQ=255)
                                                                                                                                                 ctggtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCgg       <  6:398677/148‑1 (MQ=255)
                                                                                                                                                  tggtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGt     >  6:11166/1‑149 (MQ=255)
                                                                                                                                                     tcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTgg  <  4:41312/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGG  >  NC_002947/712940‑713236

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: