Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,118,195 T→G 6.7% intergenic (+27/‑32) PP_4504 → / → PP_4505 hypothetical protein/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,118,1950TG 6.7% 208.4 / 6.6 75intergenic (+27/‑32)PP_4504/PP_4505hypothetical protein/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base G (3/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

GTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGC  >  NC_002947/5118048‑5118333
                                                                                                                                                   |                                                                                                                                          
gTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAt                                                                                                                                            >  5:182716/1‑148 (MQ=255)
      cACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcg                                                                                                                                     >  1:113927/1‑149 (MQ=255)
       aCGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgg                                                                                                                                    >  5:198578/1‑149 (MQ=255)
         ggCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgggc                                                                                                                                  <  1:289357/149‑1 (MQ=255)
         ggCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCACATGCCCGGCCCGCATCACTGcgggc                                                                                                                                  >  1:65588/1‑149 (MQ=255)
            tGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCAAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgggccgg                                                                                                                               >  6:183403/1‑149 (MQ=255)
              aCCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgggccgggc                                                                                                                             >  8:138146/1‑149 (MQ=255)
              aCCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgggccgg                                                                                                                               <  6:152744/147‑1 (MQ=255)
              aCCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTCaggtccgggc                                                                                                                             >  8:133210/1‑149 (MQ=255)
                cagcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCtt                                                                                                                           >  8:305750/1‑149 (MQ=255)
                cagcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgggcat                                                                                                                           >  1:224946/1‑147 (MQ=255)
                cagcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgggcat                                                                                                                           <  7:209292/149‑3 (MQ=255)
                cagcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGAAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCtt                                                                                                                           <  5:103994/149‑1 (MQ=255)
                 agcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCttt                                                                                                                          >  7:49530/1‑149 (MQ=255)
                 agcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgggcatg                                                                                                                          <  2:224946/149‑4 (MQ=255)
                 agcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCCGCCCGCCTCACTGcgggccggcgttt                                                                                                                          >  8:330744/1‑144 (MQ=255)
                  gcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCtttt                                                                                                                         >  8:168462/1‑149 (MQ=255)
                   cagAACAAGTTTGACCAACTCCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGCCCTGGAAAGCTGAAAAAGGCGCCAAGCCCCAGTAAGCTGTCCGCCCCTTGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTc                                                                                                                        <  6:43039/149‑1 (MQ=255)
                      aagaagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCg                                                                                                                     >  4:305154/1‑149 (MQ=255)
                        gaagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGtt                                                                                                                   >  2:196619/1‑149 (MQ=255)
                         aagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGttt                                                                                                                  <  5:184364/149‑1 (MQ=255)
                           gTTTGACGAACTGCAAACAGAACGCGCCGCCAAGAACCCCGAGGGGCAGGAGTGCCAGGCCGGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTag                                                                                                                <  5:311205/149‑1 (MQ=255)
                           gTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGc                                                                                                                                      >  3:124793/1‑127 (MQ=255)
                           gTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGc                                                                                                                                      <  4:124793/127‑1 (MQ=255)
                           gTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTGCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATCCCCGCCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTag                                                                                                                <  5:240147/149‑1 (MQ=255)
                            tttGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTagg                                                                                                               <  3:362159/149‑1 (MQ=255)
                               gACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTaggag                                                                                                             <  3:305154/148‑1 (MQ=255)
                               gACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCACGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTaggag                                                                                                             <  8:17939/148‑1 (MQ=255)
                                    aCTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCAc                                                                                                       <  2:91987/149‑1 (MQ=255)
                                       gCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCAAGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTg                                                                                                    >  3:113482/1‑149 (MQ=255)
                                       gCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCAAGCCCCGCCCGCATCACTGCGGCCCGCGCGTATCGGGTGTAGGCGATGCCCTTg                                                                                                    >  3:213979/1‑149 (MQ=255)
                                        cAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACtt                                                                                                     <  6:182716/147‑1 (MQ=255)
                                        cAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGc                                                                                                   >  6:189848/1‑149 (MQ=255)
                                        cAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGc                                                                                                   <  8:182035/149‑1 (MQ=255)
                                           aaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt                                                                                                <  7:307236/149‑1 (MQ=255)
                                           aaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt                                                                                                >  6:48204/1‑149 (MQ=255)
                                           aaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt                                                                                                <  7:168462/149‑1 (MQ=255)
                                                    gccgccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAACGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCt                                                                                       >  3:226722/1‑149 (MQ=255)
                                                     ccgccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTg                                                                                      <  8:2855/149‑1 (MQ=255)
                                                        ccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGgcg                                                                                   >  6:137698/1‑149 (MQ=255)
                                                               aGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTg                                                                            <  7:305750/149‑1 (MQ=255)
                                                               aGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTg                                                                            <  4:5116/149‑1 (MQ=255)
                                                                 ccGAGTGGCAGGAGTTCCAGGCCTCGAAAGCTGAAAAACCCGCCAAGGCCCAGTAACCAGTCCGCCCCATGCCCGGCCCGCATCACTGCGGTCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGc                                                                          <  4:300947/149‑1 (MQ=255)
                                                                   gAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCCAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACt                                                                                                      >  5:256387/1‑119 (MQ=255)
                                                                   gAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACt                                                                                                      <  6:256387/119‑1 (MQ=255)
                                                                        gCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTc                                                                   >  1:272780/1‑149 (MQ=255)
                                                                        gCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATAACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTc                                                                   >  1:343677/1‑149 (MQ=255)
                                                                         cAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTct                                                                  <  5:189848/149‑1 (MQ=255)
                                                                         cAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTct                                                                  <  5:48204/149‑1 (MQ=255)
                                                                          aGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctg                                                                 >  2:167756/1‑149 (MQ=255)
                                                                           ggAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctgg                                                                <  3:110967/149‑1 (MQ=255)
                                                                               ttCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctggctgg                                                            >  8:98346/1‑149 (MQ=255)
                                                                                tCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCtggctggc                                                           >  5:366838/1‑149 (MQ=255)
                                                                                        tGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCAccc                                                   >  7:62977/1‑149 (MQ=255)
                                                                                          gAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCt                                                  <  8:49530/148‑1 (MQ=255)
                                                                                          gAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCAAGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTg                                                 <  4:113482/149‑1 (MQ=255)
                                                                                           aaaGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCa                                                                                 >  7:292926/1‑116 (MQ=255)
                                                                                           aaaGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCa                                                                                 <  8:292926/116‑1 (MQ=255)
                                                                                           aaaGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTgg                                                <  1:167756/149‑1 (MQ=255)
                                                                                                tGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCg                                             >  3:262655/1‑147 (MQ=255)
                                                                                                  aaaaaGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAg                                          <  6:198578/148‑1 (MQ=255)
                                                                                                   aaaaGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGc                                        <  5:183403/149‑1 (MQ=255)
                                                                                                         cgcCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCg                                  <  2:343677/149‑1 (MQ=255)
                                                                                                                gCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCctgc                           <  7:133210/149‑1 (MQ=255)
                                                                                                                gCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCctgc                           <  7:138146/149‑1 (MQ=255)
                                                                                                                  ccAGTAAGCTGTCCGCCCGATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCctgctg                         <  7:330744/149‑1 (MQ=255)
                                                                                                                  ccAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCctgctg                         >  4:108135/1‑149 (MQ=255)
                                                                                                                   cAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTgg                        <  7:98346/149‑1 (MQ=255)
                                                                                                                     gTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTgggg                      <  7:213987/149‑1 (MQ=255)
                                                                                                                          cTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACAt                 <  4:226722/149‑1 (MQ=255)
                                                                                                                             tCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATg                <  2:113927/147‑1 (MQ=255)
                                                                                                                                  cccATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAAc          >  2:158030/1‑148 (MQ=255)
                                                                                                                                    cATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAg       <  1:158030/149‑1 (MQ=255)
                                                                                                                                        cccggcccgCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGAc    <  1:196619/148‑1 (MQ=255)
                                                                                                                                         ccggcccgCAGTGATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGc  >  2:200288/1‑149 (MQ=255)
                                                                                                                                         ccggcccgCAGTGATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGc  >  6:14067/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GTTGCTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGC  >  NC_002947/5118048‑5118333

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