Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,513,037 G→T 7.8% intergenic (+18/+44) sspB → / ← PP_1322 ClpXP protease specificity‑enhancing factor/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,513,0370GT 7.8% 811.9 / 40.5 243intergenic (+18/+44)sspB/PP_1322ClpXP protease specificity‑enhancing factor/lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (101/123);  new base T (12/7);  total (113/130)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 260 total aligned reads displayed.

CTGGAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGG  >  NC_002947/1512888‑1513186
                                                                                                                                                     |                                                                                                                                                      
  ggAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                      <  1:5554499/150‑1 (MQ=255)
  ggAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                      <  1:6187834/150‑1 (MQ=255)
    aGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggatgg                                                                                                                                                    <  2:1330556/150‑1 (MQ=255)
       cGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATc                                                                                                                                                 >  2:2831484/1‑150 (MQ=255)
       cGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATc                                                                                                                                                 <  1:476379/150‑1 (MQ=255)
        gTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCg                                                                                                                                                >  1:7273366/1‑150 (MQ=255)
        gTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAAGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCg                                                                                                                                                >  1:1035922/1‑150 (MQ=255)
           ttCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                            <  2:3230276/150‑1 (MQ=255)
           ttCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                            >  1:4225063/1‑150 (MQ=255)
            tCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCATCCGGATCGCCtt                                                                                                                                           <  1:3771277/150‑1 (MQ=255)
              gatgggtgGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCtttt                                                                                                                                         <  2:2576138/150‑1 (MQ=255)
              gatggatgGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCtttt                                                                                                                                         >  2:571391/1‑150 (MQ=255)
              gatggatgGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCtttt                                                                                                                                         >  2:6188634/1‑150 (MQ=255)
               atggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCttttt                                                                                                                                        >  2:1762533/1‑150 (MQ=255)
                tggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCATCCGGATCGCCtttttt                                                                                                                                       <  2:4860219/150‑1 (MQ=255)
                  gatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgt                                                                                                                                     <  1:1003034/150‑1 (MQ=255)
                  gatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCAATtttgt                                                                                                                                     <  2:4519725/150‑1 (MQ=255)
                    tggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                   >  1:2655026/1‑150 (MQ=255)
                    tggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                   <  1:4099958/150‑1 (MQ=255)
                    tggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                   >  2:6672630/1‑150 (MQ=255)
                     ggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgtttg                                                                                                                                  <  2:1544754/150‑1 (MQ=255)
                       cgacgaGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                <  2:7284992/150‑1 (MQ=255)
                         acgaGCTGGAAGACGATGATGTGCAGCCGGATGACGAGGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCt                                                                                                                              <  2:4924311/150‑1 (MQ=255)
                             gCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTcc                                                                                                                          >  1:3890101/1‑150 (MQ=255)
                               tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAg                                                                                                                        <  1:1330861/150‑1 (MQ=255)
                               tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAg                                                                                                                        <  2:7101743/150‑1 (MQ=255)
                               tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGCCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAg                                                                                                                        <  1:3635753/150‑1 (MQ=255)
                                 gAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGACAGTGGTCAAGTAACAAAAAAAGGCGATCCATCCGGATCGCCttt                                                                                                                                          >  2:3868108/1‑130 (MQ=255)
                                 gAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGc                                                                                                                      >  2:4068231/1‑150 (MQ=255)
                                   aGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCaa                                                                                                                    >  1:6354769/1‑150 (MQ=255)
                                     aCGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATc                                                                                                                  >  1:5811796/1‑150 (MQ=255)
                                      cGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCa                                                                                                                 <  2:5828063/150‑1 (MQ=255)
                                      cGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCa                                                                                                                 <  1:4573785/150‑1 (MQ=255)
                                      cGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCa                                                                                                                 >  2:7166047/1‑150 (MQ=255)
                                      cGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCa                                                                                                                 <  2:2704950/150‑1 (MQ=255)
                                       gatgatGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGCTCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCAc                                                                                                                >  2:7220713/1‑150 (MQ=255)
                                       gatgatGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCAc                                                                                                                >  2:6319396/1‑150 (MQ=255)
                                        atgctgTGCAGCCGGATGACGATGGCCCGCCAGACGACGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTATTTTGTTTGGCCTGTCCAGGCAATCACt                                                                                                               <  2:6124518/150‑1 (MQ=255)
                                        atgatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACt                                                                                                               <  1:5325084/150‑1 (MQ=255)
                                          gatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGa                                                                                                             <  1:1820186/150‑1 (MQ=255)
                                          gatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGa                                                                                                             <  2:5075789/150‑1 (MQ=255)
                                               gCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCtttt                                                                                                                                         <  1:978519/117‑1 (MQ=255)
                                                cAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAAt                                                                                                       >  2:2064691/1‑150 (MQ=255)
                                                 aGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGTCTGTCCAGGCAATCACTGAGTCAATc                                                                                                      >  1:7354732/1‑150 (MQ=255)
                                                 aGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATc                                                                                                      >  2:591368/1‑150 (MQ=255)
                                                 aGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATc                                                                                                      <  1:1456934/150‑1 (MQ=255)
                                                    cGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGAt                                                                                                   <  1:5101488/150‑1 (MQ=255)
                                                      gATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGt                                                                                                 >  1:1941075/1‑150 (MQ=255)
                                                       aTGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTa                                                                                                <  2:438858/150‑1 (MQ=255)
                                                            gATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTAGGTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGa                                                                                           >  2:1053785/1‑150 (MQ=255)
                                                             aTGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGaa                                                                                          <  1:3486124/150‑1 (MQ=255)
                                                                 cccGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGc                                                                                      >  1:3106139/1‑150 (MQ=255)
                                                                 cccGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGc                                                                                      >  1:2381255/1‑150 (MQ=255)
                                                                 cccGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGc                                                                                      >  2:3295919/1‑150 (MQ=255)
                                                                  ccGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCt                                                                                     >  2:1905556/1‑150 (MQ=255)
                                                                     ccAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCACCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGa                                                                                  <  2:4621249/150‑1 (MQ=255)
                                                                      cAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGAc                                                                                 >  1:3323076/1‑150 (MQ=255)
                                                                          cGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtc                                                                             <  1:2462148/150‑1 (MQ=255)
                                                                           gATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtct                                                                            <  2:1268038/150‑1 (MQ=255)
                                                                              ggAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCACTCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtcttct                                                                         >  1:4648928/1‑150 (MQ=255)
                                                                                 gggCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCa                                                                      >  2:518760/1‑150 (MQ=255)
                                                                                    cAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAATGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGc                                                                   >  1:5427589/1‑150 (MQ=255)
                                                                                              gTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACAcc                                                         <  1:1545940/150‑1 (MQ=255)
                                                                                               tCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACAccc                                                        <  1:2075377/150‑1 (MQ=255)
                                                                                                 cGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTg                                                      <  1:570760/150‑1 (MQ=255)
                                                                                                  gACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGc                                                     >  2:1813429/1‑150 (MQ=255)
                                                                                                   aCTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGca                                                    <  1:3724968/150‑1 (MQ=255)
                                                                                                        ccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACg                                               >  1:3855772/1‑150 (MQ=255)
                                                                                                        ccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGAAATCACTGAGTCAATCGATGTACTCGAAGATCTTGACGATCTTCTGCACTCGCGAAACACCCTGCACCACg                                               >  2:5012893/1‑150 (MQ=255)
                                                                                                              aaGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTg                                         >  1:5412242/1‑150 (MQ=255)
                                                                                                                   tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGa                                    >  1:2020389/1‑150 (MQ=255)
                                                                                                                   tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGa                                    >  1:1188175/1‑150 (MQ=255)
                                                                                                                   tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGa                                    >  1:4058581/1‑150 (MQ=255)
                                                                                                                   tGAAAGTGGTAAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTATTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGa                                    <  1:706454/150‑1 (MQ=255)
                                                                                                                         tGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGc                              <  2:5274499/150‑1 (MQ=255)
                                                                                                                           gTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGctt                            >  1:2317234/1‑150 (MQ=255)
                                                                                                                              aaGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttctt                         <  2:3591640/150‑1 (MQ=255)
                                                                                                                              aaGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttctt                         <  2:3014160/150‑1 (MQ=255)
                                                                                                                                   aCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg                    <  2:622820/150‑1 (MQ=255)
                                                                                                                                    cAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACggg                   <  1:6127413/150‑1 (MQ=255)
                                                                                                                                    cAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACggg                   <  2:2501199/150‑1 (MQ=255)
                                                                                                                                     caaaaaaGGCGATCCATCCGGATCGCCTTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACggg                   >  1:4211950/2‑150 (MQ=39)
                                                                                                                                     aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGt                  >  2:172464/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGt                  >  1:5950042/1‑150 (MQ=255)
                                                                                                                                        aaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGAc               >  2:1127396/1‑150 (MQ=255)
                                                                                                                                         aaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACCGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGAcc              >  2:3475639/1‑150 (MQ=255)
                                                                                                                                         aaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGAcc              >  2:6447020/1‑150 (MQ=255)
                                                                                                                                          aaGGCGCTCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCa             <  2:3932719/150‑1 (MQ=255)
                                                                                                                                           aGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAg            <  2:5953441/150‑1 (MQ=255)
                                                                                                                                           aGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAg            <  1:2206519/150‑1 (MQ=255)
                                                                                                                                               gATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGccc        >  2:907408/1‑150 (MQ=255)
                                                                                                                                                 tCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGc                                                     <  1:3162609/103‑1 (MQ=255)
                                                                                                                                                 tCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcag      <  1:4432574/150‑1 (MQ=255)
                                                                                                                                                 tCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcag      >  2:2373230/1‑150 (MQ=255)
                                                                                                                                                  ccGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagc     >  2:5054039/1‑150 (MQ=255)
                                                                                                                                                   cGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagca    <  2:5958915/150‑1 (MQ=255)
                                                                                                                                                    ggatggatCGCGTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag   >  2:3038890/1‑150 (MQ=255)
                                                                                                                                                    ggatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag   >  1:2870108/1‑150 (MQ=255)
                                                                                                                                                    ggatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag   <  1:5954313/150‑1 (MQ=255)
                                                                                                                                                     gatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAgg  <  1:6672630/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGGAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGG  >  NC_002947/1512888‑1513186

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: