Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,621,920 G→A 7.0% L264L (CTG→CTA) ‡ PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,621,9200GA 7.0% 279.0 / 7.5 95L264L (CTG→CTA) ‡PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (41/47);  new base A (5/2);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAG  >  NC_002947/1621771‑1622061
                                                                                                                                                     |                                                                                                                                             
gacTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTg                                                                                                                                               <  1:1039385/150‑1 (MQ=255)
gacTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTg                                                                                                                                               <  2:1248230/150‑1 (MQ=255)
 acTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTgg                                                                                                                                              <  2:584555/150‑1 (MQ=255)
   ttGCGAGCGTTGCCGGTAATGCAGCAGGAGCGCGAGCTGGGCCCGCTGGTGCGGGCCTGGAACCACTTTACCGAGCGCTTGCGCGGCAAGTTCGCGCGCCAGTCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTggcg                                                                                                                                            >  2:1728205/1‑150 (MQ=255)
    tGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTggcgg                                                                                                                                           <  2:905377/150‑1 (MQ=255)
        aGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCCAGCGGCGtt                                                                                                                                       <  1:631684/150‑1 (MQ=255)
           gTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                    <  1:1840163/150‑1 (MQ=255)
           gTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCCCGACTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                    <  1:1467952/150‑1 (MQ=255)
            ttGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCCAGCGGCGttcgca                                                                                                                                   <  2:2361781/150‑5 (MQ=255)
               ccGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCa                                                                                                                                >  2:1493255/1‑150 (MQ=255)
                cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                               >  2:1154527/1‑150 (MQ=255)
                cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                               >  1:1204038/1‑150 (MQ=255)
                cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                               <  1:1493255/150‑1 (MQ=255)
                cGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                               <  1:776014/150‑1 (MQ=255)
                     gtgcaggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTc                                                                                                                          <  1:972271/150‑1 (MQ=255)
                          ggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGATCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCt                                                                                                                     <  1:1069316/150‑1 (MQ=255)
                          ggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGATGGCACGTGTCGAGCt                                                                                                                     <  1:1174486/150‑1 (MQ=255)
                          ggtgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCAACGCACGCCGCTAGTTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCt                                                                                                                     <  1:1855732/150‑1 (MQ=255)
                            tgcagCGCGAGCTTGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTgg                                                                                                                   <  2:1624905/150‑1 (MQ=255)
                            tgcagCGCGAGCTGGGCCCGCTGGTTCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTgg                                                                                                                   <  2:1340553/150‑1 (MQ=255)
                            tgcagCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTgg                                                                                                                   <  1:1625209/150‑1 (MQ=255)
                            tgcagCGCGAGCTGGGACCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTgg                                                                                                                   <  2:952596/150‑1 (MQ=255)
                                gcgcGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTGTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCtt                                                                                                               >  1:193075/1‑150 (MQ=255)
                                 cgcgAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTg                                                                                                              <  1:79158/150‑1 (MQ=255)
                                 cgcgAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTg                                                                                                              >  2:1516308/1‑150 (MQ=255)
                                     aGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCt                                                                                                          <  2:1286336/150‑1 (MQ=255)
                                      gCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTc                                                                                                         >  1:2272603/1‑150 (MQ=255)
                                        tGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCgg                                                                                                       <  2:1260688/150‑1 (MQ=255)
                                          ggCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGcc                                                                                                     >  2:2157691/1‑150 (MQ=255)
                                            cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                               >  2:1764213/1‑122 (MQ=255)
                                            cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCAc                                                                                                                               <  1:1764213/122‑1 (MQ=255)
                                            cccGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGa                                                                                                   <  2:90314/150‑1 (MQ=255)
                                               gCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGa                                                                                                   <  1:1572940/147‑1 (MQ=255)
                                               gCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGa                                                                                                   >  2:1572940/1‑147 (MQ=255)
                                                   gTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCGCGCAGGGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAg                                                                                            >  1:151571/1‑150 (MQ=255)
                                                   gTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAg                                                                                            <  2:2220296/150‑1 (MQ=255)
                                                   cagcgggcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAg                                                                                            <  1:207800/144‑1 (MQ=255)
                                                     gcgcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGa                                                                                          <  2:852737/150‑1 (MQ=255)
                                                     gcgcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGa                                                                                          <  2:1204038/150‑1 (MQ=255)
                                                      cgcgcCTAGAACCACTCCACCCAGGGCTTGCGCCGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAg                                                                                         <  2:1124528/150‑1 (MQ=255)
                                                       gcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGt                                                                                        >  2:2344172/1‑150 (MQ=255)
                                                       gcgcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGt                                                                                        >  2:2117325/1‑150 (MQ=255)
                                                       gagcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTTCGCTCGGCCGAGCCGCAGGAGt                                                                                        >  1:2020722/3‑150 (MQ=255)
                                                         gcCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTgg                                                                                      >  2:1174259/1‑150 (MQ=255)
                                                          ccTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCtt                                                                                                               >  2:499414/1‑124 (MQ=255)
                                                          ccTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCtt                                                                                                               <  1:499414/124‑1 (MQ=255)
                                                          ccTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTggc                                                                                     >  2:60930/1‑150 (MQ=255)
                                                           cTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCGCGAACTGCGAACGCCGCGGGCGGCGTTGAAGGAAGGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTggcg                                                                                    >  1:748613/1‑150 (MQ=255)
                                                                    cTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCtt                                                                           >  2:2032483/1‑150 (MQ=255)
                                                                     tttACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTg                                                                          <  2:468989/150‑1 (MQ=255)
                                                                      ttACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTgg                                                                         <  2:969743/150‑1 (MQ=255)
                                                                       tACCGAGCGCTTGCGCGGCCAGTTCGCGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGa                                                                        <  2:521906/150‑1 (MQ=255)
                                                                          cGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGCGCAGGTCGTCCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTc                                                                     >  2:63544/1‑150 (MQ=255)
                                                                           gatcgcTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCg                                                                    >  2:2192119/4‑150 (MQ=255)
                                                                           gAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                    <  1:1311404/86‑1 (MQ=255)
                                                                           gAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCg                                                                    <  1:2344172/150‑1 (MQ=255)
                                                                           gAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACCGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                    >  2:1311404/1‑86 (MQ=255)
                                                                            aGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCgg                                                                   <  1:1436869/150‑1 (MQ=255)
                                                                                cTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGc                                                               <  1:212174/150‑1 (MQ=255)
                                                                                 ttGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCa                                                                                 >  2:2255201/1‑131 (MQ=255)
                                                                                 ttGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCa                                                                                 <  1:2255201/131‑1 (MQ=255)
                                                                                 ttGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGcc                                                              >  2:1180865/1‑150 (MQ=255)
                                                                                    cgcgGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAg                                                           <  1:2159963/150‑1 (MQ=255)
                                                                                     gcgGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGGCGGCAGCCCAgg                                                          >  1:1948567/1‑150 (MQ=255)
                                                                                          cAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCAcc                                                     >  1:816640/1‑150 (MQ=255)
                                                                                           aGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCg                                                    >  1:1815559/1‑150 (MQ=255)
                                                                                                gAGCGCCAGGGGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCGGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAgg                                               <  2:496728/150‑1 (MQ=255)
                                                                                                gAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAgg                                               <  2:2210093/150‑1 (MQ=255)
                                                                                                    gcCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGa                                           >  1:321690/1‑150 (MQ=255)
                                                                                                    gcCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGa                                           >  1:321709/1‑150 (MQ=255)
                                                                                                      cAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCa                                                            <  2:162007/131‑1 (MQ=255)
                                                                                                      cAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCa                                                            >  1:162007/1‑131 (MQ=255)
                                                                                                      cAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGAcc                                         >  1:134575/1‑150 (MQ=255)
                                                                                                      cAGGCGCAGTTCATTGCCGACGCGGCCAAAGAACTGCGAACGCCTCTGGCGGCGTTGAAGGCACGTGTCGAGGTTGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGGcc                                         >  2:277940/1‑150 (MQ=255)
                                                                                                       aGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGCGGCGGCAGACCTGGGAGTCTGCAGCCCAGGGCACCGCACGGCTGaaaa                                        >  2:1000508/1‑147 (MQ=255)
                                                                                                        ggCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAg                                                                                                                                   >  2:891990/1‑58 (MQ=255)
                                                                                                        ggCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAg                                                                                                                                   <  1:891990/58‑1 (MQ=255)
                                                                                                             aGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTgg                                  >  1:1355969/1‑150 (MQ=255)
                                                                                                              gTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGc                                 <  1:1516308/150‑1 (MQ=255)
                                                                                                              gTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGc                                 >  1:1580068/1‑150 (MQ=255)
                                                                                                               ttCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGcc                                >  2:409658/1‑150 (MQ=255)
                                                                                                                   ttGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGACAAca                            >  1:1173634/1‑149 (MQ=255)
                                                                                                                      ccGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                    >  2:2288094/1‑43 (MQ=255)
                                                                                                                      ccGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGaa                                                                                                                                    <  1:2288094/43‑1 (MQ=255)
                                                                                                                            cgGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCtgttgt                   <  1:409658/150‑1 (MQ=255)
                                                                                                                               cccACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTccc                <  1:2117325/150‑1 (MQ=255)
                                                                                                                                ccACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGc                                                                  >  2:1123794/1‑99 (MQ=255)
                                                                                                                                ccACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGc                                                                  <  1:1123794/99‑1 (MQ=255)
                                                                                                                                 cACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTg              <  1:1928252/150‑1 (MQ=255)
                                                                                                                                 cACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTg              <  1:2192119/150‑1 (MQ=255)
                                                                                                                                  aCGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCTGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTCTTGTCCCTgg             >  2:1374982/1‑150 (MQ=255)
                                                                                                                                       ccttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGTCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCCGGCCAACCAGCTGTTGTCCCTGGCGCgg        >  2:352/6‑150 (MQ=255)
                                                                                                                                       cTGCGAACGCCGCTGGCGGCTTTGAAGGCACGTGTCGAGCTGGGCTTGCTCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCgg        <  1:1166472/150‑1 (MQ=255)
                                                                                                                                        cttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCggg       >  2:1663175/5‑150 (MQ=255)
                                                                                                                                         ttcaaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCAGCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCgggg      >  2:1618000/4‑149 (MQ=255)
                                                                                                                                         gCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGt      <  2:2287634/150‑1 (MQ=255)
                                                                                                                                            aaCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTc                                                                     >  1:1586677/1‑84 (MQ=255)
                                                                                                                                            aaCGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTc                                                                     <  2:1586677/84‑1 (MQ=255)
                                                                                                                                            aaCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCga   >  2:1092308/1‑150 (MQ=255)
                                                                                                                                             aCGCCGCCAGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCGGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCgag  >  2:1448527/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GACTTGCGAGCGTTGCCGGTGGTGCAGGTGCAGCGCGAGCTGGGCCCGCTGGTGCGCGCCTTGAACCACTTTACCGAGCGCTTGCGCGGCCAGTTCGAGCGCCAGGCGCAGTTCATTGCCGACGCGGCCCACGAACTGCGAACGCCGCTGGCGGCGTTGAAGGCACGTGTCGAGCTGGGCTTGCGCTCGGCCGAGCCGCAGGAGTGGCGGCAGACCTTGGAGTCGGCAGCCCAGGGCACCGACAGGCTGACCCACCTGGCCAACCAGCTGTTGTCCCTGGCGCGGGTCGAG  >  NC_002947/1621771‑1622061

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