Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,508,090 T→C 10.7% intergenic (+20/+42) 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,0900TC10.7% 223.0 / 12.9 84intergenic (+20/+42)sspB/PP_1322ClpXP protease specificity‑enhancing factor/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (4/5);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GGAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGGT  >  NC_002947/1507941‑1508238
                                                                                                                                                     |                                                                                                                                                     
ggAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                       <  1:1431473/150‑1 (MQ=255)
ggAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                       <  2:598275/150‑1 (MQ=255)
ggAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTAATCAATTAACAAAAAAAGGCGATCCggat                                                                                                                                                       <  2:692081/150‑1 (MQ=255)
     cGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATc                                                                                                                                                  <  2:2235987/150‑1 (MQ=255)
         ttCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                             >  1:986383/1‑150 (MQ=255)
                  tggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                    >  1:1321089/1‑150 (MQ=255)
                  gggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                    <  1:1299097/149‑1 (MQ=255)
                   ggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCATCCGGATCGCCtttttttgtta                                                                                                                                   <  2:1312289/150‑6 (MQ=39)
                    gCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCATCCGGATCGCCtttttttgttac                                                                                                                                  <  2:1034917/150‑7 (MQ=38)
                     cgacgaGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                 <  1:607238/150‑1 (MQ=255)
                     cgacgaGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGAAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                 <  2:2193909/150‑1 (MQ=255)
                             tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAg                                                                                                                         <  1:217724/150‑1 (MQ=255)
                                  gACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCGTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTCTGTTGTGGCCTGTCCAGGCAAt                                                                                                                    >  1:121175/1‑150 (MQ=255)
                                   aCGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                       >  2:36884/1‑115 (MQ=255)
                                   aCGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggat                                                                                                                                                       <  1:36884/115‑1 (MQ=255)
                                   aCGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATc                                                                                                                   >  1:1127157/1‑150 (MQ=255)
                                      atgatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                             <  2:1050413/121‑1 (MQ=255)
                                      atgatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCt                                                                                                                                             >  1:1050413/1‑121 (MQ=255)
                                        gatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGa                                                                                                              <  2:986383/150‑1 (MQ=255)
                                        gatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGACCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGa                                                                                                              <  1:2401221/150‑1 (MQ=255)
                                           gtgCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTc                                                                                                           >  1:1327773/1‑150 (MQ=255)
                                               aGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATc                                                                                                       <  2:2269456/150‑1 (MQ=255)
                                                gCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCg                                                                                                      >  2:728608/1‑150 (MQ=255)
                                                  cgaaTGACAATGGCCCGCCAGAGGAGGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAACAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGAt                                                                                                    <  2:1321089/147‑1 (MQ=255)
                                                  cGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGAt                                                                                                    <  1:1789276/150‑1 (MQ=255)
                                                    gATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGt                                                                                                  >  1:1366667/1‑150 (MQ=255)
                                                         cGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCg                                                                                             >  2:248609/1‑150 (MQ=255)
                                                          gATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                 >  1:735077/1‑113 (MQ=255)
                                                          gATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                 <  2:735077/113‑1 (MQ=255)
                                                          gATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGa                                                                                            >  1:1405117/1‑150 (MQ=255)
                                                            ttgCCCGCCAGACGATGGAGGGCAGCCGCCGGGTCCGCCTGGCCGGCGAAGCCGGAAAGTTGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAAc                                                                                          <  2:390539/148‑1 (MQ=255)
                                                                ccGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCt                                                                                      <  1:728608/150‑1 (MQ=255)
                                                                 cGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCtt                                                                                     <  1:1939949/150‑1 (MQ=255)
                                                                 cGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCtt                                                                                     >  1:1937958/1‑150 (MQ=255)
                                                                  gCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCCGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTg                                                                                    >  1:258164/1‑150 (MQ=255)
                                                                   ccAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGa                                                                                   <  1:977211/150‑1 (MQ=255)
                                                                        cGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtc                                                                              <  1:1346559/150‑1 (MQ=255)
                                                                          aTGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtctt                                                                            <  1:904908/150‑1 (MQ=255)
                                                                            ggAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGaa                                                                                           <  1:1287235/133‑1 (MQ=255)
                                                                            ggAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGaa                                                                                           >  2:1287235/1‑133 (MQ=255)
                                                                            ggAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAGGCGATCCCTCCGGATCGCCttttttt                                                                                                                                       <  2:1919802/88‑2 (MQ=25)
                                                                            ggAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAGGCGATCCATCCGGATCGCCttttttt                                                                                                                                       >  1:1919802/1‑87 (MQ=21)
                                                                            ggAGGCCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGCTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtaccag                                                                          <  1:1264467/150‑6 (MQ=255)
                                                                             gAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTg                                                                         <  1:1156449/150‑1 (MQ=255)
                                                                               gggCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                 >  2:2110143/1‑92 (MQ=255)
                                                                               gggCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGc                                                                                                                                 <  1:2110143/92‑1 (MQ=255)
                                                                                ggCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCAc                                                                      <  2:342243/150‑1 (MQ=255)
                                                                                   aGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                    <  2:2244432/85‑1 (MQ=255)
                                                                                   aGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                    >  1:2244432/1‑85 (MQ=255)
                                                                                     ccgccgCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                    >  1:1962908/1‑83 (MQ=255)
                                                                                     ccgccgCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTtttgttt                                                                                                                                    <  2:1962908/83‑1 (MQ=255)
                                                                                             tCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACAccc                                                         <  1:248609/150‑1 (MQ=255)
                                                                                                gACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGc                                                      <  2:1366667/150‑1 (MQ=255)
                                                                                                    ggccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCACGCAATCACTGAGTCAATCGAGGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcacca                                                  <  2:1327773/150‑1 (MQ=255)
                                                                                                      ccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACACCCGAAACACCCTGCACCACg                                                >  1:430429/1‑150 (MQ=255)
                                                                                                       cggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGa                                                                                   >  1:1474285/1‑114 (MQ=255)
                                                                                                       cggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGa                                                                                   <  2:1474285/114‑1 (MQ=255)
                                                                                                           cAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGt                                           >  1:1989839/1‑150 (MQ=255)
                                                                                                             aGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTgg                                         <  2:1127157/150‑1 (MQ=255)
                                                                                                             aGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACCGCAGTgg                                         >  2:637246/1‑150 (MQ=255)
                                                                                                              gCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGc                                                                                       >  1:2202421/1‑103 (MQ=255)
                                                                                                              gCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGc                                                                                       <  2:2202421/103‑1 (MQ=255)
                                                                                                              gCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGc                                        >  2:399200/1‑150 (MQ=255)
                                                                                                                cTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCg                                      >  2:1231190/1‑150 (MQ=255)
                                                                                                                 tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGa                                     <  2:1695770/150‑1 (MQ=255)
                                                                                                                 tGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGa                                     >  1:590854/1‑150 (MQ=255)
                                                                                                                   aaaGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGt                                   <  2:1989839/150‑1 (MQ=255)
                                                                                                                    aaGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGtt                                  >  1:199477/1‑150 (MQ=255)
                                                                                                                         gTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGtt                                  <  2:2081646/145‑1 (MQ=255)
                                                                                                                         gTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGtt                                  >  1:2081646/1‑145 (MQ=255)
                                                                                                                            aaGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACGACGGCAGTGGCCGAGGTGGcttctt                          >  1:443065/1‑150 (MQ=255)
                                                                                                                             gccaaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTg                         >  2:1303837/8‑150 (MQ=38)
                                                                                                                              ccaaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTg                         >  1:1459409/7‑150 (MQ=38)
                                                                                                                               caaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGAc                       <  1:1303837/145‑1 (MQ=39)
                                                                                                                                aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttctt                          >  1:470698/5‑147 (MQ=38)
                                                                                                                                aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttctt                          <  2:470698/143‑1 (MQ=38)
                                                                                                                                 aCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtc                                                                              <  2:1462218/93‑1 (MQ=255)
                                                                                                                                 aCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtc                                                                              >  1:1462218/1‑93 (MQ=255)
                                                                                                                                 aCAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg                     <  2:1405117/150‑1 (MQ=255)
                                                                                                                                  cAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACggg                    >  1:1634159/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTgg                                         <  2:1216859/128‑1 (MQ=255)
                                                                                                                                   aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTgg                                         >  1:1216859/1‑128 (MQ=255)
                                                                                                                                   aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgggg                   >  1:2259137/1‑149 (MQ=255)
                                                                                                                                   aaaaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTATCCAGGCAAACACTGGCTCAATCGATGTGCGCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGt                   >  1:887535/1‑150 (MQ=255)
                                                                                                                                     aaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTgg                                <  2:2446463/135‑1 (MQ=255)
                                                                                                                                     aaaaaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTgg                                >  1:2446463/1‑135 (MQ=255)
                                                                                                                                        aaGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCa              <  1:923297/150‑1 (MQ=255)
                                                                                                                                              aTCCGGATGGATCGCCTTTTTTTTTTTTTTTTTTCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTTGCTTCTTGACGGGTGACCAGGCCca        >  2:643878/1‑150 (MQ=255)
                                                                                                                                              aTCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCca        <  2:887535/150‑1 (MQ=255)
                                                                                                                                                    atggatCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGGt  <  1:399200/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAGCCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCCTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAGGT  >  NC_002947/1507941‑1508238

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