Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,839 A→C 13.2% D110A (GAC→GCC)  PP_4562 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,8390AC13.2% 137.8 / 13.2 53D110A (GAC→GCC) PP_4562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base C (3/4);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

CGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTG  >  NC_002947/5182695‑5182984
                                                                                                                                                |                                                                                                                                                    
cGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGacac                                                                                                                                                   <  8:22702/148‑1 (MQ=255)
cGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGacac                                                                                                                                                   >  7:22702/1‑148 (MQ=255)
 ggAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACAccc                                                                                                                                                 >  5:211966/1‑149 (MQ=255)
     aGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTcc                                                                                                                                             >  8:163742/1‑149 (MQ=255)
       cGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                                                                                                                                           <  7:220162/149‑1 (MQ=255)
       cGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                                                                                                                                           <  7:225899/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                                                                                                                                    <  7:163742/149‑1 (MQ=255)
            aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                    <  5:9065/149‑1 (MQ=255)
             gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                    >  4:362869/1‑148 (MQ=255)
             gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                   <  5:112467/149‑1 (MQ=255)
             gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                   <  5:107108/149‑1 (MQ=255)
             gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCg                                                                                                                                   <  7:52364/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                                                                                                                               <  4:301503/132‑1 (MQ=255)
                                 ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCTCGCCACCCGGACGGGTG‑‑TCgcgc                                                                                                                               >  3:301503/1‑132 (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                                                                                                     <  3:50989/148‑1 (MQ=255)
                                           ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAg                                                                                                     <  6:217587/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                                                                                              >  1:221934/1‑139 (MQ=255)
                                                           gAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                              <  2:221934/139‑1 (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:121753/149‑1 (MQ=255)
                                                                    gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                           <  2:47904/149‑1 (MQ=255)
                                                                    gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                           <  6:360979/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          >  6:381283/1‑149 (MQ=255)
                                                                                                                                     gctgctGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc          >  8:62475/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)
                                                                                                                                                |                                                                                                                                                    
CGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTG  >  NC_002947/5182695‑5182984

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: