Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,118,196 C→T 6.7% intergenic (+28/‑31) 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,1960CT 6.7% 218.2 / 3.1 75intergenic (+28/‑31)PP_4504/PP_4505hypothetical protein/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (29/40);  new base T (3/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGC  >  NC_002947/5118054‑5118333
                                                                                                                                              |                                                                                                                                         
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                                                                                                <  7:168462/149‑1 (MQ=255)
                                     aaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt                                                                                                >  6:48204/1‑149 (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)
                                                                                                                                              |                                                                                                                                         
CACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGC  >  NC_002947/5118054‑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: ‑