Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,846 T→G 12.8% R112R (CGT→CGG PP_4562 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,8460TG12.8% 135.3 / 13.8 55R112R (CGT→CGGPP_4562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base G (3/4);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

AGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTG  >  NC_002947/5182700‑5182990
                                                                                                                                                  |                                                                                                                                                   
aGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTcc                                                                                                                                                   >  8:163742/1‑149 (MQ=255)
  cGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                                                                                                                                                 <  7:225899/149‑1 (MQ=255)
  cGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                                                                                                                                                 <  7:220162/149‑1 (MQ=255)
   gATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCggg                                                                                                                                                >  4:346727/1‑149 (MQ=255)
     tGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGTCAAGCTCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGCCACCCGGACGGGTg                                                                                                                                              <  4:373711/149‑1 (MQ=255)
     tGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTg                                                                                                                                              <  3:346727/149‑1 (MQ=255)
       aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                          <  5:9065/149‑1 (MQ=255)
       aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                          <  7:163742/149‑1 (MQ=255)
        gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                          >  4:362869/1‑148 (MQ=255)
        gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                         <  5:112467/149‑1 (MQ=255)
        gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                         <  7:52364/149‑1 (MQ=255)
        gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                         <  5:107108/149‑1 (MQ=255)
          tGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTCGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGag                                                                                                                                      >  5:121753/1‑149 (MQ=255)
           gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCGCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCCGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCgcgc                                                                                                                                     <  6:268373/149‑1 (MQ=255)
            cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCGcgca                                                                                                                                    <  6:38399/149‑1 (MQ=255)
                       ggTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCTAGCAGCTGGACTTcc                                                                                                                         <  4:326024/149‑1 (MQ=255)
                            ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCgcgc                                                                                                                                     >  3:301503/1‑132 (MQ=255)
                            ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCgcgc                                                                                                                                     <  4:301503/132‑1 (MQ=255)
                             caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACt                                                                                                                   >  4:2269/1‑149 (MQ=255)
                               ccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGAGCCGACCGGGACTCCCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGCCTTCCGCAACTTc                                                                                                                 >  3:324566/1‑149 (MQ=255)
                                caATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCg                                                                                                                >  7:287107/1‑149 (MQ=255)
                                      ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAgg                                                                                                          >  7:97911/1‑149 (MQ=255)
                                      ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAg                                                                                                           <  6:217587/148‑1 (MQ=255)
                                      ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAg                                                                                                           <  3:50989/148‑1 (MQ=255)
                                       cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGc                                                                                                         <  3:362869/149‑1 (MQ=255)
                                             gCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAgg                                                                                                   <  2:125814/149‑1 (MQ=255)
                                               tCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGcc                                                                                                 >  4:6588/1‑149 (MQ=255)
                                                      gAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                    <  2:221934/139‑1 (MQ=255)
                                                      gAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                    >  1:221934/1‑139 (MQ=255)
                                                           aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtct                                                                                     <  3:2269/149‑1 (MQ=255)
                                                           aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtct                                                                                     >  2:87915/1‑149 (MQ=255)
                                                              ggAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGa                                                                                  >  5:360979/1‑149 (MQ=255)
                                                               gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                 <  6:360979/149‑1 (MQ=255)
                                                               gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                 <  2:47904/149‑1 (MQ=255)
                                                               gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                 <  6:121753/149‑1 (MQ=255)
                                                               gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCCGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                 <  1:65188/149‑1 (MQ=255)
                                                                   tGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGa                                                                             <  8:209790/149‑1 (MQ=255)
                                                                             tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGa                                                                   >  7:31779/1‑149 (MQ=255)
                                                                               gCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGAt                                                                  <  1:87915/148‑1 (MQ=255)
                                                                                cAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGa                                                                <  8:287107/149‑1 (MQ=255)
                                                                                     gccgcTCCACCATCAGACTGGCCCTGATCGCCGACCGTGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATc                                                           <  4:324566/149‑1 (MQ=255)
                                                                                              ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                                  <  6:211966/149‑1 (MQ=255)
                                                                                                             tCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGcc                                    >  5:113925/1‑148 (MQ=255)
                                                                                                                     gACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCctt                           >  3:100654/1‑149 (MQ=255)
                                                                                                                              gtccagctgctCGCCACCCGGACGGGTG‑‑TCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt                  >  3:192586/6‑149 (MQ=38)
                                                                                                                                gctgctGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                >  8:62475/1‑149 (MQ=255)
                                                                                                                                gctgctGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                >  6:381283/1‑149 (MQ=255)
                                                                                                                                ccagctgctCGCCACCCGGACGGGTGTCGCGCAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGt                   >  3:242308/4‑149 (MQ=38)
                                                                                                                                   gctgcGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCt             >  3:137367/1‑149 (MQ=255)
                                                                                                                                     tgcGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg           >  3:112162/1‑149 (MQ=255)
                                                                                                                                     tgcGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTAGAACAACCCGCCAAGGGCCTTCTTGAAGTTGACCtt            >  7:126322/1‑148 (MQ=255)
                                                                                                                                        gcgACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTg        <  4:137367/149‑1 (MQ=255)
                                                                                                                                               ccGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                >  7:265293/1‑134 (MQ=255)
                                                                                                                                               ccGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                <  8:265293/134‑1 (MQ=255)
                                                                                                                                               ccGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGCCTATCTCGATGCGATCGTCATGGATGACGATCAGCCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTg  >  2:145153/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTG  >  NC_002947/5182700‑5182990

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: