Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,508,080 A→C 10.2% intergenic (+19/+43) 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,508,0800AC10.2% 188.1 / 11.2 67intergenic (+19/+43)sspB/PP_1322ClpXP protease specificity‑enhancing factor/lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base C (4/3);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGG  >  NC_002947/1507936‑1508228
                                                                                                                                                |                                                                                                                                                     
ccGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggatgga                                                                                                                                                   >  2:380568/1‑149 (MQ=255)
ccGTCTTCGATGGATGGCCACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggatgga                                                                                                                                                   >  7:222465/1‑149 (MQ=255)
   tCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCg                                                                                                                                                >  4:309681/1‑149 (MQ=255)
       cGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAGGCGATCCATCCGGATCGCCt                                                                                                                                            <  8:124812/147‑1 (MQ=255)
       cGATGGATGGCCACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGc                                                                                                                                               <  8:222465/146‑1 (MQ=255)
         atggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                      <  3:21600/137‑1 (MQ=255)
         atggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                      >  4:21600/1‑137 (MQ=255)
          tggatggCGACGAGCTGGAAGACGATTATGTGCAGCCGGATGACGAAGGCCCGCCACACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCttttt                                                                                                                                        <  3:241284/149‑1 (MQ=255)
          tggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCCGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCttttt                                                                                                                                        >  6:229356/1‑149 (MQ=255)
            gatggCGACGAGCTGTAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAATGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttg                                                                                                                                      <  2:234789/149‑1 (MQ=255)
            gatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCtttttt                                                                                                                                       <  3:309681/148‑1 (MQ=255)
            gacggCGACGAGCTGGAAGACGATGATGGGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGGCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCtttttt                                                                                                                                       <  8:284049/145‑1 (MQ=255)
                gCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgtttg                                                                                                                                  <  1:36775/149‑1 (MQ=255)
                  gacgaGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                <  2:105620/149‑1 (MQ=255)
                           gAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCGGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAgg                                                                                                                       <  8:39657/149‑1 (MQ=255)
                             aGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGc                                                                                                                      >  8:92957/1‑148 (MQ=255)
                                 gatgatGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCa                                                                                                                 >  8:224634/1‑149 (MQ=255)
                                  atgatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCAc                                                                                                                >  3:357018/1‑149 (MQ=255)
                                    gatgTGCAGCCGGATGACGATCGCCCGCCAGACGATGGAGGGCAGCCCCCGCGTCCGACTGGCCGGCCAAGACAGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACt                                                                                                               >  8:75365/1‑148 (MQ=255)
                                      tgtgCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTAGTTTGGCCTGTCCAGGCAATCACTGAg                                                                                                            <  3:30547/149‑1 (MQ=255)
                                         gCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAgg                                                                                                                       >  2:107451/1‑135 (MQ=255)
                                         gCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAgg                                                                                                                       <  1:107451/135‑1 (MQ=255)
                                               ggATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATc                                                                                                      <  3:169279/146‑1 (MQ=255)
                                                gATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATg                                                                                                  <  6:356405/149‑1 (MQ=255)
                                                  tGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                            <  8:176024/105‑1 (MQ=255)
                                                  tGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                            >  7:176024/1‑105 (MQ=255)
                                                           cccGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAg                                                                                       >  7:263173/1‑149 (MQ=255)
                                                              gCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCtt                                                                                    >  1:157339/1‑149 (MQ=255)
                                                                 aGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGAc                                                                                 >  5:155286/1‑149 (MQ=255)
                                                                  gACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATg                                                                                                  <  1:54557/131‑1 (MQ=255)
                                                                  gACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATg                                                                                                  >  2:54557/1‑131 (MQ=255)
                                                                          aGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTg                                                                        >  1:37474/1‑149 (MQ=255)
                                                                               aGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTcc                                                                                                                          >  7:256189/1‑94 (MQ=255)
                                                                               aGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTcc                                                                                                                          <  8:256189/94‑1 (MQ=255)
                                                                                    ccgcGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGaa                                                              <  4:184555/149‑1 (MQ=255)
                                                                                    ccgcGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGa                                                               >  3:184555/1‑148 (MQ=255)
                                                                                        gTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAacac                                                          >  6:188993/1‑149 (MQ=255)
                                                                                        gTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAacac                                                          <  6:155286/149‑1 (MQ=255)
                                                                                          ccGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACAccc                                                        >  1:273954/1‑149 (MQ=255)
                                                                                            gACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                <  3:372349/75‑1 (MQ=255)
                                                                                            gACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                >  4:372349/1‑75 (MQ=255)
                                                                                                ggccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcac                                                   >  4:88402/1‑148 (MQ=255)
                                                                                                 gccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcacca                                                 >  7:253252/1‑149 (MQ=255)
                                                                                                 gccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGAGGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcacca                                                 >  2:345145/1‑149 (MQ=255)
                                                                                                         aGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTg                                         <  8:263173/149‑1 (MQ=255)
                                                                                                          gCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTgg                                        >  1:67254/1‑149 (MQ=255)
                                                                                                             tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGcc                                      >  7:91731/1‑148 (MQ=255)
                                                                                                             tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGcc                                      >  8:183100/1‑148 (MQ=255)
                                                                                                               aaaGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAg                                   <  5:188993/149‑1 (MQ=255)
                                                                                                               aaaGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAg                                   >  2:158088/1‑149 (MQ=255)
                                                                                                                  gTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTg                                >  3:277318/1‑149 (MQ=255)
                                                                                                                      tCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGctt                            >  8:87214/1‑149 (MQ=255)
                                                                                                                       cAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttc                           >  1:322446/1‑149 (MQ=255)
                                                                                                                        aaGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttc                           <  4:277318/148‑1 (MQ=255)
                                                                                                                            aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcaccac                                                >  8:123423/5‑123 (MQ=38)
                                                                                                                            aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcaccac                                                <  7:123423/119‑1 (MQ=38)
                                                                                                                            aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGAc                      >  4:166166/5‑149 (MQ=39)
                                                                                                                            aaCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTg                                                      <  1:67044/117‑1 (MQ=255)
                                                                                                                            aaCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTg                                                      >  2:67044/1‑117 (MQ=255)
                                                                                                                             aCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACg                     <  7:92957/149‑1 (MQ=255)
                                                                                                                              cAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg                    <  3:290074/149‑1 (MQ=255)
                                                                                                                              cAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg                    <  1:158088/149‑1 (MQ=255)
                                                                                                                              cAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg                    <  3:88402/149‑1 (MQ=255)
                                                                                                                               caaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGtt                                 <  6:383369/134‑1 (MQ=255)
                                                                                                                               caaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGtt                                 >  5:383369/2‑135 (MQ=255)
                                                                                                                               aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACggg                   >  5:237554/1‑149 (MQ=255)
                                                                                                                               aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACggg                   >  4:217514/1‑149 (MQ=255)
                                                                                                                                   aaaGGCGATCCATCCGGATCGCCTTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACACTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGa                >  1:196155/1‑149 (MQ=255)
                                                                                                                                        cGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGc          >  4:372877/1‑149 (MQ=255)
                                                                                                                                         gATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGcc         <  2:157339/149‑1 (MQ=255)
                                                                                                                                         gATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATAGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGTCTTCTTGACGGGTGACCAGGcc         >  1:343030/1‑149 (MQ=255)
                                                                                                                                              ggatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCtt                                                                                    >  4:293133/1‑69 (MQ=255)
                                                                                                                                              ggatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCtt                                                                                    <  3:293133/69‑1 (MQ=255)
                                                                                                                                               gatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCTGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag   >  2:99178/1‑149 (MQ=255)
                                                                                                                                               gatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag   >  8:203334/1‑149 (MQ=255)
                                                                                                                                               gatggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag   >  8:114204/1‑149 (MQ=255)
                                                                                                                                                atggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtctt                                                                           <  4:128123/76‑1 (MQ=255)
                                                                                                                                                atggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtctt                                                                           >  3:128123/1‑76 (MQ=255)
                                                                                                                                                atggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcatg  <  2:273954/149‑3 (MQ=255)
                                                                                                                                                atggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAgg  <  8:253252/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGG  >  NC_002947/1507936‑1508228

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