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.1% 2047.2 / 1.5 534A614A (GCC→GCGpchFpyochelin non‑ribosomal peptide synthetase PchF
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base G (338/190);  total (344/190)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*NC_008463791,8630TC99.1% 2046.7 / 4.0 534W615M (TGG→CTG) pchFpyochelin non‑ribosomal peptide synthetase PchF
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base C (339/190);  total (344/190)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_008463791,8640GT98.9% 2200.0 / ‑3.2 534W615M (TGG→CTG) pchFpyochelin non‑ribosomal peptide synthetase PchF
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base T (338/190);  total (344/190)
Fisher's exact test for biased strand distribution p-value = 9.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 534 total aligned reads displayed.

TCGCGGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCCCCGAAGTGGTCGCCTGGGTCTGCCTGGAAAGCGATGTCCCGAACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTCCGG  >  NC_008463/791720‑792005
                                                                                                                                              |||                                                                                                                                             
   cgGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTCTGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCt                                                                                                                                      <  2:7674381/151‑1 (MQ=255)
   cgGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCt                                                                                                                                      <  1:4992889/151‑1 (MQ=255)
   cgGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCt                                                                                                                                      >  1:2517798/1‑151 (MQ=255)
    ggTGATGCTGCCGAAACGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTg                                                                                                                                     >  2:6755077/1‑151 (MQ=255)
      tGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGa                                                                                                                                   <  1:1180676/151‑1 (MQ=255)
      tGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGa                                                                                                                                   >  1:9356885/1‑151 (MQ=255)
      tGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGa                                                                                                                                   <  2:1561647/151‑1 (MQ=255)
      agATGCTGCCGACAAGCGCCGCGCAATTGGACGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTAGGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGa                                                                                                                                   <  1:10030177/150‑1 (MQ=255)
       gATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGaa                                                                                                                                  >  1:187051/1‑151 (MQ=255)
       gATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGaa                                                                                                                                  >  1:9105783/1‑151 (MQ=255)
        aTGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGaaa                                                                                                                                 <  2:3688085/151‑1 (MQ=255)
          gctgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGc                                                                                                                               >  2:1639762/1‑151 (MQ=255)
           ctgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCg                                                                                                                              >  1:3082343/1‑151 (MQ=255)
            tgcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGa                                                                                                                             >  1:519017/1‑151 (MQ=255)
             gcCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGAt                                                                                                                            <  2:4681213/151‑1 (MQ=255)
                gAAAAGCGCCGCGCAATTGGTCGAGGTGATCGGCATCCTCAAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTAGCCAGCGCCGAAGTGGTCGCGCTGGACTGCCTGGAAAGCGATGTc                                                                                                                         >  2:7233028/1‑151 (MQ=255)
                 aaaaGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTcc                                                                                                                        >  1:3328812/1‑151 (MQ=255)
                  aaaGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTccc                                                                                                                       <  1:2730620/151‑1 (MQ=255)
                   aaGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCg                                                                                                                      <  1:3722917/151‑1 (MQ=255)
                    aGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGTAAGCGATGTcgcgg                                                                                                                     >  1:8245575/1‑150 (MQ=255)
                     gcgcCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCgga                                                                                                                    >  2:6266154/1‑149 (MQ=255)
                     gcgcCGCGCAATTGCTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGAACCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCgga                                                                                                                    <  1:1022323/151‑3 (MQ=255)
                         cgcgCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGcc                                                                                                                >  1:4161155/1‑151 (MQ=255)
                           cgcAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCgg                                                                                                              >  1:3602861/1‑151 (MQ=255)
                           cgcAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCgg                                                                                                              >  2:5425645/1‑151 (MQ=255)
                                ttGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTgcg                                                                                                         <  1:7449061/151‑1 (MQ=255)
                                   gTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCt                                                                                                      >  1:2018392/1‑151 (MQ=255)
                                    tCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTg                                                                                                     <  2:5680901/151‑1 (MQ=255)
                                        ggTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTg                                                                                                 >  1:8882433/1‑151 (MQ=255)
                                        ggTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTg                                                                                                 <  1:7783149/151‑1 (MQ=255)
                                          tGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCTGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGCAGTGGTCGCGCTGGTCTGCCTGGTAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGc                                                                                               >  2:8084322/1‑151 (MQ=255)
                                          tGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGc                                                                                               >  2:7902505/1‑151 (MQ=255)
                                           gATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACAGCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGcc                                                                                              <  1:5647409/151‑1 (MQ=255)
                                            aTCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCa                                                                                             <  2:3552614/151‑1 (MQ=255)
                                            aTCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCa                                                                                             >  2:7001964/1‑151 (MQ=255)
                                              cGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATc                                                                                           <  1:3460206/151‑1 (MQ=255)
                                               ggCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCg                                                                                          >  1:9946005/1‑151 (MQ=255)
                                                 cATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGAc                                                                                        <  2:8480800/151‑1 (MQ=255)
                                                 cATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGAc                                                                                        <  2:2924158/151‑1 (MQ=255)
                                                 cATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGAc                                                                                        >  1:2123342/1‑151 (MQ=255)
                                                  aTCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGAcc                                                                                       >  2:8022351/1‑151 (MQ=255)
                                                     ctccAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGc                                                                                    <  2:234222/151‑1 (MQ=255)
                                                      tccAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCt                                                                                   <  1:287928/151‑1 (MQ=255)
                                                       ccAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTg                                                                                  <  2:245385/151‑1 (MQ=255)
                                                                cgcAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGAc                                                                         >  2:9734923/1‑151 (MQ=255)
                                                                cgcAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGAc                                                                         >  1:7843316/1‑151 (MQ=255)
                                                                 gcAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACa                                                                        <  2:2668153/151‑1 (MQ=255)
                                                                  cAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCCTCGACCGGCTGGCCGCCGACAg                                                                       >  1:5088672/1‑151 (MQ=255)
                                                                   aGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGAGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAgc                                                                      >  1:9036151/1‑151 (MQ=255)
                                                                    gCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAgcg                                                                     >  2:5096979/1‑151 (MQ=255)
                                                                     ccTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAgcgc                                                                    >  1:7234410/1‑151 (MQ=255)
                                                                      cTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGcc                                                                   >  2:7822028/1‑151 (MQ=255)
                                                                       tATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCt                                                                  <  1:7101682/151‑1 (MQ=255)
                                                                          gtgCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGgc                                                               <  1:9862244/151‑1 (MQ=255)
                                                                           tgCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGgcc                                                              <  1:11112319/151‑1 (MQ=255)
                                                                                gTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAc                                                         >  1:9269435/1‑151 (MQ=255)
                                                                                 tGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAcc                                                        <  1:6361077/151‑1 (MQ=255)
                                                                                   gACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgc                                                      >  2:9409428/1‑151 (MQ=255)
                                                                                   gACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgc                                                      >  1:8574762/1‑151 (MQ=255)
                                                                                    aCATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgcc                                                     >  2:2001701/1‑151 (MQ=255)
                                                                                      aTCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgcccg                                                   >  2:10746085/1‑151 (MQ=255)
                                                                                      aTCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCAccgcccg                                                   <  1:9602515/151‑1 (MQ=255)
                                                                                        ccGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCGGCGTGGCCATCGACCGGCTGGCCGCCGACTGCGCCTGGCCGCCACCTCCagcg                                                 >  2:10414563/1‑151 (MQ=255)
                                                                                           ccAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCgcggcg                                              <  1:4767809/151‑1 (MQ=255)
                                                                                            cAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGcggcgg                                             >  2:3286638/1‑151 (MQ=255)
                                                                                             aGCCGCCGCAGCGGCGCCCGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGa                                            <  1:5764754/151‑1 (MQ=255)
                                                                                              gccgccGCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAg                                           <  2:64897/151‑1 (MQ=255)
                                                                                               ccgccgCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAgg                                          >  1:10349574/1‑151 (MQ=255)
                                                                                               ccgccgCTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAgg                                          <  2:10980293/151‑1 (MQ=255)
                                                                                                cgccgcTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGt                                         <  2:5991997/151‑1 (MQ=255)
                                                                                                  ccgcTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTgg                                       >  2:7537292/1‑151 (MQ=255)
                                                                                                  ccgcTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTgg                                       >  2:4836565/1‑151 (MQ=255)
                                                                                                  ccgcTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTgg                                       >  1:5926945/1‑151 (MQ=255)
                                                                                                    gcTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTggcg                                     >  2:9029803/1‑151 (MQ=255)
                                                                                                     cTGCGGCGCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTggcgg                                    <  2:5887247/151‑1 (MQ=255)
                                                                                                      tGCGGCGCCAGGCGATCCTCGCCAGCGCCGACGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCGTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGGggcggc                                   >  1:8281047/1‑151 (MQ=255)
                                                                                                       gcggcgCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGgcggcg                                  <  1:264159/151‑1 (MQ=255)
                                                                                                         ggcgCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGga                                <  2:5604883/151‑1 (MQ=255)
                                                                                                         ggcgCCAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGga                                >  2:9069743/1‑151 (MQ=255)
                                                                                                              cAGGCGATCCTCGCCCGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACCGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGAcc                           >  2:314584/1‑151 (MQ=255)
                                                                                                              cAGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGAcc                           >  2:1226010/1‑151 (MQ=255)
                                                                                                               aGGCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCt                          >  1:4910196/1‑151 (MQ=255)
                                                                                                                cgCGATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCAGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTc                         <  1:8001206/150‑1 (MQ=255)
                                                                                                                   gATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGACAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCCTCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGcg                      >  2:9591681/1‑150 (MQ=255)
                                                                                                                   gATCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGcc                      >  1:6957230/1‑151 (MQ=255)
                                                                                                                     tCCTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTa                    >  2:10229229/1‑151 (MQ=255)
                                                                                                                      ccTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTAc                   >  2:7497294/1‑151 (MQ=255)
                                                                                                                      ccTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTAc                   <  2:1887529/151‑1 (MQ=255)
                                                                                                                       cTCGCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGTCAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACg                  >  2:10554122/1‑151 (MQ=255)
                                                                                                                          gCCAGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACCGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGa               >  2:8791829/1‑151 (MQ=255)
                                                                                                                             agcgacGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGAGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGAtat            <  2:2123203/146‑3 (MQ=255)
                                                                                                                             aGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCt            >  1:10744808/1‑151 (MQ=255)
                                                                                                                             aGCGCCGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCt            >  1:6157951/1‑151 (MQ=255)
                                                                                                                               cgccGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTac          >  2:9156168/3‑151 (MQ=255)
                                                                                                                                gccGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTaca         >  2:8290688/2‑151 (MQ=255)
                                                                                                                                  cGAAGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACAcc       >  1:2291840/1‑151 (MQ=255)
                                                                                                                                   gacgTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCt      >  1:3941931/4‑151 (MQ=255)
                                                                                                                                   gAAGTGGTCGCGCTGGTCCGCCTGGAAAGCGATGTCCCGGCCGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCCCGCCTACGTGATCTACACCt      <  2:1014187/151‑1 (MQ=255)
                                                                                                                                     aGTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTcc    >  1:1363086/1‑151 (MQ=255)
                                                                                                                                      gTGGTCGCGCTGGTCTGCCTGGAAAGCGATGTCCCGGACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTCCg   >  2:4039504/1‑151 (MQ=255)
                                                                                                                                              |||                                                                                                                                             
TCGCGGTGATGCTGCCGAAAAGCGCCGCGCAATTGGTCGCGGTGATCGGCATCCTCCAGGCCGGCGCAGCCTATGTGCCGGTGGACATCCGCCAGCCGCCGCTGCGGCGCCAGGCGATCCTCGCCAGCCCCGAAGTGGTCGCCTGGGTCTGCCTGGAAAGCGATGTCCCGAACGCCGGCTGCGCCTGCGTGGCCATCGACCGGCTGGCCGCCGACAGCGCCTGGCCGCCACCGCCCGCGGCGGAGGTGGCGGCGGACGACCTCGCCTACGTGATCTACACCTCCGG  >  NC_008463/791720‑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: