Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,824 A→G 5.3% V690A (GTT→GCT)  PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,684,8240AG 5.3% 276.8 / 4.2 95V690A (GTT→GCT) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base A (48/41);  new base G (3/2);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

GCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCCCTTGTACTCAC  >  NC_002947/5684679‑5684971
                                                                                                                                                 |                                                                                                                                                       
gcgcGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCt                                                                                                                                                      <  2:84340/149‑1 (MQ=255)
gcgcGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCAGGCCAACCt                                                                                                                                                      >  3:33972/1‑149 (MQ=255)
     acaacaACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCAAGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGCCAGCCTGcat                                                                                                                                                  >  6:1852/1‑146 (MQ=37)
      caacaaCGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGa                                                                                                                                                >  2:2687/1‑149 (MQ=255)
       aacaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                                               >  7:324263/1‑149 (MQ=255)
        acaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCGACCTGCTCGATc                                                                                                                                              <  8:186315/149‑1 (MQ=255)
        acaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATc                                                                                                                                              <  7:13748/149‑1 (MQ=255)
        acaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATc                                                                                                                                              <  2:11472/149‑1 (MQ=255)
         caacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCACCCTGCTCGATCa                                                                                                                                             >  8:241238/1‑149 (MQ=255)
         caacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAcc                                                                                                                                                       >  1:217813/1‑139 (MQ=255)
         caacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAcc                                                                                                                                                       <  2:217813/139‑1 (MQ=255)
         caacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCAAGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGCCAGCCTGcatgaaga                                                                                                                                             <  5:1852/149‑8 (MQ=37)
              tctcTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTcc                                                                                                                                        >  8:102256/1‑149 (MQ=255)
                     aTCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCATGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCAc                                                                                                                                 >  3:254805/1‑149 (MQ=255)
                      tCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCAcc                                                                                                                                >  5:45169/1‑149 (MQ=255)
                        gAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCgg                                                                                                                              <  4:178400/149‑1 (MQ=255)
                            aGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGAt                                                                                                                          <  3:137931/149‑1 (MQ=255)
                              gATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATgg                                                                                                                        >  4:275873/1‑149 (MQ=255)
                               aTGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggc                                                                                                                       >  2:355336/1‑149 (MQ=255)
                                tGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcg                                                                                                                      >  3:167013/1‑149 (MQ=255)
                                tGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcg                                                                                                                      >  3:103451/1‑149 (MQ=255)
                                 gTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcgg                                                                                                                     >  2:185913/1‑149 (MQ=255)
                                 gTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcgg                                                                                                                     >  4:114738/1‑149 (MQ=255)
                                  tCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGt                                                                                                                                           >  5:26409/1‑126 (MQ=255)
                                  tCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGt                                                                                                                                           <  6:26409/126‑1 (MQ=255)
                                   ccATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGa                                                                                                                           <  3:254802/141‑1 (MQ=255)
                                   ccATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCATGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCAACCGCCTGGCCAACCTGCTCGATCAGTACCGCGGCACCGGGGa                                                                                                                           >  4:254802/1‑141 (MQ=255)
                                     atggcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAg                                                                                                                 <  4:115159/149‑1 (MQ=255)
                                     atggcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAg                                                                                                                 >  1:137209/1‑149 (MQ=255)
                                         catggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGc                                                                                                             <  7:151114/149‑1 (MQ=255)
                                           tggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGt                                                                                                          <  1:66336/149‑1 (MQ=255)
                                             ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCATGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGc                                                                                                        <  4:254805/149‑1 (MQ=255)
                                             ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTg                                                                                                         <  5:383057/148‑1 (MQ=255)
                                              gCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCt                                                                                                       <  2:137209/149‑1 (MQ=255)
                                                cAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTAt                                                                                                     >  8:94025/1‑149 (MQ=255)
                                                cAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTAt                                                                                                     >  1:53269/1‑149 (MQ=255)
                                                   gAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                                  <  7:102256/149‑1 (MQ=255)
                                                   gAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                                  <  7:94025/149‑1 (MQ=255)
                                                   gAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                                                  <  4:293684/149‑1 (MQ=255)
                                                      gCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGtt                                                                                               <  7:241238/149‑1 (MQ=255)
                                                       cTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTg                                                                                             >  5:236526/1‑149 (MQ=255)
                                                       cTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTg                                                                                             >  1:186732/1‑149 (MQ=255)
                                                          gCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTgc‑‑gc                                                                                        <  4:123809/149‑1 (MQ=255)
                                                          gCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGCTGCGGGGCTTCCTGCAGGCTGGACAGGAGGAGGTCCACGGCCTGTCGAACCGGCTCGATCAGTACCTGGCCTCCGGGGATGGCGGCCAGGGGCGTGCTGTCGAGTTgc‑‑gc                                                                                        >  8:174984/1‑149 (MQ=255)
                                                               cAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGc                                                                                   >  3:86003/1‑149 (MQ=255)
                                                               cAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGc                                                                                   <  7:28933/149‑1 (MQ=255)
                                                                aaCTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGcc                                                                                  >  5:208407/1‑149 (MQ=255)
                                                                  cTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCa                                                                                <  3:275873/149‑1 (MQ=255)
                                                                   tGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGcc                                                                                  >  2:7122/1‑146 (MQ=255)
                                                                     cGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGAt                                                                             <  8:324263/149‑1 (MQ=255)
                                                                       gATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATg                                                                            >  1:101015/1‑148 (MQ=255)
                                                                         tCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATg                                                                                                                         >  4:279663/1‑105 (MQ=255)
                                                                         tCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATg                                                                                                                         <  3:279663/105‑1 (MQ=255)
                                                                             gCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGaacaac                                                                     >  6:363397/1‑149 (MQ=255)
                                                                             gCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGaacaac                                                                     <  5:363397/149‑1 (MQ=255)
                                                                             gCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGaacaac                                                                     >  1:179442/1‑149 (MQ=255)
                                                                             gCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTAGCGAGTTGC‑‑GCCAAGCCCAGATGGaaccac                                                                     >  8:166189/1‑149 (MQ=255)
                                                                                   tgcGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCAAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTAACGCGGCACCGGGGAAGGCGGCCAGGGGCGTGCAATCGAGTTGC‑‑CCGAAGCCAAGATGGAACAACGACcgc                                                               >  6:380734/1‑149 (MQ=255)
                                                                                      ggccggcATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCa                                                            <  4:167013/149‑1 (MQ=255)
                                                                                          ggcATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCt                                                        >  5:358231/1‑149 (MQ=255)
                                                                                            cATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTc                                                       >  7:208037/1‑148 (MQ=255)
                                                                                            aaTGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGACTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTc                                                       >  1:23590/2‑148 (MQ=255)
                                                                                                 tGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGcca                                                 <  8:159141/149‑1 (MQ=255)
                                                                                                  gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccat                                                >  2:350261/1‑149 (MQ=255)
                                                                                                    gggTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccattc                                              <  6:208407/149‑1 (MQ=255)
                                                                                                       tgatcgagCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccattccat                                           >  4:68267/8‑149 (MQ=37)
                                                                                                            aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCgg                                      >  3:253933/1‑149 (MQ=255)
                                                                                                            aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCgg                                      >  3:278551/1‑149 (MQ=255)
                                                                                                            aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCAGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACCACGACCCCGCAGCCTCCAGCCATTCCATTTCgg                                      >  3:280524/1‑149 (MQ=255)
                                                                                                             agCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGc                                     >  8:34884/2‑149 (MQ=38)
                                                                                                              gCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGaa                                                                         >  4:192620/1‑112 (MQ=255)
                                                                                                              gCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGaa                                                                         <  3:192620/112‑1 (MQ=255)
                                                                                                               ccggtgggCCCGGCGGTGGCCCGCCGCCGGGCCAGCCGGCGCGTTCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCa                                   <  6:379617/143‑1 (MQ=255)
                                                                                                               cAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCa                                   >  5:245309/1‑149 (MQ=255)
                                                                                                                aGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAt                                  >  5:166201/1‑149 (MQ=255)
                                                                                                                 ggCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAtt                                 <  1:7122/149‑1 (MQ=255)
                                                                                                                  gCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAttt                                >  5:322548/1‑149 (MQ=255)
                                                                                                                   cTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg                               <  4:98125/149‑1 (MQ=255)
                                                                                                                     ggCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg                               >  1:119065/1‑147 (MQ=255)
                                                                                                                      gCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGc                              <  6:322548/147‑1 (MQ=255)
                                                                                                                      gCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAg                            >  3:158513/1‑149 (MQ=255)
                                                                                                                       ccAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGt                           <  1:2687/149‑1 (MQ=255)
                                                                                                                       ccAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGt                           <  2:119065/149‑1 (MQ=255)
                                                                                                                        cAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAg                                                    <  8:294383/123‑1 (MQ=255)
                                                                                                                        cAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAg                                                    >  7:294383/1‑123 (MQ=255)
                                                                                                                         aGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGTCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTgg                         >  3:218774/1‑149 (MQ=255)
                                                                                                                         aGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTgg                         >  2:316102/1‑149 (MQ=255)
                                                                                                                           gcggcgGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCa                       <  4:33972/149‑1 (MQ=255)
                                                                                                                                ggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGa                  <  1:355336/149‑1 (MQ=255)
                                                                                                                                    cACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATgg               <  1:185913/148‑1 (MQ=255)
                                                                                                                                     aCGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGcc             >  2:354499/1‑149 (MQ=255)
                                                                                                                                     aCGGCCGGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTTGCTATGTAGTTCGCTTCCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGAt                 <  1:303153/149‑1 (MQ=255)
                                                                                                                                        gCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCCCtt          <  2:179442/149‑1 (MQ=255)
                                                                                                                                        gCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCCCtt          <  4:103451/149‑1 (MQ=255)
                                                                                                                                             gCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCCCTTGTACt     >  3:129424/1‑149 (MQ=255)
                                                                                                                                                aaCCGGCTCGCTCAGTTCCGCGGCCCCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GGCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCCGCCATTTGCAGTGGCAGTTGATGGCCCTTGTACTCAc  <  4:280524/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
GCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGC‑‑GCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCCCTTGTACTCAC  >  NC_002947/5684679‑5684971

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