Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,890,848 A→T 10.7% intergenic (+80/+40) PP_5164 → / ← plpB penicillin amidase family protein/NLPA lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,890,8480AT10.7% 269.0 / 21.4 92intergenic (+80/+40)PP_5164/plpBpenicillin amidase family protein/NLPA lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (37/45);  new base T (4/6);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGT  >  NC_002947/5890700‑5890994
                                                                                                                                                    |                                                                                                                                                  
cAGTATAAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                   <  2:516615/150‑1 (MQ=255)
   tATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGgtg                                                                                                                                                >  2:568972/1‑150 (MQ=255)
    aTCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGgtgt                                                                                                                                               <  2:2303253/150‑1 (MQ=255)
       aGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCgg                                                                                                                                            >  2:1006790/1‑150 (MQ=255)
         cGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCg                                                                                                                                          <  1:568972/150‑1 (MQ=255)
          gggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGt                                                                                                                                         <  2:948438/150‑1 (MQ=255)
          gggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGt                                                                                                                                         <  2:1916506/150‑1 (MQ=255)
            gcgcAGATGGGGGTGATTCCGGCGTAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                       <  2:234381/150‑1 (MQ=255)
            gcgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGttt                                                                                                                                       <  2:1453146/150‑2 (MQ=255)
            gcgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                       >  1:1359408/1‑150 (MQ=255)
                 gATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGcaca                                                                                                                                  <  2:1700463/150‑1 (MQ=255)
                  aTGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACat                                                                                                                                 >  2:288264/1‑150 (MQ=255)
                        gTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAgg                                                                                                                           >  1:1494448/1‑150 (MQ=255)
                                ggCGCAAAGTACCTCACGGCTGGGACCCCCCAAGTAAGAGGCTGTTTGGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAGCGCCGACACCAGTTGACACCGTTGGCACATATCAGGCGTTGCTa                                                                                                                   <  1:551112/150‑1 (MQ=255)
                                  cgcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                                 <  2:928225/150‑1 (MQ=255)
                                           ccTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTc                                                                                                                                              >  1:998493/1‑112 (MQ=255)
                                           ccTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTc                                                                                                                                              <  2:998493/112‑1 (MQ=255)
                                           ccTTACGGCTGCTACCCCCGCAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGCGCAGGGCTGAAGAAAAACGCCGACACAAGGGGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGa                                                                                                        <  1:1900038/150‑1 (MQ=255)
                                            cTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGc                                                                                                                                           >  1:334768/1‑114 (MQ=255)
                                            cTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGc                                                                                                                                           <  2:334768/114‑1 (MQ=255)
                                            cTTACGGATGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGaa                                                                                                       <  2:818297/150‑1 (MQ=255)
                                                  gCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                 <  2:762512/150‑1 (MQ=255)
                                                    tGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCa                                                                                               <  1:1006790/150‑1 (MQ=255)
                                                    tGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCa                                                                                               >  2:1498980/1‑150 (MQ=255)
                                                          ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCATGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACgg                                                                                         >  2:2194764/1‑150 (MQ=255)
                                                          ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                   >  2:1188336/1‑92 (MQ=255)
                                                          ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                   <  1:1188336/92‑1 (MQ=255)
                                                          ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACgg                                                                                         <  1:2115782/150‑1 (MQ=255)
                                                          ccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACgg                                                                                         >  1:1207068/1‑150 (MQ=255)
                                                          ccccGAAGTAAGAGGCTGCTTTGCGGACGCTCCGCAGCTAAACGCTGCTACGGTGCTGTGCGCGCAGGGATGAAGAAAAACGCCCCCACAAGGTGTCGGCGTTGGCACATTTCAGGCGTTGCTAGCGAAAGCAGAACGCCGGCAGTACgg                                                                                         <  2:112775/150‑1 (MQ=255)
                                                           cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGc                                                                                        >  2:2305127/1‑150 (MQ=255)
                                                           cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGc                                                                                        >  1:10203/1‑150 (MQ=255)
                                                           cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGc                                                                                        <  1:1826239/150‑1 (MQ=255)
                                                            ccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTgg                                                                                                                                      >  2:141895/1‑103 (MQ=255)
                                                            ccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTgg                                                                                                                                      <  1:141895/103‑1 (MQ=255)
                                                               aaGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCt                                                                                    <  1:1498980/150‑1 (MQ=255)
                                                                 gTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTg                                                                                  <  2:1644463/150‑1 (MQ=255)
                                                                  tAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAt                                                                                                            <  2:1502832/123‑1 (MQ=255)
                                                                  tAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAt                                                                                                            >  1:1502832/1‑123 (MQ=255)
                                                                     gagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGtttttc                                                                                                                                    <  2:361242/96‑5 (MQ=38)
                                                                     gagGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGtttttc                                                                                                                                    >  1:361242/1‑92 (MQ=38)
                                                                        gctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                   <  2:1509698/78‑1 (MQ=255)
                                                                        gctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                   >  1:1509698/1‑78 (MQ=255)
                                                                        gctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                 <  2:732913/128‑1 (MQ=255)
                                                                        gctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTActtc                                                                           >  2:378112/1‑150 (MQ=255)
                                                                        gctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATAAGGCGTTGCTAGCGAAATCAGAACGCCgg                                                                                                 >  1:732913/1‑128 (MQ=255)
                                                                            ctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTgg                                                                       <  1:2305127/150‑1 (MQ=255)
                                                                             tttGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCATGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGc                                                                      <  1:2194764/150‑1 (MQ=255)
                                                                             tttGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTActtc                                                                           <  2:1706667/145‑1 (MQ=255)
                                                                             tttGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTActtc                                                                           >  1:1706667/1‑145 (MQ=255)
                                                                                 cGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCg                                                                                                                                          >  1:961911/1‑78 (MQ=39)
                                                                                 cGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCg                                                                                                                                          <  2:961911/78‑1 (MQ=38)
                                                                                   gCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGTGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTa                                                                               <  1:1927740/135‑1 (MQ=255)
                                                                                   gCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGTGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTa                                                                               >  2:1927740/1‑135 (MQ=255)
                                                                                        ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCt                                                           >  2:2292434/1‑150 (MQ=255)
                                                                                           gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGcac                                                                                                                                   <  1:2509/75‑1 (MQ=255)
                                                                                           gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGcac                                                                                                                                   >  2:2509/1‑75 (MQ=255)
                                                                                               caAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTc                                                    <  1:378112/150‑1 (MQ=255)
                                                                                               caAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTc                                                    >  1:1139227/1‑150 (MQ=255)
                                                                                                aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTgg                                                                       >  1:231179/1‑130 (MQ=255)
                                                                                                aaaaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTgg                                                                       <  2:231179/130‑1 (MQ=255)
                                                                                                 aaaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCg                                                  <  2:1365793/150‑1 (MQ=255)
                                                                                                    ggCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTATCGAAATCAGa                                                                                                        <  2:1132443/93‑1 (MQ=255)
                                                                                                    ggCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTATCGAAATCAGa                                                                                                        >  1:1132443/1‑93 (MQ=255)
                                                                                                       ttCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTg                                                                                  >  2:800983/1‑112 (MQ=255)
                                                                                                       ttCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTg                                                                                  <  1:800983/112‑1 (MQ=255)
                                                                                                          tACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGa                                                                                                        >  2:1081425/1‑87 (MQ=255)
                                                                                                          tACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGa                                                                                                        <  1:1081425/87‑1 (MQ=255)
                                                                                                          tACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCa                                         >  2:1022128/1‑150 (MQ=255)
                                                                                                          tACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCa                                         >  2:1787030/1‑150 (MQ=255)
                                                                                                             ggTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGTTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAggg                                      >  1:661389/1‑150 (MQ=255)
                                                                                                                  tgtgCGGGCAGGGCTGAAGAAAAACGCAGACACAAAGCCAAGCCGTTGGCACATATCAGGCCTTGCTTGCGGAATCAGAAAGCCCGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCTCTGACTTCCGCGCTGGTCTggagatta                                 >  2:1560688/1‑144 (MQ=255)
                                                                                                                       gggCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAg                            >  2:2009918/1‑150 (MQ=255)
                                                                                                                          cAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTgg                                                                       <  1:2001009/104‑1 (MQ=255)
                                                                                                                          cAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTgg                                                                       >  2:2001009/1‑104 (MQ=255)
                                                                                                                             ggCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGc                                                                                                                          <  2:1731762/50‑1 (MQ=255)
                                                                                                                             ggCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGc                                                                                                                          >  1:1731762/1‑50 (MQ=255)
                                                                                                                                tGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGAt                   <  1:1022128/150‑1 (MQ=255)
                                                                                                                                tGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGAt                   <  2:1207068/150‑1 (MQ=255)
                                                                                                                                 gaagaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAGCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCGTGACTTCCGGGCTGGGCAGGGCCTTGGCCAGTTTCTGGATg                  <  1:511082/150‑1 (MQ=255)
                                                                                                                                   agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc                <  1:2292434/150‑1 (MQ=255)
                                                                                                                                    tgccaaCGCCGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaatgaat                                                               <  2:666388/98‑1 (MQ=38)
                                                                                                                                    tgccaaCGCCGACACCTTGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaatgaat                                                               >  1:666388/5‑102 (MQ=38)
                                                                                                                                            ccGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTc       <  2:1139227/150‑1 (MQ=255)
                                                                                                                                            ccGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTc       <  2:10203/150‑1 (MQ=255)
                                                                                                                                              gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCt     >  2:1742417/1‑150 (MQ=255)
                                                                                                                                              gACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCt     >  1:82666/1‑150 (MQ=255)
                                                                                                                                               acacAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTAc                                                                              <  2:1818769/76‑1 (MQ=255)
                                                                                                                                               acacAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTAc                                                                              >  1:1818769/1‑76 (MQ=255)
                                                                                                                                               acacAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCtt    <  2:82666/150‑1 (MQ=255)
                                                                                                                                                cacaAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCttg   <  2:1359408/150‑1 (MQ=255)
                                                                                                                                                 acaAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCttgt  >  1:915372/1‑150 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                                  >  2:2039654/1‑116 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                                  <  1:2039654/116‑1 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc                <  2:1769995/134‑1 (MQ=255)
                                                                                                                                                   aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc                >  1:1769995/1‑134 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGT  >  NC_002947/5890700‑5890994

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