Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,995,334 A→T 8.2% intergenic (‑99/+51) PP_3521 ← / ← PP_3522 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,995,3340AT 8.2% 167.5 / 3.6 58intergenic (‑99/+51)PP_3521/PP_3522membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base T (3/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCGCTGGTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCG  >  NC_002947/3995191‑3995476
                                                                                                                                                |                                                                                                                                               
ggCGCTGGTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCg                                                                                                                                            <  8:194407/149‑1 (MQ=255)
ggCGCTGGTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCg                                                                                                                                            <  2:413308/149‑1 (MQ=255)
 gcgcTGGTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCgg                                                                                                                                           <  1:411183/149‑1 (MQ=255)
  cgcTGGTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCGTAGCCGgt                                                                                                                                         <  7:237487/149‑1 (MQ=255)
       gTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCGTAGCCGgtgcca                                                                                                                                     <  3:298254/148‑4 (MQ=255)
            gCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCgg                                                                                                                                 <  4:378843/147‑1 (MQ=255)
              ggTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGttt                                                                                                                             <  7:327834/149‑3 (MQ=255)
               gTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTa                                                                                                                              >  1:102108/1‑148 (MQ=255)
                                    tGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGAcc                                                                                                         >  7:161380/1‑147 (MQ=255)
                                        cTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGcc                                                                                                    <  8:16595/148‑1 (MQ=255)
                                         tCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCgg                                                                                                  >  4:72339/1‑149 (MQ=255)
                                                gtgGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCtttt                                                                                             >  7:175747/1‑147 (MQ=255)
                                                   gTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCa                                                                                           <  6:2405/146‑1 (MQ=255)
                                                   gTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCa                                                                                           >  5:2405/1‑146 (MQ=255)
                                                   gTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAgg                                                                                         <  3:72339/148‑1 (MQ=255)
                                                      aCCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGtttt                                                                                     >  8:55917/1‑149 (MQ=255)
                                                      aCCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGtttt                                                                                     <  8:207629/149‑1 (MQ=255)
                                                         tctGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATTTACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAg                                                                                  >  8:138408/1‑149 (MQ=255)
                                                               gTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGa                                                                             >  3:95529/1‑149 (MQ=255)
                                                                  ggATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCAGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGAcgg                                                                          >  6:384518/1‑149 (MQ=255)
                                                                       gttgttATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACggcggc                                                                      >  5:290739/1‑148 (MQ=255)
                                                                        ttgttATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTa                                                                    >  2:246566/1‑149 (MQ=255)
                                                                         tgttATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTAc                                                                   <  7:406676/149‑1 (MQ=255)
                                                                              taGCGGCATTTAGTCAGGACCGGCACCGATTTACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCtt                                                              <  7:138408/149‑1 (MQ=255)
                                                                              taGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCt                                                               <  1:84871/148‑1 (MQ=255)
                                                                                 cGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGa                                                           >  4:388947/1‑149 (MQ=255)
                                                                                  ggCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGAt                                                          <  5:321526/149‑1 (MQ=255)
                                                                                         aGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGATCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTccac                                                    <  1:52940/148‑1 (MQ=255)
                                                                                          gTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTccac                                                    >  2:391882/1‑147 (MQ=255)
                                                                                           tCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACgg                                                                                                               >  5:110925/1‑87 (MQ=255)
                                                                                           tCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACgg                                                                                                               <  6:110925/87‑1 (MQ=255)
                                                                                           tCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCtt                                                              <  3:218701/136‑1 (MQ=255)
                                                                                           tCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCtt                                                              >  4:218701/1‑136 (MQ=255)
                                                                                             aGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGc                                                <  8:175747/148‑1 (MQ=255)
                                                                                              ggACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCCAGCCGGCTTTTCAGGTTTTGAGCAAGCAGGCGGCTACTGCTTCGCTGTCCAGCAGCAt                                              >  8:400795/1‑149 (MQ=38)
                                                                                               gACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATg                                             >  5:64347/1‑149 (MQ=255)
                                                                                                aCCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGc                                            <  4:95529/149‑1 (MQ=255)
                                                                                                       ccGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAg                                     <  8:161380/149‑1 (MQ=255)
                                                                                                               cAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGTCGTCCAGTGCCAGcc                             >  6:302916/1‑149 (MQ=255)
                                                                                                               cAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCAACTGCTTAGATCTCCACCAGCATGCCGTCCAGTGCCAGcc                             >  1:280720/1‑149 (MQ=255)
                                                                                                               cAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGACGGTGTTCGCGGGTAAACCCGCGCACACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCACCACAAGCATGCCGTCCAGTGCCAGcc                             >  5:124225/1‑149 (MQ=255)
                                                                                                                 aGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCACACGGGCGACCGCGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCGACTGCGTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCgg                           >  2:369086/1‑149 (MQ=255)
                                                                                                                   gaagaaCCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCgggg                         >  1:369780/1‑149 (MQ=255)
                                                                                                                     agaaCCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCa                       >  8:47464/1‑149 (MQ=255)
                                                                                                                            ccGTGAGGATGGCACCGGCTACGCCGGTTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTa                 >  1:247875/1‑149 (MQ=255)
                                                                                                                                 aGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCg           >  1:127305/1‑149 (MQ=255)
                                                                                                                                  ggATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCg                                          >  4:301030/1‑117 (MQ=255)
                                                                                                                                  ggATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCg                                          <  3:301030/117‑1 (MQ=255)
                                                                                                                                  cgcgaacACCGGCGTAGCCGGTGCTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGt          >  7:292146/7‑149 (MQ=255)
                                                                                                                                  cgcgaacACCGGCGTAGCCGGTGCTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGt          >  7:286287/7‑149 (MQ=255)
                                                                                                                                    aTGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGTCGGGGCACCGGTATCAGCGTGc        <  1:391882/149‑1 (MQ=255)
                                                                                                                                     gaacACCGGCGTAGCCGGTGCTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTg         >  4:321687/4‑147 (MQ=255)
                                                                                                                                       gCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAg      <  2:127305/148‑1 (MQ=255)
                                                                                                                                        cACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGt    >  4:231224/1‑149 (MQ=255)
                                                                                                                                         aCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATc               <  8:204314/137‑1 (MQ=255)
                                                                                                                                         aCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATc               >  7:204314/1‑137 (MQ=255)
                                                                                                                                         aCCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTc   >  6:316427/1‑149 (MQ=255)
                                                                                                                                          ccGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCg  <  2:102108/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGCGCTGGTCTTGCGGTCAGTAAGCGTGGTCGTGCCTGTACTCATGGTGTGGTCACCTCTGAGTGGGATTGTTGTTATAGCGGCATTTAGTCAGGACCGGCACCGATATACAAGGAAGAACCACCGTGAGGATGGCACCGGCTACGCCGGTGTTCGCGGGTAAACCCGCGCCCACGGGCGACCGAGCCGGCTTTTCAGGTTTTGAGCAAGACGGCGGCTACTGCTTCGATCTCCACCAGCATGCCGTCCAGTGCCAGCCGGGGCACCGGTATCAGCGTGCAGGTCG  >  NC_002947/3995191‑3995476

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: