Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 621,820 C→T A124T (GCC→ACC)  fepD ← ferric enterobactin ABC transporter membrane subunit FebD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,8200CT99.6% 1025.9 / NA 253A124T (GCC→ACC) fepDferric enterobactin ABC transporter membrane subunit FebD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (118/134);  minor base G (1/0);  total (119/134)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

CACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGG  >  NC_000913/621722‑621917
                                                                                                  |                                                                                                 
cACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCa                                                                                                 <  2:824578/101‑1 (MQ=255)
 aCGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCAt                                                                                                >  2:824987/1‑101 (MQ=255)
  cgcccgccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATc                                                                                               <  1:824856/101‑1 (MQ=255)
     ccgccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGcc                                                                                            >  1:824991/1‑101 (MQ=255)
     ccgccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGcc                                                                                            >  2:824990/1‑101 (MQ=255)
      cgccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCa                                                                                           <  2:824533/101‑1 (MQ=255)
       gccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAg                                                                                          <  1:824626/101‑1 (MQ=255)
       gccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAg                                                                                          >  2:824994/1‑101 (MQ=255)
        ccAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAgt                                                                                         <  1:824761/101‑1 (MQ=255)
         cAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAgtt                                                                                        >  1:824995/1‑101 (MQ=255)
          aGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAgttg                                                                                       <  2:824737/101‑1 (MQ=255)
          aGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAgttg                                                                                       >  1:824998/1‑101 (MQ=255)
           gggTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAgttgt                                                                                      >  1:825000/1‑101 (MQ=255)
            ggTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAgttgtt                                                                                     <  2:824545/101‑1 (MQ=255)
             gTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTc                                                                                    <  1:824488/101‑1 (MQ=255)
               tAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCt                                                                                  >  2:825007/1‑101 (MQ=255)
               tAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCt                                                                                  >  1:825005/1‑101 (MQ=255)
               tAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCt                                                                                  <  1:824338/101‑1 (MQ=255)
                    gCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGcgcg                                                                             <  1:824895/101‑1 (MQ=255)
                      aCCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGa                                                                           <  2:824922/101‑1 (MQ=255)
                      aCCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGa                                                                           >  1:825010/1‑101 (MQ=255)
                           aCTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGt                                                                      >  2:825014/1‑101 (MQ=255)
                           aCTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGt                                                                      <  2:824644/101‑1 (MQ=255)
                           aCTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGt                                                                      <  2:824729/101‑1 (MQ=255)
                            cTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTa                                                                     >  2:825018/1‑101 (MQ=255)
                            cTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTa                                                                     >  1:825016/1‑101 (MQ=255)
                            cTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTa                                                                     >  1:825019/1‑101 (MQ=255)
                              tAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAAc                                                                   <  2:824871/101‑1 (MQ=255)
                                aCTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCa                                                                 >  1:825022/1‑101 (MQ=255)
                                 cTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCaa                                                                <  2:824431/101‑1 (MQ=255)
                                   gCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAAc                                                              <  2:824590/101‑1 (MQ=255)
                                   gCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAAc                                                              <  1:824573/101‑1 (MQ=255)
                                     ccgccgCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAg                                                            <  2:824611/101‑1 (MQ=255)
                                      cgccgcCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAgc                                                           <  1:824520/101‑1 (MQ=255)
                                      cgccgcCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAgc                                                           >  1:825027/1‑101 (MQ=255)
                                        ccgccCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAgcgc                                                         >  2:825031/1‑101 (MQ=255)
                                        ccgccCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAgcgc                                                         >  1:825029/1‑101 (MQ=255)
                                         cgccCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGcgcc                                                        >  1:825034/1‑101 (MQ=255)
                                         cgccCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGcgcc                                                        <  1:824689/101‑1 (MQ=255)
                                          gccCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGcgccg                                                       >  2:825035/1‑101 (MQ=255)
                                           cccTGACTGCCGGCAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACa                                                             >  1:825038/1‑94 (MQ=255)
                                            ccTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCa                                                     <  1:824801/101‑1 (MQ=255)
                                            ccTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCa                                                     >  2:825041/1‑101 (MQ=255)
                                              tGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAcc                                                   >  1:825043/1‑101 (MQ=255)
                                               gACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAccc                                                  <  2:824574/101‑1 (MQ=255)
                                               gACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAccc                                                  <  1:824808/101‑1 (MQ=255)
                                                aCTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCa                                                 <  1:824728/101‑1 (MQ=255)
                                                 cTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAg                                                <  1:824678/101‑1 (MQ=255)
                                                 cTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAg                                                <  2:824627/101‑1 (MQ=255)
                                                 cTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAg                                                >  1:825045/1‑101 (MQ=255)
                                                   gCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGca                                              >  2:825047/1‑101 (MQ=255)
                                                    ccGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGcac                                             <  1:824685/101‑1 (MQ=255)
                                                     cGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGcaca                                            <  2:824767/101‑1 (MQ=255)
                                                      ggTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACaa                                           <  2:824677/101‑1 (MQ=255)
                                                      ggTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACaa                                           >  2:825051/1‑101 (MQ=255)
                                                      ggTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACaa                                           >  2:825053/1‑101 (MQ=255)
                                                        tAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATg                                         <  1:824653/101‑1 (MQ=255)
                                                        tAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATg                                         >  1:825058/1‑101 (MQ=255)
                                                        tAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATg                                         >  2:825057/1‑101 (MQ=255)
                                                         aaaGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATgg                                        <  1:824630/101‑1 (MQ=255)
                                                            ggCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCaa                                     >  2:825061/1‑101 (MQ=255)
                                                             gCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCaaa                                    <  2:824659/101‑1 (MQ=255)
                                                              caacaaTCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAg                                   >  1:825063/1‑101 (MQ=255)
                                                               aacaaTCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGc                                  <  2:824743/101‑1 (MQ=255)
                                                                acaaTCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCt                                 >  1:825067/1‑101 (MQ=255)
                                                                acaaTCAATGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCt                                 <  2:824896/101‑1 (MQ=255)
                                                                   atcaatGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGc                              >  1:825069/1‑101 (MQ=255)
                                                                    tcaatGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCt                             <  2:824507/101‑1 (MQ=255)
                                                                     caatGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTc                            >  2:825071/1‑101 (MQ=255)
                                                                      aatGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTcc                           <  1:824584/101‑1 (MQ=255)
                                                                      aatGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTcc                           >  1:825073/1‑101 (MQ=255)
                                                                       atGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCg                          <  1:824618/101‑1 (MQ=255)
                                                                        tGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCgg                         <  1:824606/101‑1 (MQ=255)
                                                                        tGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCgg                         <  1:824707/101‑1 (MQ=255)
                                                                        tGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCgg                         >  1:825075/1‑101 (MQ=255)
                                                                        tGAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCgg                         <  2:824624/101‑1 (MQ=255)
                                                                         gAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGc                        <  2:824847/101‑1 (MQ=255)
                                                                         gAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGc                        <  1:824664/101‑1 (MQ=255)
                                                                         gAGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGc                        >  2:825078/1‑101 (MQ=255)
                                                                          aGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCg                       <  2:824804/101‑1 (MQ=255)
                                                                          aGGCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCg                       <  2:824705/101‑1 (MQ=255)
                                                                           ggCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGt                      >  1:825081/1‑101 (MQ=255)
                                                                            gCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGtt                     >  2:825084/1‑101 (MQ=255)
                                                                            gCCACCAGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGtt                     <  2:824698/101‑1 (MQ=255)
                                                                             ccaccaGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTc                    >  1:825085/1‑101 (MQ=255)
                                                                             ccaccaGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTc                    <  2:824758/101‑1 (MQ=255)
                                                                               accaGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCAc                  <  1:824646/101‑1 (MQ=255)
                                                                                ccaGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACg                 >  1:825088/1‑101 (MQ=255)
                                                                                ccaGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACg                 <  2:824716/101‑1 (MQ=255)
                                                                                 caGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGc                >  1:825090/1‑101 (MQ=255)
                                                                                 caGCGCCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGc                <  1:824550/101‑1 (MQ=255)
                                                                                     gcCCCGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCcaag            >  1:825094/1‑101 (MQ=255)
                                                                                       cccGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCcaagca          <  2:824770/101‑1 (MQ=255)
                                                                                       cccGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCcaagca          >  2:825097/1‑101 (MQ=255)
                                                                                        ccGGCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCcaagcaa         >  1:825098/1‑101 (MQ=255)
                                                                                          ggCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCaagcaagc       <  1:824735/101‑1 (MQ=255)
                                                                                           gCGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGcc      <  1:824781/101‑1 (MQ=255)
                                                                                            cGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCg     <  2:824682/101‑1 (MQ=255)
                                                                                            cGAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCg     >  2:825102/1‑101 (MQ=255)
                                                                                             gAAGGTCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCgg    >  1:825104/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGG  >  NC_000913/621722‑621917

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