Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,104,943 G→C 6.3% intergenic (‑62/‑143) PP_1875 ← / → PP_1876 two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,104,9430GC 6.3% 527.1 / 19.7 158intergenic (‑62/‑143)PP_1875/PP_1876two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (58/90);  new base C (4/6);  total (62/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 164 total aligned reads displayed.

ACCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAC  >  NC_002947/2104794‑2105091
                                                                                                                                                     |                                                                                                                                                    
aCCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCgggg                                                                                                                                                      <  1:1038680/150‑1 (MQ=255)
aCCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCgggg                                                                                                                                                      >  1:3006915/1‑150 (MQ=255)
       ttgcttgTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCtt                                                                                                                                               <  2:3441324/150‑1 (MQ=255)
        tgcttgTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTg                                                                                                                                              <  1:3827412/150‑1 (MQ=255)
         gCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGG                                                                                                                                             >  1:1095923‑M1/1‑150 (MQ=255)
          cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGC                                                                                                                                            <  1:5386101‑M1/150‑1 (MQ=255)
          cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                            <  2:4801604‑M1/150‑1 (MQ=255)
          cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                            <  2:3630588‑M1/150‑1 (MQ=255)
          cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                            <  1:1695292‑M1/150‑1 (MQ=255)
          cTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                            <  1:5313562‑M1/150‑1 (MQ=255)
           ttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCG                                                                                                                                           >  1:6829524‑M1/1‑150 (MQ=255)
                tGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCT                                                                                                                                      <  2:5680843‑M1/150‑1 (MQ=255)
                    cggCGTTTGGTTAATTTGATATCCTTGAGTCCCGTGCGGGTTTCCCGTTCCTTGCGCTGCGCAAGCTTAGCCGATTAGCCCCGGGTGGTGGTGGTTCTGTCTGTTAATTTCTGTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCcgc                                                                                                                                  <  1:3203528/149‑1 (MQ=255)
                    aGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATGGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGC                                                                                                                                  >  2:5526526‑M1/1‑150 (MQ=255)
                    aGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGC                                                                                                                                  >  2:3203528‑M1/1‑150 (MQ=255)
                        gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                              <  2:2604116‑M1/150‑1 (MQ=255)
                        gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGC                                                                                                                              >  2:6663319‑M1/1‑150 (MQ=255)
                        gATTGGTGAACTTGATATCCATGAGGCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                              <  1:3846748‑M1/150‑1 (MQ=255)
                          ttGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTT                                                                                                                            <  1:2161737‑M1/150‑1 (MQ=255)
                            ggTGACCTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGC                                                                                                                          <  2:1163001‑M1/150‑1 (MQ=255)
                             gTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCC                                                                                                                         >  1:1263806‑M1/1‑150 (MQ=255)
                             gTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCAGAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCC                                                                                                                         <  2:64458‑M1/150‑1 (MQ=255)
                              tGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCC                                                                                                                        <  2:694499‑M1/150‑1 (MQ=255)
                                 aCTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGct                                                                                                                     <  2:4517706‑M1/150‑3 (MQ=255)
                                  cTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTc                                                                                                                    <  2:2721867/150‑1 (MQ=255)
                                   ttGATATCCATGAGTCCCGAGCCGGGTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGACTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCAAGCCCGcgga                                                                                                                   <  2:5036644‑M1/150‑5 (MQ=255)
                                    ttatatccATGAGTCCCGAGCGGGGTTCCGGTTCCTTGCGCTGCGCACGCATCGCCGAATCGCCCGGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCGGCTTGGCGCCCCTCCGCGGGCAAGCCCGcggag                                                                                                                  <  1:24855‑M1/148‑6 (MQ=255)
                                    tGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTccc                                                                                                                  <  2:6239756/150‑1 (MQ=255)
                                     gATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctccca                                                                                                                 <  2:4774357‑M1/150‑7 (MQ=255)
                                     gATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCca                                                                                                                 <  2:1177191/150‑1 (MQ=255)
                                       tatCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTTCCCCTTACcaca                                                                                                               >  2:685850/1‑150 (MQ=255)
                                       tatCCAGGAGTCCCGAGCCGGTTTCCCGGTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCcaca                                                                                                               <  1:6227833/150‑1 (MQ=255)
                                        ctCCATGAGTCCCGAGCCGGTGTCCCCCTCCTCGCGCAGCGCAAGCAGAGCCGCATAGCCCCGGCCGCTGGTGGGTCTGTCTGTTCATTACAGAGGATCGCCAAGCGGCGGCCGCTGGGCGCCCCTCCGCGGGCTTGCCCGTCCCcaccg                                                                                                              <  2:1974659/149‑3 (MQ=255)
                                         tCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAgg                                                                                                             <  2:409589/150‑1 (MQ=255)
                                          ccATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCGTGCCCGGTACCACCGGt                                                                                                            >  2:6364231/1‑150 (MQ=255)
                                           cATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCGTGGCGCCCCTCCGCGGGCGGGCCCCGTCCcccacgga                                                                                                           >  2:1863141/1‑146 (MQ=255)
                                            aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTACCACAGGTAc                                                                                                          >  2:2101382/1‑150 (MQ=255)
                                            aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGCGatc                                                                                                                                        >  1:5476614‑M1/1‑117 (MQ=255)
                                            aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGCGatc                                                                                                                                        <  2:5476614‑M1/120‑4 (MQ=255)
                                             tGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACa                                                                                                         <  2:4018044/150‑1 (MQ=255)
                                                gTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCt                                                                                                      <  2:1511500/150‑1 (MQ=255)
                                                gTCCCGAGCCGGTTTCCCGTTCCTTGCGCGGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCt                                                                                                      <  2:6115396/150‑1 (MQ=255)
                                                  cccGAGCCGTTTTCCCGTTCCTTGCGCTGCGCAAGCATCGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                                    <  1:2476182/150‑1 (MQ=255)
                                                  cccGAGCCGGTTTCCCGGTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                                    <  2:4250688/150‑1 (MQ=255)
                                                    cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGGGTTCCCCGTGCCACAGGGACAGCGTAGc                                                                                                  >  1:2721723/1‑150 (MQ=255)
                                                    cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                                  <  2:1555107/150‑1 (MQ=255)
                                                        ccGGTTTCCCGTTCCTGGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCa                                                                                              <  2:905133/150‑1 (MQ=255)
                                                          ggTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTTGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAgg                                                                                            <  2:6792135/150‑1 (MQ=255)
                                                                ccGGTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                                      <  1:5078031/150‑1 (MQ=255)
                                                                  gTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTTCCCCTTCCCACAGGTACACCGTAGCCCCACGCGGAACGc                                                                                    >  2:6487131/1‑150 (MQ=255)
                                                                    tCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                                  <  1:2244911/150‑1 (MQ=255)
                                                                       ttGCGCTGCGCAACCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCGCGctcccccggggcgggcggtcagccgggggctccgctggcc                                                                               >  2:6628991‑M1/1‑110 (MQ=255)
                                                                            ctgcgcAAGCTTAGCCGAATCGCCCCGGCTGGTGGTGGGTCTGTCTTTTAATTTCAGTTGATCGCCACGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                          <  1:590780/150‑1 (MQ=39)
                                                                            ctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                          <  2:4819672/150‑1 (MQ=255)
                                                                            ctgcgcAAGCAGAGCCGAATAGCCCCGGTTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                          <  2:6017735/150‑1 (MQ=255)
                                                                             ggcgcAAGCATAGCCGAATCGCCCCGGCTGGTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTgg                                                                         <  1:3187213/149‑1 (MQ=255)
                                                                               cgcAAGCATAGCCGAATCGCCCCGGCTGCTGGTGGGTCTGTCTTTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGa                                                                       <  1:1999121/150‑1 (MQ=255)
                                                                               cgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGa                                                                       <  1:4689421/150‑1 (MQ=255)
                                                                                 catgcatAGCCGTATAGCCCCGGTTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtccctgtgggaac                                                                     <  2:3376508‑M1/147‑51 (MQ=255)
                                                                                 cacgcatAGCCGAATCGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCTGTTGATCGCCACGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaac                                                                     <  1:783913‑M1/147‑51 (MQ=255)
                                                                                 cAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTTTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGttcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaac                                                                     <  1:6034371‑M1/150‑51 (MQ=255)
                                                                                 cAAGCACTGCCGAATAGCCCCGGCTGCGGGTGGGTCTGTCTGTTCATAACAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGc                                                                     >  1:2429799‑M2/1‑149 (MQ=255)
                                                                                  aaGCATCGCCGAATCGCCCCGGCTGGTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCCGCTGGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcg                                                                    >  1:2558036‑M2/1‑148 (MQ=255)
                                                                                  aaGCATCGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAACg                                                                    <  1:1311042/150‑1 (MQ=255)
                                                                                  aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTTTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcg                                                                    >  2:2379018‑M2/1‑148 (MQ=255)
                                                                                  aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacg                                                                    <  2:4536919‑M1/150‑52 (MQ=255)
                                                                                  aaGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcg                                                                    >  2:2678489‑M2/1‑148 (MQ=255)
                                                                                   aGCATCCCCGAATCGCCCCGGCTGCTGGGGGGTGTGTCTTTTAATTTCAGTTGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                   >  1:2449023‑M2/1‑147 (MQ=255)
                                                                                   aGCATAGCCGAATCGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgg                                                                   <  1:33738‑M1/150‑53 (MQ=255)
                                                                                   aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCCGCTGGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                   >  1:4663804‑M2/1‑147 (MQ=255)
                                                                                   aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                   >  2:1263806‑M2/1‑147 (MQ=255)
                                                                                   aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                   >  2:3793799‑M2/1‑147 (MQ=255)
                                                                                   aGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                   >  1:4353967‑M2/1‑147 (MQ=255)
                                                                                    gcctAGCCGACTCGCCCCGGCTGCTGGTGGGTCTGTCTGGTCATGTCAGTGGATCGCCACGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCggg                                                                  <  2:4986641/147‑1 (MQ=39)
                                                                                    gCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGttcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacggg                                                                  <  1:5415900‑M1/150‑54 (MQ=255)
                                                                                    gCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcggg                                                                  >  2:3173349‑M2/1‑146 (MQ=255)
                                                                                     cATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgggc                                                                 <  2:294505‑M1/150‑55 (MQ=255)
                                                                                      aTAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCACGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggca                                                                >  2:6031349‑M2/1‑144 (MQ=255)
                                                                                       tAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGTGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaa                                                               >  2:6583477‑M2/1‑143 (MQ=255)
                                                                                       tAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCGTGCCCCGTCCCCCCGGGACCGCGTAGCCCCCGGCGGAACGCCCCCCCGGGGGGCGCGGGccc                                                               >  2:6588084/1‑148 (MQ=255)
                                                                                          ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCC                                                                                                                        >  1:6444618‑M1/1‑90 (MQ=255)
                                                                                          ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCC                                                                                                                        <  2:6444618‑M1/90‑1 (MQ=255)
                                                                                            gAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccg                                                          >  2:4325042‑M2/1‑138 (MQ=255)
                                                                                              atcgCCCGGGCTGCTGGTGGTTCTGTCTTTTCTTTTCTGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTGGCCCGTTCCCACAGTTACAGCTTAGCCCCAGGCTTAACGCAAGCCCTGTGGGAGCGGGCAAGCCcgcg                                                        <  1:5116532/147‑1 (MQ=255)
                                                                                                    ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCGTGCCCCGTGCCCCAGGTACAGGGTATCCCCCGGCGGAGCGCAACCCCGGGTGGAGCGGGCCAGCCCCCGGCgggg                                                  >  2:3447973/1‑150 (MQ=255)
                                                                                                         tgctgGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGGTTTCCCGGTGCCACACGTAAAGCTTATCCCCCCGCGGAACCCCAGCCCCGGGGGCGCGGGCCACCCCCCggcgcggcgccc                                             >  2:5343163/1‑148 (MQ=255)
                                                                                                         tgctgGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGCGTGCCCCGTGCCCCAGCGACACCCTAGCCCCAACCAGAACGCAACCCCCGGGGGGGCCGGGAAGCCCCCGGGGGGgcgccc                                             >  1:1669217/1‑148 (MQ=255)
                                                                                                                gggggtgtgtgtTCTTTTCTGTTGATGGCCATGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTGAGCCCCAGGCTGAACGCACGCCCTGTGGGAGCGGGCACGCCCGCGGAGGGGCGCGCAGCGGcc                                      <  2:5044854/145‑1 (MQ=38)
                                                                                                                          gtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                                   <  1:1878297‑M2/111‑4 (MQ=255)
                                                                                                                           ttCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTa                                                                                                           >  2:4380788/1‑70 (MQ=255)
                                                                                                                             cATTTCGGTGGAGCGCCACGGGGGGGCCGCGGGGCCCCCCCCCGCGGGGTTGCCCGTCCCCCCCGGCACCGCTGCGCCCCGGGCTGAAGGCACGCCCGGTGGGAGGGGGCACGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                         <  2:5582104/150‑1 (MQ=255)
                                                                                                                                   aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTgg                                                                         <  2:4053379/96‑1 (MQ=255)
                                                                                                                                   aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGC                                                                                                                          >  2:2726464‑M1/1‑47 (MQ=255)
                                                                                                                                   aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGC                                                                                                                          <  1:2726464‑M1/47‑1 (MQ=255)
                                                                                                                                    gtggATCGCAACGCGGGGGCCGCTGGGCGCCCCTCCGCGGGCTTGCCCGTGCCCACCGGTACCGCTTCGCCCCCGGCGGAACGCAAGCCCTGTGGGAGCGGGCACGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGtt                  <  2:3858465/147‑1 (MQ=37)
                                                                                                                                           gCCAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGaa                                                                                       >  1:5735871/1‑74 (MQ=37)
                                                                                                                                            ccAAGCGGCCCCCGCTTGGCGCCCCTCCGCGGGCGTGCCCCTTGCCACAGGGACAGCGGAGCCCCACGCGGGACGCCACCCCCGGGGGGGCGGGCCCCCCCCCGGCGGGGCGCCCACCGGCCCCCCAAAAACAGAGCGTGTTAGGAgtag          >  2:5232382/1‑148 (MQ=255)
                                                                                                                                             cAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCt                                                                                          >  1:1565985/1‑69 (MQ=255)
                                                                                                                                             cAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCt                                                                                          <  2:1565985/69‑1 (MQ=255)
                                                                                                                                                    gggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACCGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGTGCGGGCACGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAc  <  2:6829524/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGGGCTTGCTTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAC  >  NC_002947/2104794‑2105091

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