Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,511,675 G→A 12.3% T19M (ACG→ATG) ‡ copA‑I ← copper resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,511,6750GA12.3% 146.8 / 9.3 57T19M (ACG→ATG) ‡copA‑Icopper resistance protein A
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base A (3/4);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTTGACGGTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGTGTTGTCAGCTTGCTG  >  NC_002947/2511536‑2511819
                                                                                                                                            |                                                                                                                                                
gttgACGGTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCGCCGCCAGCACGCTggcg                                                                                                                                          >  8:106492/1‑149 (MQ=255)
 ttgACGGTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCACGCTggcgg                                                                                                                                         <  6:354138/148‑1 (MQ=255)
 ttgACGGTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCACGCTggcgg                                                                                                                                         <  8:80359/148‑1 (MQ=255)
   gACGGTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCGCCGCCAGCACGCTggcggg                                                                                                                                        <  7:106492/148‑2 (MQ=255)
      ggTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGcc                                                                                                                                     >  8:249431/1‑147 (MQ=255)
       gTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGACGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAg                                                                                                                                  >  6:145880/1‑149 (MQ=255)
                gTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTggcggc                                                                                                                                        <  1:25595/134‑1 (MQ=255)
                gTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTggcggc                                                                                                                                        >  2:25595/1‑134 (MQ=255)
                        tCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTAc                                                                                                                 <  1:318781/149‑1 (MQ=255)
                          ggCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTAcgg                                                                                                               >  8:150933/1‑149 (MQ=255)
                              tgtgATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTAAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTAcggcgg                                                                                                            >  5:324901/1‑148 (MQ=255)
                                tgATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAg                                                                                                         <  5:231136/149‑1 (MQ=255)
                                       gACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTgc                                                                                                  >  1:83685/1‑149 (MQ=255)
                                        aCCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTgc                                                                                                  <  6:324901/148‑1 (MQ=255)
                                        aCCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGACGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTgc                                                                                                  <  5:145880/148‑1 (MQ=255)
                                           ggAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCaa                                                                                              <  8:306708/149‑1 (MQ=255)
                                               tCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCa                                                                                                                                   <  8:288373/108‑1 (MQ=255)
                                               tCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCa                                                                                                                                   >  7:288373/1‑108 (MQ=255)
                                               tCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATg                                                                                          >  4:208326/1‑149 (MQ=255)
                                                      aaTGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAg                                                                                     <  3:161136/147‑1 (MQ=255)
                                                      aaTGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAAAGTACGGCGGGAGGGGTTGAGCAACATGCCGAg                                                                                     >  4:161136/1‑147 (MQ=255)
                                                         gAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGcat                                                                                 >  3:61015/1‑146 (MQ=255)
                                                           aCAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCa                                                                               <  3:208326/148‑1 (MQ=255)
                                                                  gAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAAc                                                                                             >  6:299428/1‑127 (MQ=255)
                                                                  gAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAAc                                                                                             <  5:299428/127‑1 (MQ=255)
                                                                    aaTGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAg                                                                     >  7:100332/1‑149 (MQ=255)
                                                                         tGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCGGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCAt                                                                              >  7:89660/1‑135 (MQ=255)
                                                                         tGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCGGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCAt                                                                              <  8:89660/135‑1 (MQ=255)
                                                                           gccggccAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGAt                                                              >  8:309730/1‑149 (MQ=255)
                                                                             cggccAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCCGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGt                                                            >  1:309287/1‑149 (MQ=255)
                                                                             cggccAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGt                                                            <  7:150933/149‑1 (MQ=255)
                                                                              cgccAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGTCGCCACAGACCCAGGCCCGCCAGCGTGCGGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGtt                                                           <  8:176239/148‑1 (MQ=255)
                                                                               gccAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTc                                                          >  1:163633/1‑149 (MQ=255)
                                                                                     tctcGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTgc                                                     <  8:292868/148‑1 (MQ=255)
                                                                                     tctcGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTgc                                                     >  7:292868/1‑148 (MQ=255)
                                                                                         ggCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCACGCCCGCCAGCGTGCTGGCGGCGACCACGCCTTTGACAAACGAACCGCGGGACGGGGTGCGCAACATGACGAGGCCTCATGCGGAATTGCCACGATGTTCGCTGCCGCct                                                >  4:8831/1‑149 (MQ=255)
                                                                                             ttCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCctgctg                                            >  3:215943/1‑149 (MQ=255)
                                                                                              tCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGt                                           >  3:202647/1‑149 (MQ=255)
                                                                                                 gCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAg                                        <  2:163633/149‑1 (MQ=255)
                                                                                                    gCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAgg                                       >  5:237642/1‑147 (MQ=255)
                                                                                                     ccTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAg                                    >  8:123369/1‑149 (MQ=255)
                                                                                                         gACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCt                                  >  4:53056/1‑147 (MQ=255)
                                                                                                         gACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGa                                <  8:100332/149‑1 (MQ=255)
                                                                                                           ccAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCTTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAgg                              <  4:61015/149‑1 (MQ=255)
                                                                                                                 ggCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGcaca                        <  7:123369/149‑1 (MQ=255)
                                                                                                                      cacaGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACa                   <  1:106681/149‑1 (MQ=255)
                                                                                                                        caGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAgt                 >  8:145572/1‑149 (MQ=255)
                                                                                                                         aGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAgtg                <  7:309730/149‑1 (MQ=255)
                                                                                                                          gACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAgtgt               <  7:249431/149‑1 (MQ=255)
                                                                                                                           aCCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGtgtt              <  6:237642/149‑1 (MQ=255)
                                                                                                                             ctgggcgccGCCAGCACGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGtgttg             >  7:30916/8‑149 (MQ=255)
                                                                                                                             ccAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGtgttgt            >  2:66144/1‑149 (MQ=255)
                                                                                                                              cAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGtgttgt            >  4:336921/1‑148 (MQ=255)
                                                                                                                                   gccGCCAGCACGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGTGTTGTCAGCt       >  1:279232/2‑149 (MQ=255)
                                                                                                                                    ccGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGTGTTGTCAGCtt      >  2:299569/1‑149 (MQ=255)
                                                                                                                                     cGCCAGCACGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGTGTTGTCAGCTtg     <  8:30916/149‑1 (MQ=255)
                                                                                                                                        cAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGTGTTGTCAGCTtgctg  <  7:145572/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GTTGACGGTCAGTGCGGTACGCGGTCGGCCTGTGATGTTGACCGGAGTCTGGCCAATGAACAGCTCGAAATGCTGGCCGGCCAGGTCTCGGCCTTCAGCGGCCTGGACCAGCGGGCGCCACAGACCCAGGCCCGCCAGCGTGCTGGCGGCGCCCAGGCCTTTGACAAACGTACGGCGGGAGGGGTTGCGCAACATGGCGAGGCCTCATGCGGAATAGCCACGATGTTCGCTGCCGCCTGCTGTCAGGTAGCTGAGGTGCACATTACAGTGTTGTCAGCTTGCTG  >  NC_002947/2511536‑2511819

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: