Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,351,980 A→C 6.7% intergenic (‑233/+31) PP_4707 ← / ← pnp OsmY‑like protein/polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,351,9800AC 6.7% 172.3 / 5.1 60intergenic (‑233/+31)PP_4707/pnpOsmY‑like protein/polyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base C (2/2);  total (31/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 = 4.11e-01

TTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGG  >  NC_002947/5351838‑5352126
                                                                                                                                              |                                                                                                                                                  
ttCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGAAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGcc                                                                                                                                              >  5:265010/1‑149 (MQ=255)
   tGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCtt                                                                                                                                           >  2:82024/1‑149 (MQ=255)
           ccTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCttttttcat                                                                                                                                    <  3:309118/148‑7 (MQ=255)
           ccTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTAt                                                                                                                                    <  1:183928/148‑1 (MQ=255)
              tttAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCg                                                                                                                                >  7:151354/1‑149 (MQ=255)
              tttAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCg                                                                                                                                >  4:56825/1‑149 (MQ=255)
               ttAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGc                                                                                                                                 <  2:110766/147‑1 (MQ=255)
               ttAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGt                                                                                                                               >  5:317537/1‑149 (MQ=255)
                  aaTCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTGTCTCTATTGCGTGAc                                                                                                                            >  8:172237/1‑149 (MQ=255)
                    tCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGt                                                                                                                          >  4:2640/1‑149 (MQ=255)
                      aaTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAGTCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAGGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTc                                                                                                                        >  5:356984/1‑149 (MQ=255)
                      aaTAAGTTAGCAAATATCTATAACTTCGGTGGTCTTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGt                                                                                                                          >  6:133802/1‑147 (MQ=255)
                       aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCg                                                                                                                       >  4:97390/1‑149 (MQ=255)
                             tcgCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGCTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAg                                                                                                                  <  5:209785/146‑1 (MQ=255)
                                    tatctatcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGa                                                                                                          <  6:137932/149‑1 (MQ=255)
                                            aCTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTg                                                                                                   <  3:56825/148‑1 (MQ=255)
                                             cTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGc                                                                                                 >  3:19663/1‑149 (MQ=255)
                                                cGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGcc                                                                                                <  3:34910/147‑1 (MQ=255)
                                                 ggtggtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCa                                                                                             <  7:67648/149‑1 (MQ=255)
                                                  gtggtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCg                                                                                                                       <  6:367950/122‑1 (MQ=255)
                                                     gtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGGTTAGACGCCCGAGGCCTTGGCCGCAGGGa                                                                                         <  5:133802/149‑1 (MQ=255)
                                                     gtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAggg                                                                                          <  4:19663/148‑1 (MQ=255)
                                                     gtCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAggg                                                                                          <  6:265010/148‑1 (MQ=255)
                                                      tCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGAc                                                                                        >  5:303461/1‑149 (MQ=255)
                                                            aTTTTGCAATCAGCGCGGTAACCCTGCATGCAAGATGCAAGTCCGCAAAAGGTGCCCGTCCCCATGAAAAACGGGCCCGGCGACGGGCCCTTGCTCTATGGCGTGACGTTCGCTTAGCCGCCCGAGGACTTGGCCGCAGGGCCggaggt                                                                                  >  4:271943/1‑144 (MQ=255)
                                                            aTTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCttt                                                                                  <  3:122390/149‑1 (MQ=255)
                                                                 gCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGa                                                                                                                             >  5:56294/1‑101 (MQ=255)
                                                                 gCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGa                                                                                                                             <  6:56294/101‑1 (MQ=255)
                                                                    aTCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGa                                                                            >  3:142734/1‑147 (MQ=255)
                                                                     tCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCGTTGCTCGACAg                                                                         >  8:294172/1‑149 (MQ=255)
                                                                          gcgGTAAACCTGCATGCAAGATGCAAGGCCGCAAAAGGTTCCCGTCCCGATGAAATACGGGCCCTGCGTAGGGCCCTTTCTCTTTTGCGTGACGTTGGGTTAGCCGCGGGAGGCCTTGGCGGCAGGGTCGTCTTTGATCGAGAGCTTgt                                                                    <  4:29791/149‑2 (MQ=255)
                                                                           cgtgaaaCCTGCATGCAAGATGCAAGGCCGCAAAAGGTGCCCGTCCCCATGACAAAAGGGCCCTGCGAAGGGCCCTTGCTCTATTGCGTGACGTTCGCTTAGCCGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGAt                                                                   <  3:271943/145‑1 (MQ=255)
                                                                                     gcatgcaAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGc                                                                                                                                 <  5:164873/77‑1 (MQ=255)
                                                                                     gcatgcaAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTGTCTCTATTGc                                                                                                                                 >  6:164873/1‑77 (MQ=255)
                                                                                          caagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTc                                                    <  7:120867/149‑1 (MQ=255)
                                                                                           aagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACa                                                                          >  4:191690/1‑126 (MQ=255)
                                                                                           aagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACa                                                                          <  3:191690/126‑1 (MQ=255)
                                                                                                    gTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACTTCCAg                                           >  4:329698/1‑148 (MQ=255)
                                                                                                       cGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTtct                                                                                                                                        <  3:138378/52‑1 (MQ=255)
                                                                                                       cGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTtct                                                                                                                                        >  4:138378/1‑52 (MQ=255)
                                                                                                        gCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACaa                                      <  6:317537/149‑1 (MQ=255)
                                                                                                          aaaaGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAg                                     >  5:349502/1‑148 (MQ=255)
                                                                                                                ttCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTc                               >  7:10937/1‑148 (MQ=255)
                                                                                                                     gTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTc                         >  7:145974/1‑149 (MQ=255)
                                                                                                                        cccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTg                      <  4:142734/149‑1 (MQ=255)
                                                                                                                        cccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTg                      <  8:151354/149‑1 (MQ=255)
                                                                                                                         ccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                  >  8:216910/1‑72 (MQ=255)
                                                                                                                         ccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                  <  7:216910/72‑1 (MQ=255)
                                                                                                                           aTGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGAt                                                                               >  4:327719/1‑89 (MQ=255)
                                                                                                                           aTGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGAt                                                                               <  3:327719/89‑1 (MQ=255)
                                                                                                                           aTGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAg                                           <  8:218796/125‑1 (MQ=255)
                                                                                                                           aTGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAg                                           >  7:218796/1‑125 (MQ=255)
                                                                                                                              aaaaaaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGccttc                >  3:37759/1‑149 (MQ=255)
                                                                                                                                aaaaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGccttcct              <  7:172237/149‑1 (MQ=255)
                                                                                                                                  aaGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGccttcctt             >  3:152925/1‑148 (MQ=255)
                                                                                                                                   aGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGa           <  3:97390/149‑1 (MQ=255)
                                                                                                                                     ggCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTa                                         >  3:112571/1‑117 (MQ=255)
                                                                                                                                     ggCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTa                                         <  4:112571/117‑1 (MQ=255)
                                                                                                                                            cGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCgg  >  4:142664/1‑149 (MQ=255)
                                                                                                                                            cGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCgg  >  4:116729/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCGG  >  NC_002947/5351838‑5352126

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