Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,819,717 (C)5→4 intergenic (‑115/+55) PA14_RS22115 ← / ← lepB DUF4845 domain‑containing protein/signal peptidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,819,7130C.98.4% 1838.8 / 14.4 386intergenic (‑111/+59)PA14_RS22115/lepBDUF4845 domain‑containing protein/signal peptidase I
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base . (218/160);  total (218/168)
Fisher's exact test for biased strand distribution p-value = 6.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

AGCAGGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTGCC  >  NC_008463/4819565‑4819858
                                                                                                                                                    |                                                                                                                                                 
 gcagGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGATGTGAAAGCGGGGCTGCCCTCGGGAGCAGccc                                                                                                                                                 <  2:869249/150‑1 (MQ=255)
 gcagGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGccc                                                                                                                                                 <  2:9667198/150‑1 (MQ=255)
 gcagGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGccc                                                                                                                                                 <  1:2106491/150‑1 (MQ=255)
  cagGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTCGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGcccc                                                                                                                                                >  2:7033716/1‑150 (MQ=255)
   agGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGccccg                                                                                                                                               >  2:5311876/1‑149 (MQ=255)
    cggACATCCCCTTCTGCGTACGCGCGTCCCTCATGTCGGTCCTCCTGTCAGTGTAGTTTTTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGccccgc                                                                                                                                              <  2:1742539/149‑3 (MQ=255)
          tCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGGTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTcgcg                                                                                                                                       >  2:10908444/1‑150 (MQ=255)
           ccccTTCTGCGAACGCGCGAACCTCATGTCGGACATCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGTATGTTAAAGCGGGGCTGCCCTCGGGAGAAGCCCCGCTTCGCGa                                                                                                                                      <  2:3620632/150‑1 (MQ=255)
              cTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGatta                                                                                                                                   >  1:10155525/1‑150 (MQ=255)
               ttCTGCGAACGCGCGAACCCCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGAAACGCCATCGGGCTGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCGGCCCTCGGGAGCAGCCCCGCTTCGCGattat                                                                                                                                  <  2:5512659/150‑1 (MQ=255)
                tCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTGCGCGattatt                                                                                                                                 >  1:1521953/1‑150 (MQ=255)
                  tGCGAACGCGCGACCCTCATGTCGGACCTCAGTTCAGTGTAGTTTTTCCGTCCAGCGAGCGGTAACGCCATCGGGCTGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTca                                                                                                                               <  2:4511371/150‑1 (MQ=255)
                  tGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTct                                                                                                                               >  2:9443210/1‑149 (MQ=255)
                  tGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCGTCGCGATTATTca                                                                                                                               >  1:1542055/1‑150 (MQ=255)
                           gcgAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAgcg                                                                                                                      <  2:7207721/150‑1 (MQ=255)
                             gAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGcc                                                                                                                    >  1:402018/1‑150 (MQ=255)
                               aacagaTGTCGCACCTAATGTCTTTGTTGTGATTCAGTCCATCGAGCGGAATCGCCATGGGTTTGGGCCCTGGTTCTCCTTTCGCGTCTGTTACAGCGGGGCGTCCCTGGGGAGCAGGGCCGCTTCGCTTTTATTCAGCAGGAGCGCCCt                                                                                                                  <  2:2523346/145‑1 (MQ=255)
                               aCCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCt                                                                                                                  >  2:8436434/1‑150 (MQ=255)
                               aCCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCt                                                                                                                  >  2:9225047/1‑150 (MQ=255)
                                  ccctggcGGACCTCTGGCCTGTGTTGTTATGCCGCCCCGCGCGCGGTAACGCCAGCGGGCCGGGCCCGGGGGCGGCGTGCGCGGATGTTACTGCGGGGCTGCCCTGGGGACCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTcc                                                                                                               <  1:1198423/144‑1 (MQ=255)
                                   cATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCcg                                                                                                              >  1:1570993/1‑150 (MQ=255)
                                      gTCGGACCTCATGTCCGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCCGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGt                                                                                                           <  2:667584/150‑1 (MQ=255)
                                      gTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGACGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGt                                                                                                           >  1:3718806/1‑150 (MQ=255)
                                         ggACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGc                                                                                                        >  2:3881896/1‑150 (MQ=255)
                                          gACCTCATGTCCGGGTAGTTATTCCGTCCAGCGAGCGGTACCGCCATCGGGCAGGGCCCTGCGGCGGCGGTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGCGCAGCCCCGCGTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCg                                                                                                       >  2:5915530/1‑150 (MQ=255)
                                           aCCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGa                                                                                                      >  2:4425648/1‑150 (MQ=255)
                                              tCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAAcc                                                                                                   >  2:10483177/1‑150 (MQ=255)
                                                 tGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTgc                                                                                                >  2:7908018/1‑150 (MQ=255)
                                                 tGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGAAACTCGATTATTCACCAGGAGCGCCCTACCACATGGCCAACCTgc                                                                                                >  2:7959645/1‑150 (MQ=255)
                                                     aGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGcac                                                                                            >  1:7781877/1‑150 (MQ=255)
                                                      gtgtAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGcaca                                                                                           >  1:10624032/1‑150 (MQ=255)
                                                      gtgtAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGcaca                                                                                           >  2:1251253/1‑150 (MQ=255)
                                                      gtgtAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGcaca                                                                                           >  1:9520569/1‑150 (MQ=255)
                                                         tAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCTGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGt                                                                                        <  1:7171052/150‑1 (MQ=255)
                                                              aTTCCGTCCAGCGAGCGTTAACGCCATCGGGCTGGGCCCGGCGGCGGCGTTCGCGGAAGTGAATGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTg                                                                                   <  1:784204/150‑1 (MQ=255)
                                                                tCCGTCCAGCTAGCGGTAACGCCATCGGGCTGGGCCGTGCGGGGGCTTTCGGGGATGTGAATGGGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGaa                                                                                 <  1:9835486/150‑1 (MQ=255)
                                                                 ccGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAAt                                                                                >  1:3671020/1‑150 (MQ=255)
                                                                  cGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATc                                                                               >  1:4788989/1‑150 (MQ=255)
                                                                  cGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATc                                                                               >  2:2309713/1‑150 (MQ=255)
                                                                   gTCCCGCGAGCGGTACCGCCATGGGGCAGGGCCCTGCGGCGGCGTTGGCGGATGTGAAAGGGGGGCTGCCCTCGGGTGCTGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCa                                                                              <  1:10553365/150‑1 (MQ=255)
                                                                   gTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCa                                                                              <  2:6458501/150‑1 (MQ=255)
                                                                    tCCAGCGAGCGGTACCGCAATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCAc                                                                             <  1:6797089/150‑1 (MQ=255)
                                                                     ccAGCGCGCGTTAACGCCTTGGGGCTGGGCCCTGGGGGGGCGTTCGCGGTAGTGAATGCGGGGCTGCCCTGGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACg                                                                            <  1:5361789/150‑1 (MQ=255)
                                                                     ccACCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGAGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACg                                                                            <  2:5158370/150‑1 (MQ=255)
                                                                      cAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGATGTGAATGGGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGc                                                                           <  2:6120178/150‑1 (MQ=255)
                                                                        gcgagcgGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCg                                                                         >  1:9066971/1‑150 (MQ=255)
                                                                          gagcgGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGAc                                                                       >  1:7624366/1‑150 (MQ=255)
                                                                           agcgGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGAcc                                                                      >  2:7815857/1‑150 (MQ=255)
                                                                            gcgGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGATGTGAATGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGAccc                                                                     <  2:3497419/150‑1 (MQ=255)
                                                                              gttAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGTAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCgg                                                                   <  1:5959345/148‑1 (MQ=255)
                                                                               gTAACGCCATCGGGCAGGGCCCTGGGGCGGCGTTCGCGTAAGTGAAAGCGGGGCTGCCCGCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCggg                                                                  <  2:7535706/150‑1 (MQ=255)
                                                                               gTAACGCCAGCGGGCAGGGCCCTGCGGCGTCGTTCGCGGAAGTGATAGCGGGGCTGCCCTCGTGAGCTGCCACGCTTCGCGATTATTCAGCAGTAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCTTTGTATCACGCCGACCgggg                                                                  <  2:10048762/150‑1 (MQ=255)
                                                                                tAACGCCATCGGGCAGGGCCCGGCGGCGGCGTTCGCGGACGTGAACGCGGGGCTGCCCGCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGga                                                                 <  1:8940351/150‑1 (MQ=255)
                                                                                   cGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGaa                                                              >  2:9561163/1‑150 (MQ=255)
                                                                                     ccATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGt                                                            >  1:11038685/1‑150 (MQ=255)
                                                                                      cagcgggccgggcCCGGCGGGGGCGTTCGCGGAAGTGACAGCGGGGCTGCCCGCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGtt                                                           <  1:6574642/147‑1 (MQ=255)
                                                                                        tCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCg                                                         >  1:2592363/1‑150 (MQ=255)
                                                                                          gggCAGGGCCCGGCGGGGGCGTTCGCGGATGTTAATGCGGGGCTGCCCGCGGGAGCACCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGc                                                       <  1:4932799/150‑1 (MQ=255)
                                                                                           ggCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCa                                                      >  2:11066769/1‑150 (MQ=255)
                                                                                            gCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAg                                                     <  2:10366173/150‑1 (MQ=255)
                                                                                             cAGGGCCCTGCGGCGGCGTTCGCGTAAGTTAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAgg                                                    <  1:5231864/150‑1 (MQ=255)
                                                                                              tgggCCCTGGGGCGGCTTTCGCGGAAGTGAAGGGGGGGCTCCCCTGGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGt                                                   <  1:8398336/149‑1 (MQ=255)
                                                                                              aGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGt                                                   >  1:9259507/1‑150 (MQ=255)
                                                                                              aGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGt                                                   >  2:5023482/1‑150 (MQ=255)
                                                                                              aGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGt                                                   <  2:4089415/150‑1 (MQ=255)
                                                                                                ggCCCTGCGGCGGCGTTCGCGGAAGTTAAAGCGGGGCTGCCCTCGGGTGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTg                                                 <  1:5162000/150‑1 (MQ=255)
                                                                                                ggCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTg                                                 >  2:10802206/1‑150 (MQ=255)
                                                                                                 gCCCTGCGGCGGCGTTCGCGGAAGTCAAAGGGGGGCTGCCCCCGGCAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGc                                                <  1:2245998/150‑1 (MQ=255)
                                                                                                  cccTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCt                                               <  2:6564334/150‑1 (MQ=255)
                                                                                                   ccTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTc                                              >  1:1289538/1‑150 (MQ=255)
                                                                                                              gTTCGCGGAAGTGAAAGCGGGGCTGACCTCGGGAGCAGCCCCGCTACGCGATTATTCAGCAGGAGCGCCATTCCGCGTGGCGAACCTGAGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTACGGCAGGTTGCTCATCTTCGGGTc                                   >  2:9086390/1‑150 (MQ=255)
                                                                                                              gTTCGCGGAAGTGAAACCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATGATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATc                                   <  2:1570221/150‑1 (MQ=255)
                                                                                                               ttCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCg                                  >  2:300208/1‑150 (MQ=255)
                                                                                                               ttCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCGTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTGCGGATCg                                  >  2:1377497/1‑150 (MQ=255)
                                                                                                                   cgGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGACTCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGcc                              >  2:2709753/1‑150 (MQ=255)
                                                                                                                   cgGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGcc                              <  2:8237313/150‑1 (MQ=255)
                                                                                                                    ggAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCa                             >  2:5774827/1‑150 (MQ=255)
                                                                                                                     gAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAg                            >  2:10577612/1‑150 (MQ=255)
                                                                                                                       aGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCt                          >  1:10352491/1‑150 (MQ=255)
                                                                                                                           aaaGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATc                      <  2:6233172/150‑1 (MQ=255)
                                                                                                                            aaGCGGGGCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATcc                     >  2:9621796/1‑150 (MQ=255)
                                                                                                                                ggggCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCacac                 <  1:6173901/150‑1 (MQ=255)
                                                                                                                                 gggCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACg                >  1:743971/1‑150 (MQ=255)
                                                                                                                                 gggCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACg                >  2:4851929/1‑150 (MQ=255)
                                                                                                                                  ggCTGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACgg               >  1:4483025/1‑150 (MQ=255)
                                                                                                                                     tGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGa            >  2:3078082/1‑150 (MQ=255)
                                                                                                                                     tGCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGa            >  1:1807086/1‑150 (MQ=255)
                                                                                                                                     ggCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGa            <  1:8757827/149‑1 (MQ=255)
                                                                                                                                      gCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTTCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGaa           <  2:3712540/150‑1 (MQ=255)
                                                                                                                                      gCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGaa           <  2:7661623/150‑1 (MQ=255)
                                                                                                                                      gCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGaa           <  1:5387032/150‑1 (MQ=255)
                                                                                                                                      gCCCTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGaa           >  1:4835388/1‑150 (MQ=255)
                                                                                                                                         cTCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGc        >  1:8246139/1‑150 (MQ=255)
                                                                                                                                          tCGGGAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGcc       <  1:3009398/150‑1 (MQ=255)
                                                                                                                                            gggAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCtt     >  1:10887319/1‑150 (MQ=255)
                                                                                                                                             ggAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTg    >  2:11185880/1‑150 (MQ=255)
                                                                                                                                              gAGCAGCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTGc   >  1:7959030/1‑150 (MQ=255)
                                                                                                                                               agcagcCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTGcc  >  2:5715198/1‑150 (MQ=255)
                                                                                                                                               agcagcCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTGcc  >  2:8839318/1‑150 (MQ=255)
                                                                                                                                               agcagcCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTGcc  >  1:5688475/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGCAGGGACATCCCCTTCTGCGAACGCGCGAACCTCATGTCGGACCTCATGTCAGTGTAGTTATTCCGTCCAGCGAGCGGTAACGCCATCGGGCAGGGCCCTGCGGCGGCGTTCGCGGAAGTGAAAGCGGGGCTGCCCTCGGGAGCAGCCCCCGCTTCGCGATTATTCAGCAGGAGCGCCCTTCCGCGTGGCGAACCTGCGCACAGGTCAGTGAATCACGCCGACCCGGGAGAAGTTCGGCAGGTTGCTCATCTTCGGATCGGGCCAGCTCATCCACACGGCGAAGGCCTTGCC  >  NC_008463/4819565‑4819858

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