Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,100,099 C→A 6.7% intergenic (‑167/‑38) 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,100,0990CA 6.7% 380.7 / 11.2 120intergenic (‑167/‑38)PP_1875/PP_1876two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (53/59);  new base A (3/5);  total (56/64)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 125 total aligned reads displayed.

GCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTGGGGCACCACACTGGCC  >  NC_002947/2099951‑2100247
                                                                                                                                                     |                                                                                                                                                     
gctgGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGTCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGgcgcgcc                                                                                                                                                      >  1:610517/1‑149 (MQ=255)
    gtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCgg                                                                                                                                                  <  2:1606334/150‑1 (MQ=255)
      gggTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGTTACAGCTTAGCCCCAGGCTGAACGCAAGCACTGTGGGAGCGGGCAAGCCCGCGGAAGGGCGCGCAGCGGcc                                                                                                                                                >  1:1040488/1‑150 (MQ=255)
        gtctgtctGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGcccc                                                                                                                                              <  1:635751/150‑1 (MQ=255)
             tctgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaaa                                                                                                                                        >  1:1138443/1‑150 (MQ=255)
             tctgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaaa                                                                                                                                        <  2:1138443/150‑1 (MQ=255)
             tctgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCACACCCGTTCCCACAGGCACGGCTTAGCCCCAGGCTGAACGCAAGCCCTGCGGGAGCGGGCAAGCCCGCGGGGGGGCGCGCAGCGGCCCCCaaaa                                                                                                                                        >  1:1467659/1‑150 (MQ=255)
               tgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCACCaaaaaa                                                                                                                                      <  2:1195647/150‑1 (MQ=255)
                 ttCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACt                                                                                                                                    >  1:864300/1‑150 (MQ=255)
                   cATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGc                                                                                                                                                    <  2:1772168/133‑1 (MQ=255)
                   cATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGc                                                                                                                                                    >  1:1772168/1‑133 (MQ=255)
                   cATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                                                                                                                                  <  2:385802/150‑1 (MQ=255)
                   cATTCCAGTCGGTCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                                                                                                                                  <  2:1109095/150‑1 (MQ=255)
                    aTTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCAGCAGGCGGAACGCAAGCCCTGTGAGAGCGGGCAAACCCGCGGTGGGGCGCGCGGCGGCCCCCAAAAAACTgag                                                                                                                                 >  1:619376/1‑150 (MQ=255)
                      ttCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGAt                                                                                                                               >  1:1824062/1‑150 (MQ=255)
                         cgTTGCTCGCCAAGCGGGGGCCGCTTGGCCCCCCTCCGGGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGt                                                                                                                            <  1:1372358/149‑1 (MQ=255)
                         aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGt                                                                                                                            <  2:1784284/150‑1 (MQ=255)
                         aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGt                                                                                                                            >  2:906751/1‑150 (MQ=255)
                           ttGATCGCCAAGGGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTa                                                                                                                          <  2:1467659/150‑1 (MQ=255)
                               tCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAg                                                                                                                      >  2:105299/1‑150 (MQ=255)
                                   cAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAgtcgt                                                                                                                  >  1:1546172/1‑150 (MQ=255)
                                     aGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCa                                                                                                                <  2:868134/150‑1 (MQ=255)
                                     aGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGCACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCa                                                                                                                <  1:1050621/150‑1 (MQ=255)
                                         ggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGa                                                                                                            >  2:618365/1‑150 (MQ=255)
                                         ggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGa                                                                                                            <  1:105299/150‑1 (MQ=255)
                                          gggCCCCTGGGCGCCCCTCCGCGGGTTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGGGCGGGAAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAc                                                                                                           <  2:1417210/150‑1 (MQ=255)
                                           ggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGCTACAGCTTAGCCCCAGGCTGAACGCAAGCACTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAcg                                                                                                          >  1:1660965/1‑150 (MQ=255)
                                            gCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaaa                                                                                                                                        >  2:810441/1‑119 (MQ=255)
                                            gCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaaa                                                                                                                                        <  1:810441/119‑1 (MQ=255)
                                             ccGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTa                                                                                                                          <  1:788667/132‑1 (MQ=255)
                                             ccGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTa                                                                                                                          >  2:788667/1‑132 (MQ=255)
                                              cGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAAACTGGGATTGTTAGGAGTCGTCATGGAcgct                                                                                                        >  1:1427611/1‑150 (MQ=255)
                                                      gcCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATc                                                                                               >  1:615319/1‑150 (MQ=255)
                                                          cTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacagacccaccagcagccggggc                                                                                           <  2:541393‑M2/150‑49 (MQ=255)
                                                           tCCGCGGGCTTGCCCGTTCCCACAGCTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCttt                                                                                          >  2:1487843/1‑150 (MQ=255)
                                                              gcgTGCTTGCCCGTTCCCACAGGTACAGCTTACCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCt                                                                                       <  1:301058/150‑1 (MQ=255)
                                                              gcgGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCActtggcgatcaactgaaatgaacagacagacccaccagcagccggggctatt                                                                                       >  1:1578895‑M2/1‑98 (MQ=255)
                                                               cgGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTg                                                                                      <  2:864300/150‑1 (MQ=255)
                                                               cgGGCTTGCCCGTTCCCACAGGTACAGCGTAGCCCCAGGCTGAACGCAAGCACTGTGGGACAGACGTGTGCTCTTGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTg                                                                                      <  1:1790667/150‑1 (MQ=255)
                                                                  gCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCtt                                                                                                   >  2:669737/1‑134 (MQ=255)
                                                                  gCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCtt                                                                                                   <  1:669737/134‑1 (MQ=255)
                                                                     tGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacagacccaccagcagccggggctattcggctat                                                                                >  1:541393‑M2/1‑91 (MQ=255)
                                                                      gCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACcgcg                                                                               >  1:1764909/1‑150 (MQ=255)
                                                                        cgggtcccaCAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCActtggcgatcaactgaaatgaacagacagacccaccagcagccggggctattcg                                                                                     >  1:1051064‑M2/5‑88 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCActtggcgatcaactgaaatgaacagacagacccaccagcagccggggctattcg                                                                                     <  2:1051064‑M2/142‑55 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAaaaaactgaaatgaacagacagacccaccagcagccggggctattcggctatgcttgcgcag                                                                             <  1:348960‑M2/150‑63 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAcgctcgc                                                                                                     >  1:1758768/1‑126 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAcgctcgc                                                                                                     >  1:1786068/1‑126 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGc                                                                                             <  2:959143/134‑1 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGc                                                                                             >  1:959143/1‑134 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacagacccaccagcagccggggctattcggctatgct                                                                             <  2:1723994‑M2/150‑63 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacagacccacca                                                                                                      <  1:131323‑M2/125‑38 (MQ=255)
                                                                        ccGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacagacccacca                                                                                                      >  2:131323‑M2/1‑88 (MQ=255)
                                                                          gTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacaga                                                                                                             <  1:1632925‑M2/116‑31 (MQ=255)
                                                                          gTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCCAAGCGGCCCCCGcttggcgatcaactgaaatgaacagacaga                                                                                                             >  2:1632925‑M2/1‑86 (MQ=255)
                                                                            tCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAAt                                                                         <  2:1764909/150‑1 (MQ=255)
                                                                              ccACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAc                                                                                                           >  2:595750/1‑114 (MQ=255)
                                                                              ccACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAc                                                                                                           <  1:595750/114‑1 (MQ=255)
                                                                              ccACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCg                                                                       <  1:971313/150‑1 (MQ=255)
                                                                              ccACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCg                                                                       >  2:1126390/1‑150 (MQ=255)
                                                                              ccACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCg                                                                       <  1:1145031/150‑1 (MQ=255)
                                                                                  aGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTAc                                                                   <  1:1126390/150‑1 (MQ=255)
                                                                                    gTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                                                                                                                                  >  2:45146/1‑85 (MQ=255)
                                                                                    gTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAg                                                                                                                         >  2:44268/1‑94 (MQ=255)
                                                                                    gTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAg                                                                                                                         <  1:44268/94‑1 (MQ=255)
                                                                                     tACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTgg                                                                <  1:1344678/150‑1 (MQ=255)
                                                                                      aCAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGgc                                                               <  2:657010/150‑1 (MQ=255)
                                                                                          cTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGgcgcgc                                                           <  2:1540669/150‑1 (MQ=255)
                                                                                          cTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGgcgcgc                                                           >  1:1542901/1‑150 (MQ=255)
                                                                                                 ccAGGCGGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGGGCCGAATCGGTACTGGCGCGCCTTGAAc                                                    <  1:1809378/150‑1 (MQ=255)
                                                                                                   aGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAAccc                                                  <  1:1244385/150‑1 (MQ=255)
                                                                                                     gCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATgg                                                                                                             <  2:1162887/89‑1 (MQ=255)
                                                                                                     gCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCt                                                <  1:1894982/150‑1 (MQ=255)
                                                                                                        gAACGCACGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCCCGCAGCGGCCCAACACTACCTATGTCGGCGAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGtt                                             <  1:442144/150‑1 (MQ=255)
                                                                                                        gAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGtt                                             <  2:1824062/150‑1 (MQ=255)
                                                                                                         aaCGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTg                                            >  1:716513/1‑150 (MQ=255)
                                                                                                          aCGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCGTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGc                                           >  2:840628/1‑150 (MQ=255)
                                                                                                           cGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAaaaaactgaaatgaacagacagacccaccagcagccggggctattcggctatgcttgcgcagcgcaaggaacgggaaaccggctcgggactcatgga                                          >  2:348960‑M2/1‑53 (MQ=255)
                                                                                                           cGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGcc                                          <  2:1542901/150‑1 (MQ=255)
                                                                                                            gCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATc                                                                                               >  2:1105348/1‑96 (MQ=255)
                                                                                                            gCAAGCCCTGTGGGAGCGGGCAAGCACGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATc                                                                                               <  1:1105348/96‑1 (MQ=255)
                                                                                                                 cccTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCttt                                                                                          <  2:1780024/96‑1 (MQ=255)
                                                                                                                 cccTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGcaca                                    >  2:760302/1‑150 (MQ=255)
                                                                                                                 cccTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGcaca                                    >  1:845004/1‑150 (MQ=255)
                                                                                                                  ccTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGa                                                                                                                       <  2:1407372/66‑1 (MQ=255)
                                                                                                                  ccTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGa                                                                                                                       >  1:1407372/1‑66 (MQ=255)
                                                                                                                   cTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAAc                                  >  2:1521376/1‑150 (MQ=255)
                                                                                                                           gCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTAGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCTCCTTGAACCACTGTTGCCGGCATCACACCCGcac                          >  1:401094/1‑150 (MQ=255)
                                                                                                                              ggCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGcc                                          <  1:299555/131‑1 (MQ=255)
                                                                                                                              ggCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGcc                                          >  2:299555/1‑131 (MQ=255)
                                                                                                                                  aGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACt                   <  1:1487843/150‑1 (MQ=255)
                                                                                                                                   gCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTg                  >  2:5894/1‑150 (MQ=255)
                                                                                                                                     ccGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTggg                <  2:1660965/150‑1 (MQ=255)
                                                                                                                                      cgcgGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTgggg               >  1:1941210/1‑150 (MQ=255)
                                                                                                                                          ggAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTGGGGcac            <  2:1546172/150‑1 (MQ=255)
                                                                                                                                            aGGGGCGCCAAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATAGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTGGGGcacca          <  2:126023/150‑1 (MQ=255)
                                                                                                                                                  gcgcAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTGGGGCACCACACTgg    <  2:615319/150‑1 (MQ=255)
                                                                                                                                                   cgcAGCGGCCCCCActtggcgatcaactgaaatgaacagacagacccaccagcagccggggctattcggctatgcttgcgcagcgcaaggaacgggaaaccggctcgggactcatggatatcaagttcaccaatcgcctgtcatacaagc   >  2:1849543‑M2/1‑14 (MQ=255)
                                                                                                                                                   cgcAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTGGGGCACCACACTGGc   >  2:1796101/1‑150 (MQ=255)
                                                                                                                                                    gcAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTATGGACCGCGCCGAATCGGTACTAGCGCGCCTTGAAACCGTGTTGCCGCCACACCGCCCGCACATCAACTGGGACCCCCCACTGGcc  >  2:545349/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAACCCCTGTTGCCGGCACAACGCCCGCACATCGACTGGGGCACCACACTGGCC  >  NC_002947/2099951‑2100247

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