Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 791,862 3 bp→GCT coding (1842‑1844/5430 nt) pchF → pyochelin non‑ribosomal peptide synthetase PchF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463791,8620CG99.8% 2093.2 / ‑5.8 538A614A (GCC→GCGpchFpyochelin non‑ribosomal peptide synthetase PchF
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (336/201);  total (337/201)
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
*NC_008463791,8630TC99.6% 2109.2 / ‑1.8 538W615M (TGG→CTG) pchFpyochelin non‑ribosomal peptide synthetase PchF
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (335/200);  total (337/201)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
*NC_008463791,8640GT99.8% 2283.1 / ‑4.9 538W615M (TGG→CTG) pchFpyochelin non‑ribosomal peptide synthetase PchF
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (336/201);  total (337/201)
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.84e-01

Only 100 of 538 total aligned reads displayed.

CGCGGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCCCCGAAGTGGTCGCCTGGGTCTGCCTGGAAAGCGATGTCCCGAACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTCCGG  >  NC_008463/791721‑792005
                                                                                                                                             |||                                                                                                                                             
cgcgGTGCTGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGCTCGGCATCCTCCCGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGc                                                                                                                                        <  2:4112353/151‑1 (MQ=255)
cgcgGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCGTATGTGCCGGTGGACAGCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGCAGTGGTCGCGCTGGTCTGc                                                                                                                                        >  1:2494948/1‑151 (MQ=255)
  cgGTGATGCTGCCGAAAAGCGCCGCGCACTTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCt                                                                                                                                      >  2:801537/1‑151 (MQ=255)
  cgGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCt                                                                                                                                      >  2:6736094/1‑151 (MQ=255)
  cgGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCt                                                                                                                                      <  1:9091441/151‑1 (MQ=255)
   ggTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTg                                                                                                                                     >  2:2599719/1‑151 (MQ=255)
   ggTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTg                                                                                                                                     >  2:2213050/1‑151 (MQ=255)
    gTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCTGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTgg                                                                                                                                    >  1:4597891/1‑151 (MQ=255)
      gATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGaa                                                                                                                                  <  2:5751658/151‑1 (MQ=255)
      gATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGaa                                                                                                                                  <  2:3318702/151‑1 (MQ=255)
       aTGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGaaa                                                                                                                                 >  1:10678568/1‑151 (MQ=255)
         gctgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGc                                                                                                                               >  2:8674615/1‑151 (MQ=255)
         gctgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGc                                                                                                                               <  2:229791/151‑1 (MQ=255)
          ctgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCg                                                                                                                              >  2:3697472/1‑151 (MQ=255)
           tgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGa                                                                                                                             >  2:9438440/1‑151 (MQ=255)
            gcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCTGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGAt                                                                                                                            <  1:9756023/151‑1 (MQ=255)
              cGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGt                                                                                                                          <  1:3356864/151‑1 (MQ=255)
              cGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGt                                                                                                                          >  2:7533323/1‑151 (MQ=255)
               gAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTc                                                                                                                         >  1:1660585/1‑151 (MQ=255)
                aaaaGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTcc                                                                                                                        >  2:9282558/1‑151 (MQ=255)
                 aaaGCGCCGCGCAAGTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTccc                                                                                                                       <  1:8054703/151‑1 (MQ=255)
                  aaGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGAGGTCCCg                                                                                                                      <  1:1131727/151‑1 (MQ=255)
                   aGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCgg                                                                                                                     >  1:544673/1‑150 (MQ=255)
                   aGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCgg                                                                                                                     >  2:5339841/1‑150 (MQ=255)
                        cgcgCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGcc                                                                                                                <  2:2707794/151‑1 (MQ=255)
                          cgcAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCgg                                                                                                              <  1:2945729/151‑1 (MQ=255)
                             aaTTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTg                                                                                                           >  2:9796297/1‑151 (MQ=255)
                             aaTTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTg                                                                                                           <  2:6924702/151‑1 (MQ=255)
                              aTTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTgc                                                                                                          <  1:1879329/151‑1 (MQ=255)
                                  gTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGACGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCt                                                                                                      >  1:2058508/1‑151 (MQ=255)
                                  gTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCt                                                                                                      >  1:4259630/1‑151 (MQ=255)
                                  gTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCt                                                                                                      >  1:2266968/1‑151 (MQ=255)
                                  gTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTAGGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCt                                                                                                      <  2:1267555/151‑1 (MQ=255)
                                    cgcgGTGATCGGCATCCTCCAGGCCGGCGCAGCCTAGGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGc                                                                                                    <  2:8499748/151‑1 (MQ=255)
                                        gTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTgg                                                                                                >  1:6712153/1‑151 (MQ=255)
                                        gTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTgg                                                                                                <  1:6910920/151‑1 (MQ=255)
                                        gTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTgg                                                                                                >  1:7042641/1‑151 (MQ=255)
                                          gATCGGCATCCTCCAGGCCGGCGCAGCCCATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGcc                                                                                              <  2:9094842/151‑1 (MQ=255)
                                           aTCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCa                                                                                             >  1:7227284/1‑151 (MQ=255)
                                           aTCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCa                                                                                             <  1:8477587/151‑1 (MQ=255)
                                           aTCGGCAACCTCCACGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCa                                                                                             <  2:4920343/151‑1 (MQ=255)
                                            tCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCAt                                                                                            <  1:3984835/151‑1 (MQ=255)
                                            tCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCAt                                                                                            <  2:10079883/151‑1 (MQ=255)
                                            tCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCAt                                                                                            <  1:5148634/151‑1 (MQ=255)
                                              ggCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCg                                                                                          >  1:2689052/1‑151 (MQ=255)
                                                cATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGAc                                                                                        >  2:10610158/1‑151 (MQ=255)
                                                 aTCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGAcc                                                                                       >  1:10200245/1‑151 (MQ=255)
                                                    ctccAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGc                                                                                    >  1:3572621/1‑151 (MQ=255)
                                                    ctccAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGc                                                                                    >  1:6148410/1‑151 (MQ=255)
                                                       cAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTgg                                                                                 >  1:7779779/1‑151 (MQ=255)
                                                            cGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGgccgcc                                                                            <  1:1046723/151‑1 (MQ=255)
                                                            cGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGgccgcc                                                                            >  1:6486734/1‑151 (MQ=255)
                                                            cGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGgccgcc                                                                            >  2:1920345/1‑151 (MQ=255)
                                                             ggCGCAGCCTATGTGCCGGTGGACCTCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGccgccg                                                                           <  2:5575151/151‑1 (MQ=255)
                                                              gcgcAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCGGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGa                                                                          >  1:7393140/1‑151 (MQ=255)
                                                                 cAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAg                                                                       >  1:9559087/1‑151 (MQ=255)
                                                                  aGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAgc                                                                      >  2:8192213/1‑151 (MQ=255)
                                                                      tATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCt                                                                  >  2:1810636/1‑151 (MQ=255)
                                                                       aTGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTg                                                                 <  1:4080427/151‑1 (MQ=255)
                                                                         gtgCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGgc                                                               <  2:605553/151‑1 (MQ=255)
                                                                         gtgCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGgc                                                               <  1:5943977/151‑1 (MQ=255)
                                                                             cGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGgccgcc                                                           <  2:1556857/151‑1 (MQ=255)
                                                                              ggTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGACa                                                          <  1:8817712/151‑1 (MQ=255)
                                                                               gTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAc                                                         >  2:3896400/1‑151 (MQ=255)
                                                                                tGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAcc                                                        <  1:6785608/151‑1 (MQ=255)
                                                                                 ggACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccg                                                       <  2:6597005/151‑1 (MQ=255)
                                                                                  gACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgc                                                      >  2:7647671/1‑151 (MQ=255)
                                                                                    cATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgccc                                                    >  1:7641063/1‑151 (MQ=255)
                                                                                    cATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgccc                                                    <  2:6252950/151‑1 (MQ=255)
                                                                                     aTCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgcccg                                                   >  1:558590/1‑151 (MQ=255)
                                                                                     aTCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgcccg                                                   >  1:6378839/1‑151 (MQ=255)
                                                                                       ccGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCcgcg                                                 <  2:1088587/151‑1 (MQ=255)
                                                                                        cGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCgcgg                                                >  2:2426942/1‑151 (MQ=255)
                                                                                         gCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCgcggc                                               >  2:3068342/1‑151 (MQ=255)
                                                                                         gCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCgcggc                                               >  1:6015495/1‑151 (MQ=255)
                                                                                         gCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCgcggc                                               >  2:2953137/1‑151 (MQ=255)
                                                                                                   gcTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTggcg                                     >  1:2717983/1‑151 (MQ=255)
                                                                                                     tGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTggcggc                                   <  2:10368682/151‑1 (MQ=255)
                                                                                                     tGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTggcggc                                   >  1:1378897/1‑151 (MQ=255)
                                                                                                      gcggcgCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGgcggcg                                  <  1:49907/151‑1 (MQ=255)
                                                                                                       cggcgCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGcggcgg                                 <  2:5907688/151‑1 (MQ=255)
                                                                                                            ccACGCGATCCTCGCCAGCGCCGAAGAGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGAGGCGGCGgacgac                            >  2:4610828/1‑151 (MQ=255)
                                                                                                                  gATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGcc                      <  1:2813026/151‑1 (MQ=255)
                                                                                                                  gATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGcc                      <  1:9289539/151‑1 (MQ=255)
                                                                                                                   ttCCGCGCCAGCGCGGAAGGGGTCGGGCTGGTCTGCCTGGAAAGCGATGTTCCGGTCGGCGGCTGCGGCTGCGTGGCCAGCGCCCGGCGGGGCGCCGCCAGCGCGTGGCCGGCAGAGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCt                     <  2:4472239/150‑1 (MQ=255)
                                                                                                                   aTCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCt                     >  2:893983/1‑151 (MQ=255)
                                                                                                                     ccTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTAc                   <  2:7415230/151‑1 (MQ=255)
                                                                                                                      cTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCATGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACg                  <  1:10392875/151‑1 (MQ=255)
                                                                                                                         gCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACgaga               >  2:7178598/1‑151 (MQ=255)
                                                                                                                         gCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGa               >  2:6484442/1‑151 (MQ=255)
                                                                                                                         gCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGa               >  2:10528070/1‑151 (MQ=255)
                                                                                                                          ccAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGAt              <  1:3802750/151‑1 (MQ=255)
                                                                                                                          ccAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGAt              >  1:8396328/1‑151 (MQ=255)
                                                                                                                           cAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATc             >  2:198633/1‑151 (MQ=255)
                                                                                                                               gccGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTaca         >  1:4767210/2‑151 (MQ=255)
                                                                                                                               gccGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTaca         >  1:5330011/2‑151 (MQ=255)
                                                                                                                               gccGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTaca         <  1:9938244/150‑1 (MQ=255)
                                                                                                                                   aaGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTc     >  1:10005095/1‑151 (MQ=255)
                                                                                                                                    cgTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCCGCGCCTGCGTGGCCCCCGACCGGCTGGCCGCCGACAGCGCCGGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCTCCTACGTGATCTACACCTcc    <  2:3617699/150‑1 (MQ=255)
                                                                                                                                    aGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTcc    >  2:7898720/1‑151 (MQ=255)
                                                                                                                                             |||                                                                                                                                             
CGCGGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCCCCGAAGTGGTCGCCTGGGTCTGCCTGGAAAGCGATGTCCCGAACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTCCGG  >  NC_008463/791721‑792005

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