Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,643 T→A 6.6% intergenic (‑90/+18) ybcI ← / ← ybcJ DUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,6430TA 6.6% 211.7 / 6.0 61intergenic (‑90/+18)ybcI/ybcJDUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (2/2);  total (30/31)
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.36e-01

AACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTG  >  NC_000913/556496‑556789
                                                                                                                                                   |                                                                                                                                                  
aacccaacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                   >  2:226945/1‑149 (MQ=255)
     aacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                   <  2:329233/144‑1 (MQ=255)
     aacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                   >  1:329233/1‑144 (MQ=255)
     aacccTAAACCAATGCATAGGGGAACGGAGGGGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGGAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGACGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGc                                                                                                                                              <  4:464334/149‑1 (MQ=255)
       cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCg                                                                                                                                            <  3:531247/149‑1 (MQ=255)
        ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                   <  4:207631/141‑1 (MQ=255)
        ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATa                                                                                                                                                   >  3:207631/1‑141 (MQ=255)
         cTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCggg                                                                                                                                          >  3:655856/1‑149 (MQ=255)
         cTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCggg                                                                                                                                          >  3:333616/1‑149 (MQ=255)
               cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAg                                                                                                                                      >  1:576890/1‑147 (MQ=255)
               cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAg                                                                                                                                      <  2:576890/147‑1 (MQ=255)
               cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggtt                                                                                                                                         >  4:156562/1‑142 (MQ=255)
               cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggtt                                                                                                                                         <  1:272029/144‑3 (MQ=255)
               cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCggtt                                                                                                                                         >  2:272029/1‑142 (MQ=255)
                aaTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGa                                                                                                                                   <  1:288017/149‑1 (MQ=255)
                aaTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGa                                                                                                                                   <  4:704639/149‑1 (MQ=255)
                  tGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAg                                                                                                                                                >  1:813563/1‑134 (MQ=255)
                  tGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTAGCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAg                                                                                                                                                <  2:813563/134‑1 (MQ=255)
                    cATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAg                                                                                                                               <  4:655856/149‑1 (MQ=255)
                      tAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGc                                                                                                                             <  4:333616/149‑1 (MQ=255)
                            aaCGGCGGCGTGCGTAATAGCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCaacaac                                                                                                                       <  3:29718/149‑1 (MQ=255)
                               ggcggcGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTg                                                                                                                    <  4:57444/149‑1 (MQ=255)
                                       gcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCaacaac                                                                                                                       <  1:206529/138‑1 (MQ=255)
                                       gcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCaacaac                                                                                                                       >  2:206529/1‑138 (MQ=255)
                                        cgtAATAACGGTTGGCATGTTGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                                                                           <  2:286461/149‑1 (MQ=255)
                                        cgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                                                                           <  4:542715/149‑1 (MQ=255)
                                            ataaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCt                                                                                                       <  4:215312/149‑1 (MQ=255)
                                              aaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGc                                                                                                     >  3:21945/1‑149 (MQ=255)
                                                   ttGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAAc                                                                                                >  3:152389/1‑149 (MQ=255)
                                                      gCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCa                                                                                             >  4:622006/1‑149 (MQ=255)
                                                       cATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCa                                                                                                                                <  2:141992/113‑1 (MQ=255)
                                                       cATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCa                                                                                                                                >  1:141992/1‑113 (MQ=255)
                                                            tgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACACCCTGTACGATGTGACCGGCAAAACTCACTGTCt                                                                                       <  2:142847/149‑1 (MQ=255)
                                                                         aaaaTGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGAt                                                                          >  2:389687/1‑149 (MQ=255)
                                                                          aaaTGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgac                                                                             >  3:2970/1‑145 (MQ=255)
                                                                          aaaTGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgac                                                                             <  4:2970/145‑1 (MQ=255)
                                                                                     cTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGttt                                                              <  4:152389/149‑1 (MQ=255)
                                                                                            cATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGtt                                                       >  3:346348/1‑149 (MQ=255)
                                                                                                 aTGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGGGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTTACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAc                                                  >  3:532690/1‑149 (MQ=255)
                                                                                                             tAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCa                                                                                             >  4:68422/1‑94 (MQ=255)
                                                                                                             tAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCa                                                                                             <  3:68422/94‑1 (MQ=255)
                                                                                                                   gCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAc                                <  1:226945/149‑1 (MQ=255)
                                                                                                                    cAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAcc                               <  1:389687/149‑1 (MQ=255)
                                                                                                                      aaaGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTg                             <  4:21945/149‑1 (MQ=255)
                                                                                                                            tGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTc                       >  4:20378/1‑149 (MQ=255)
                                                                                                                             gAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                      <  4:346348/149‑1 (MQ=255)
                                                                                                                              aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTg                                                                                                      <  4:324871/68‑1 (MQ=255)
                                                                                                                              aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTg                                                                                                      <  2:118364/68‑1 (MQ=255)
                                                                                                                              aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTg                                                                                                      >  1:118364/1‑68 (MQ=255)
                                                                                                                              aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTg                                                                                                      >  3:324871/1‑68 (MQ=255)
                                                                                                                              aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCa                                              >  3:559931/1‑124 (MQ=255)
                                                                                                                              aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCa                                              <  4:559931/124‑1 (MQ=255)
                                                                                                                                  gagaACCGGCCTGTTGATACAGGCCGGGAAAGGGATGAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAc                                                  >  4:500929/1‑116 (MQ=255)
                                                                                                                                  gagaACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGAGTTTGCAGCGTTTGCGCGTTTCAAc                                                  <  3:500929/116‑1 (MQ=255)
                                                                                                                                    gaACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcg               >  1:131455/1‑149 (MQ=255)
                                                                                                                                     aaCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcgc              >  1:809604/1‑149 (MQ=255)
                                                                                                                                            cTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTc       >  1:225254/1‑149 (MQ=255)
                                                                                                                                             tGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCg      >  4:739750/1‑149 (MQ=255)
                                                                                                                                                 gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAc                                                  >  2:560784/1‑101 (MQ=255)
                                                                                                                                                 gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAc                                                  <  1:560784/101‑1 (MQ=255)
                                                                                                                                                 gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTg  >  4:642793/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTG  >  NC_000913/556496‑556789

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: