Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 713,091 Δ1 bp 6.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,0910T. 6.5% 193.2 / 6.1 61coding (77/531 nt)lspAlipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base T (26/31);  new base . (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGCCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGG  >  NC_002947/712943‑713234
                                                                                                                                                    |                                                                                                                                                 
gcgcCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtt                                                                                                                                                  >  1:573202/1‑149 (MQ=255)
   ccGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctg                                                                                                                                               <  1:501829/150‑1 (MQ=255)
    cGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctgg                                                                                                                                              <  1:489119/150‑1 (MQ=255)
     gACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctggt                                                                                                                                             <  1:1312295/150‑1 (MQ=255)
       cGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCa                                                                                                                                          <  2:1634224/150‑1 (MQ=255)
           ggCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGa                                                                                                                                      >  2:1059853/1‑150 (MQ=255)
                 cATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctgg                                                                                                                                              >  1:711483/1‑137 (MQ=255)
                 cATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtcctgg                                                                                                                                              <  2:711483/137‑1 (MQ=255)
                  aTCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcc                                                                                                                              >  1:860356/1‑150 (MQ=255)
                   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcca                                                                                                                             <  1:1467705/150‑1 (MQ=255)
                   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcca                                                                                                                             >  2:395243/1‑150 (MQ=255)
                   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcca                                                                                                                             >  1:1136853/1‑150 (MQ=255)
                   tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcca                                                                                                                             >  1:201662/1‑150 (MQ=255)
                       gAGCACCCGGAAATCTGCTGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGCTCCTGGTCATCGACCAGGCCACCaa                                                                                                                         <  1:122837/150‑1 (MQ=255)
                          cACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCt                                                                                                                      <  2:1348903/150‑1 (MQ=255)
                                gAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACtt                                                                                                                <  1:1059853/150‑1 (MQ=255)
                                 aaaTCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTc                                                                                                               <  2:201662/150‑1 (MQ=255)
                                   aTCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaa                                                                                                             <  2:573202/150‑1 (MQ=255)
                                     cTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaaca                                                                                                           >  1:1522970/1‑150 (MQ=255)
                                      tGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                    >  1:800519/1‑138 (MQ=37)
                                      tGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                    <  2:800519/140‑3 (MQ=37)
                                          gCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGcc                                                                                                      <  2:860356/150‑1 (MQ=255)
                                             gTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtc                                                                                                                                                  >  2:136174/1‑105 (MQ=255)
                                             gTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGctggtc                                                                                                                                                  <  1:136174/105‑1 (MQ=255)
                                                   tCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATg                                                                                             <  2:37906/150‑1 (MQ=255)
                                                       cAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCaa                                                                                                                         >  1:1498061/1‑118 (MQ=255)
                                                       cAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCaa                                                                                                                         <  2:1498061/118‑1 (MQ=255)
                                                               ccGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCaa                                                                                                                         <  2:1492350/110‑1 (MQ=255)
                                                               ccGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCaa                                                                                                                         >  1:1492350/1‑110 (MQ=255)
                                                               ccGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagca                                                                                     >  2:776698/1‑146 (MQ=255)
                                                               ccGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagca                                                                                     <  1:776698/146‑1 (MQ=255)
                                                               ccGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATc                                                                                 >  1:548539/1‑150 (MQ=255)
                                                               ccGGCGCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATc                                                                                 <  2:360644/150‑1 (MQ=255)
                                                                     cGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTc                                                                           <  2:297504/150‑1 (MQ=255)
                                                                         tGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGTACGACTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                       >  1:380212/1‑150 (MQ=255)
                                                                             gTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                   <  1:632946/150‑1 (MQ=255)
                                                                             gTGCGCAACTAAGCCCTACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                   <  2:1275322/150‑1 (MQ=255)
                                                                              tGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCGTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGac                                                                  <  2:157577/150‑1 (MQ=255)
                                                                                 gcCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAc                                                                                                 <  2:423007/116‑1 (MQ=255)
                                                                                 gcCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAc                                                                                                 >  1:423007/1‑116 (MQ=255)
                                                                                 gcCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGactac                                                               >  2:971300/1‑150 (MQ=255)
                                                                                  ccACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                    <  2:1473781/96‑3 (MQ=18)
                                                                                  ccACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                    >  1:1473781/1‑94 (MQ=18)
                                                                                       aTGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGc                                                         <  2:1522970/150‑1 (MQ=255)
                                                                                        tGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCt                                                        <  1:180591/150‑1 (MQ=255)
                                                                                        tGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCt                                                        >  2:1580893/1‑150 (MQ=255)
                                                                                          ccTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTgg                                                      >  2:231624/1‑150 (MQ=255)
                                                                                           cTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCGTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGa                                                     >  2:1098277/1‑150 (MQ=255)
                                                                                           cTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGa                                                     >  2:1552767/1‑150 (MQ=255)
                                                                                                ccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTg                                                <  1:1098277/150‑1 (MQ=255)
                                                                                                 cTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTgg                                               >  2:29507/1‑150 (MQ=255)
                                                                                                      gCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTac                                          <  2:18931/150‑1 (MQ=255)
                                                                                                      gCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAAATACCCAAAATGCTTACTATTTACGAGCATTTCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTac                                          >  1:603611/1‑150 (MQ=255)
                                                                                                           gcgcTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAacac                                     >  1:230736/1‑150 (MQ=255)
                                                                                                               ttCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGt                                 >  2:965183/1‑150 (MQ=255)
                                                                                                               ttCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGt                                 >  1:383305/1‑150 (MQ=255)
                                                                                                                       ccTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGtt                         <  2:1136853/150‑1 (MQ=255)
                                                                                                                        cTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                   >  1:854778/1‑107 (MQ=255)
                                                                                                                        cTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                   <  2:854778/107‑1 (MQ=255)
                                                                                                                         ttGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCa                       <  2:230736/150‑1 (MQ=255)
                                                                                                                         ttGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCACGTTCa                       <  1:231624/150‑1 (MQ=255)
                                                                                                                             aTGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                       >  2:425586/1‑98 (MQ=255)
                                                                                                                             aTGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                       <  1:425586/98‑1 (MQ=255)
                                                                                                                                      gctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGa          >  1:809991/1‑150 (MQ=255)
                                                                                                                                      gctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGa          <  1:354379/150‑1 (MQ=255)
                                                                                                                                            cTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGc    <  1:395243/150‑1 (MQ=255)
                                                                                                                                              tGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCgg  <  1:965183/150‑1 (MQ=255)
                                                                                                                                               gCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTc                                                            >  1:845427/1‑91 (MQ=255)
                                                                                                                                               gCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTc                                                            <  2:845427/91‑1 (MQ=255)
                                                                                                                                                   gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCt                                                        <  1:464676/91‑1 (MQ=255)
                                                                                                                                                   gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCt                                                        >  2:464676/1‑91 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGCCGACGTCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGG  >  NC_002947/712943‑713234

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