Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,095,235 G→C 12.5% intergenic (‑243/+74) htpX ← / ← alaA protease HtpX/glutamate‑pyruvate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,095,2350GC12.5% 233.2 / 21.7 78intergenic (‑243/+74)htpX/alaAprotease HtpX/glutamate‑pyruvate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (37/31);  new base C (5/5);  total (42/36)
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.03e-01

GTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAG  >  NC_002947/2095089‑2095381
                                                                                                                                                  |                                                                                                                                                   
gtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGTTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                                                  <  1:1778691/150‑1 (MQ=255)
gtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                                                  <  1:306359/150‑1 (MQ=255)
gtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                                                  <  2:2757710/150‑1 (MQ=255)
 tACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                 >  1:180402/1‑150 (MQ=255)
        gAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          <  2:1766423/150‑1 (MQ=255)
         aGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         <  2:241046/150‑1 (MQ=255)
         aGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                         <  2:1041219/150‑1 (MQ=255)
          gTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATc                                                                                                                                        <  1:499051/150‑1 (MQ=255)
          gTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATc                                                                                                                                        <  1:278184/150‑1 (MQ=255)
                 agagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGccc                                                                                                                                 >  1:1293911/1‑150 (MQ=255)
                 agagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGccc                                                                                                                                 >  1:360012/1‑150 (MQ=255)
                       tGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCccaccac                                                                                                                           <  2:960957/150‑1 (MQ=255)
                               gggCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCcg                                                                                                                   >  1:2496781/1‑150 (MQ=255)
                               gggCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATGCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCGCGGCGCTCGGCGCTCGATCTGCGCCCCACCACCAcccccc                                                                                                                   >  2:1781791/1‑149 (MQ=255)
                                ggCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGcccc                                                                                                                                >  2:347136/1‑136 (MQ=255)
                                ggCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGcccc                                                                                                                                <  1:347136/136‑1 (MQ=255)
                                   tgtgGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTc                                                                                                              <  2:2496781/150‑1 (MQ=255)
                                         aTAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  1:2587300/1‑115 (MQ=255)
                                         aTAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            <  2:2587300/115‑1 (MQ=255)
                                          tAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACctc                                                                                                       >  1:2679937/1‑150 (MQ=255)
                                           agagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                 <  2:1487740/108‑1 (MQ=255)
                                           agagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                 >  1:1487740/1‑108 (MQ=255)
                                             agCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCg                                                                                                    >  2:1487689/1‑150 (MQ=255)
                                                   gCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  2:1696809/1‑105 (MQ=255)
                                                   gCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            <  1:1696809/105‑1 (MQ=255)
                                                      ggATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGGCGAACCTCTCGGTGCCCCCg                                                                                           >  2:23791/1‑150 (MQ=255)
                                                       gATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCCTCGAAACCCCCGTAGCCCCCgg                                                                                          >  2:555752/1‑150 (MQ=255)
                                                       gATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            <  2:925628/101‑1 (MQ=255)
                                                       gATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  1:925628/1‑101 (MQ=255)
                                                        aTGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGgc                                                                                         >  1:973424/1‑150 (MQ=255)
                                                           tAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGcc                                                                                      >  1:2769843/1‑150 (MQ=255)
                                                             aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAAc                                                                                                          <  1:1007154/128‑1 (MQ=255)
                                                             aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAAc                                                                                                          >  2:1007154/1‑128 (MQ=255)
                                                             aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                  >  1:578707/1‑136 (MQ=255)
                                                             aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                  <  2:578707/136‑1 (MQ=255)
                                                             aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCCCGTCGGACCTCCCGTAGCCCCCGGcccccc                                                                                    >  1:1813604/1‑150 (MQ=255)
                                                               aaaCATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGgc                                                                                         <  2:2508780/143‑1 (MQ=255)
                                                               aaaCATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGgc                                                                                         >  1:2508780/1‑143 (MQ=255)
                                                                 aCATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGATCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCACCGTGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTc                                                                                <  1:1781791/150‑1 (MQ=255)
                                                                  cATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         >  2:318194/1‑93 (MQ=255)
                                                                  cATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         <  1:318194/93‑1 (MQ=255)
                                                                   aTTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAAc                                                                                                          <  1:964060/122‑1 (MQ=255)
                                                                   aTTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAAc                                                                                                          >  2:964060/1‑122 (MQ=255)
                                                                       ggCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAAccc                                                                          >  1:1372308/1‑150 (MQ=255)
                                                                        gCTTTGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGGTCTGCGCCCCACCCCCCCCCCCGCCTCCGAACCCTCCTAGCCCCCCCCGCCCCCCTCCAAcccc                                                                         >  2:1538747/1‑149 (MQ=37)
                                                                                gCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCTCCCGTCGAACCCTACTGCGaa                                                                 <  2:360012/150‑1 (MQ=255)
                                                                                      tCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCCCGTCGAACCCCCCGTAGCCCCCCGCGCCCCCCTCCAACCCCACTCCGAATAAccc                                                           >  2:3062856/1‑150 (MQ=38)
                                                                                         aGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTCCTGCGAATAACCCTTc                                                        <  2:1293911/150‑1 (MQ=255)
                                                                                         aGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTc                                                        >  2:2738969/1‑150 (MQ=255)
                                                                                          ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       <  1:896753/150‑1 (MQ=255)
                                                                                           gCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAg                                                      <  2:1679890/150‑1 (MQ=255)
                                                                                           gCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAg                                                      <  2:2679937/150‑1 (MQ=255)
                                                                                              gCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  2:2432212/1‑62 (MQ=21)
                                                                                              gCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            <  1:2432212/62‑1 (MQ=21)
                                                                                               cTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCACCCCCCCGGCCAAACCCTCGTAGACCCCGGCGCCCCCGTCGAACCCTACCGCGGAAAACACTTCAGGgaa                                                  >  1:733521/1‑150 (MQ=37)
                                                                                               cTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGaaa                                                  <  1:1487689/150‑1 (MQ=255)
                                                                                                      tttGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGa                                           >  2:587056/1‑150 (MQ=255)
                                                                                                        tGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATc                                         <  1:1538747/150‑1 (MQ=255)
                                                                                                            gCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCCCCACCACCCCCGCGTCCAACCCCTCGTAGCCCCCGGCGCCCCCCTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCggg                                     >  2:1734993/1‑149 (MQ=255)
                                                                                                            gCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       >  2:717993/1‑132 (MQ=255)
                                                                                                            gCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       <  1:717993/132‑1 (MQ=255)
                                                                                                              aaGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       <  1:2417793/130‑1 (MQ=255)
                                                                                                              aaGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       >  2:2417793/1‑130 (MQ=255)
                                                                                                                gtggtgtggtGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAAt                                 <  1:2888331/150‑1 (MQ=255)
                                                                                                                 tggtgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCACCCCCGCGTCCAACCCCTCCTAACCCCCGGCGCCCCCGCCGAACCCTACTGCGAATAACCCTTTTGGAAATTGGCCATCCGGGCCAtt                                >  1:2629493/1‑150 (MQ=37)
                                                                                                                 tggtgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAAtt                                <  2:973424/150‑1 (MQ=255)
                                                                                                                   gtgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGc                              <  1:3062856/150‑1 (MQ=255)
                                                                                                                        gtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCACCCCCGCGTCCAACCTCTCGTAGCCCCCCGCGCCCCCGGCGAACCCTGCTGCGGAGAACCCCTCAGGAAATTGCCGATGCCGCCAATTGCCACtt                         >  1:931640/1‑150 (MQ=37)
                                                                                                                        gtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                  <  1:321553/77‑1 (MQ=255)
                                                                                                                        gtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                  >  2:321553/1‑77 (MQ=255)
                                                                                                                         tgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGaa                                                                 <  1:164119/109‑1 (MQ=255)
                                                                                                                         tgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGaa                                                                 >  2:164119/1‑109 (MQ=255)
                                                                                                                              gtgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCCGCCCATTGCCCCTTCCAGGt                   >  2:883526/1‑150 (MQ=255)
                                                                                                                                gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCa                 <  2:2524862/150‑1 (MQ=255)
                                                                                                                                gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAAACCTTCAGGAAATTGCCGATTCGGCCAATTGCCGCTTTCAgggta                 >  2:2663785/1‑147 (MQ=255)
                                                                                                                                gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCCCGTCGAACCTCTCGTAGCCCCCCGCGCCCCCCTCGAAACCTACTGCCAAAAACCCCTCAGGAAAATGCCGACCCGGCCAATTTCCCCTTCCCgggcc                 >  2:1995182/1‑147 (MQ=255)
                                                                                                                                    tCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCCAACCCGCCGCGGAAATCTCGCCTCCCCCGCCCCCCCCCCCGAAACCTACTGCGCATAACCCTTCAGGAAATTGCCGATCCCGCCAAATCCCGCCTCCAGGGCATccc             >  2:2816814/1‑149 (MQ=255)
                                                                                                                                     tgAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCAc            >  2:2113121/2‑150 (MQ=255)
                                                                                                                                      gAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       <  1:978931/106‑1 (MQ=255)
                                                                                                                                      gAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                       >  2:978931/1‑106 (MQ=255)
                                                                                                                                        gcgcCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCcg         >  2:1787617/1‑150 (MQ=255)
                                                                                                                                        gcgcCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCcg         >  1:1419930/1‑150 (MQ=255)
                                                                                                                                         cgccgcGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCAGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCcgc        <  2:2769843/150‑1 (MQ=255)
                                                                                                                                              cgcgGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGTaa   <  1:555752/150‑1 (MQ=255)
                                                                                                                                               gcgggcggCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAg  <  2:1813604/150‑1 (MQ=255)
                                                                                                                                               gcgggcggCGCTCGATCTGCGCCCCACCACCACCCCCGCGTAGCACGATTTGCTTCGGCTCGCTGTTGCTTGCAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAg  <  1:2663785/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAG  >  NC_002947/2095089‑2095381

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