Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,887,061 A→G 10.3% intergenic (‑46/‑84) PP_5161 ← / → PP_5162 lysophospholipase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,887,0610AG10.3% 154.0 / 7.3 58intergenic (‑46/‑84)PP_5161/PP_5162lysophospholipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/4);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

CGGGGTGAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCC  >  NC_002947/5886919‑5887208
                                                                                                                                              |                                                                                                                                                   
cGGGGTGAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAg                                                                                                                                               >  7:66388/1‑149 (MQ=255)
     tGAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGc                                                                                                                                          >  2:105824/1‑149 (MQ=255)
       aaaGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTc                                                                                                                                        >  7:56879/1‑149 (MQ=255)
         aGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTcc                                                                                                                                       >  6:24359/1‑148 (MQ=255)
          ggCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCca                                                                                                                                     <  2:253303/149‑1 (MQ=255)
           gCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCcac                                                                                                                                    <  5:116663/149‑1 (MQ=255)
             aGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGTGGGAGCGGCCTTGCGTCGCGAAAGGCCGCTCCcgcaa                                                                                                                                  <  8:203202/149‑2 (MQ=38)
             aGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTGCGTCGCGAAAGGCCGCTCCcgcaa                                                                                                                                  <  2:118272/149‑2 (MQ=39)
             aGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCGCTTGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAg                                                                                                                                  <  8:146885/149‑1 (MQ=255)
                    aTCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTGCGTCGCGAAAGGcc                                                                                                                                            >  1:65642/1‑132 (MQ=255)
                    aTCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTGCGTCGCGAAAGGcc                                                                                                                                            <  2:65642/132‑1 (MQ=255)
                             aTCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCtt                                                                                                                  <  3:141048/149‑1 (MQ=255)
                             aTCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCtt                                                                                                                  <  2:41983/149‑1 (MQ=255)
                              tCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCttt                                                                                                                 <  8:66388/149‑1 (MQ=255)
                               ccACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCttt                                                                                                                 >  1:152393/1‑148 (MQ=255)
                                cacaTTGATGCAAGTATTGATCTGGTATCTTCAGCTCTGGCCGGCGCCATGTCAAGCTTGCTCACCTCGAACAACCAGCATAGAACTAATTGCGGGAGCGGCCTTTCGCGACGGAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTa                                                                                                               <  6:37708/149‑1 (MQ=255)
                                cacaTTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCttt                                                                                                                 >  6:227896/1‑147 (MQ=255)
                                   aTTGATGCTAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACAAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGCCGCCAGCGTTTAGa                                                                                                             >  2:320069/1‑148 (MQ=255)
                                        tGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGa                                                                                                       <  1:29784/149‑1 (MQ=255)
                                          cAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAg                                                                                                      <  1:105824/148‑1 (MQ=255)
                                           aaGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGc                                                                                                     >  2:98015/1‑148 (MQ=255)
                                                  gATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGt                                                                                             >  5:84348/1‑149 (MQ=255)
                                                   aTCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGAATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGtt                                                                                            >  7:299293/1‑149 (MQ=255)
                                                       ggAATCTTCAGCTCTGGCCGGTGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCtt                                                                                        >  6:163209/1‑149 (MQ=255)
                                                       ggAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCtt                                                                                        >  7:218593/1‑149 (MQ=255)
                                                             ttCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                           <  2:161702/124‑1 (MQ=255)
                                                             ttCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                           >  1:161702/1‑124 (MQ=255)
                                                             ttCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCcac                                                                                    >  3:172234/1‑147 (MQ=255)
                                                                aGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACa                                                                                               >  5:114817/1‑133 (MQ=255)
                                                                aGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACa                                                                                               <  6:114817/133‑1 (MQ=255)
                                                                  ctctGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGGATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGcc                                                                             >  1:68935/1‑149 (MQ=255)
                                                                       gCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAc                                                                           <  8:56879/146‑1 (MQ=255)
                                                                         cGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGAct                                                                      >  6:80242/1‑149 (MQ=255)
                                                                           gcgcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGActct                                                                    <  6:84348/149‑1 (MQ=255)
                                                                              ccATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCgg                                                                                                                           >  8:134881/1‑91 (MQ=255)
                                                                              ccATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCgg                                                                                                                           <  7:134881/91‑1 (MQ=255)
                                                                               cATGGCAAGATTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCAt                                                                <  3:67445/149‑1 (MQ=255)
                                                                                 tGGCAAGCATGCTCACCTCGAACAACCAGAATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGACACCGACTCTTCATg                                                               >  4:52164/1‑148 (MQ=255)
                                                                                    cAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATa                                                           <  1:37374/149‑1 (MQ=255)
                                                                                                 cTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAAtcatca                                               >  4:17306/1‑148 (MQ=255)
                                                                                                 cTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAAtcatca                                               >  8:23241/1‑148 (MQ=255)
                                                                                                         acCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCtgt                                      <  5:227896/149‑1 (MQ=255)
                                                                                                             gCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCAc                                  >  4:143213/1‑149 (MQ=255)
                                                                                                               aTCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACgg                                <  2:68935/149‑1 (MQ=255)
                                                                                                                            cGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgt                    <  8:218593/148‑1 (MQ=255)
                                                                                                                            cGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgt                    <  3:17306/148‑1 (MQ=255)
                                                                                                                              ggAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGATATTGAGCAAGTCAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTgg                 <  3:52164/149‑1 (MQ=255)
                                                                                                                              ggAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTgg                 >  6:77634/1‑149 (MQ=255)
                                                                                                                                 gCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCg              >  7:209656/1‑149 (MQ=255)
                                                                                                                                 gCGGCCTTGCGTCGCGAAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgtt                   <  6:296192/144‑1 (MQ=38)
                                                                                                                                 gCGGCCTTGCGTCGCGAAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgtt                   >  5:296192/1‑144 (MQ=38)
                                                                                                                                   ggCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGgc            <  1:98015/149‑1 (MQ=255)
                                                                                                                                   ggCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGgc            >  3:70677/1‑149 (MQ=255)
                                                                                                                                     ccTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGgcgc          >  1:176704/1‑149 (MQ=255)
                                                                                                                                           gcgACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTaa     <  5:163209/148‑1 (MQ=255)
                                                                                                                                            cgACGCAAGGCCGCTCCCACAGGAATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATc   <  8:299293/149‑1 (MQ=255)
                                                                                                                                             gACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATcc  >  7:27632/1‑149 (MQ=255)
                                                                                                                                             gACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATcc  >  5:9078/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CGGGGTGAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCC  >  NC_002947/5886919‑5887208

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: