Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA62.7% ‑4.6 / 45.6 51L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base A (0/32);  total (9/42)
Fisher's exact test for biased strand distribution p-value = 3.04e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-03

AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357745‑358037
                                                                                                                                                      |                                                                                                                                                     
aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑c                                                                                                                                                 <  1:527466/151‑1 (MQ=255)
            tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctctc                                                                                                                                     >  2:1002227/1‑151 (MQ=255)
               gcaaaTCCCAGCCGAGGCCTTGTGTCAGCGGCTCCCGTGA‑TTCATTCCCGCAAA‑CACCCGCGCTTTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTC‑CGCGGCGGCTGTCTGCAAGTTATTTTTCCGCCGCAAAAAAA‑ccgccg                                                                                                                                            <  1:1185898/138‑1 (MQ=17)
                    tcccACCCGCAGCCTGGTGCCGGGGGCTACGGTGT‑TTATTTCCCGCAAA‑AACCCGCGCTTTTA‑CCCAACCCC‑TTGTAAAAAAATGCTTC‑CGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCagtagt                                                                                                                              <  2:1211666/148‑1 (MQ=255)
                     tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                            >  1:851415/1‑151 (MQ=255)
                      ccggCCGGAGCCTTGTGCCGGCGGCTACCTTGA‑TTTTTTCCCCCAAA‑CACCCGCACATTTA‑CCCACCCCC‑TTTGGAAAAATTGCCTC‑CGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTg                                                                                                                             <  2:1015942/146‑1 (MQ=255)
                       cAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                         >  2:998059/1‑151 (MQ=255)
                        aGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gg                                                                                                                        >  2:97630/1‑151 (MQ=255)
                          ccGCAGCCCTGGGCCAGCGCCTACCGTAT‑TTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                                       <  2:642170/150‑1 (MQ=255)
                          ccGCAGACCTGT‑TCAGCCGGAACCCTTATTTCATTCCCGCAAA‑CACCCGCG‑ATTTACCCCAACGCA‑AGGTAC‑AAAATGCCTTGCGCGGGGGCTGTCTGCAACCTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                      <  2:1198551/151‑1 (MQ=255)
                           cgcaacctggtgcCAGCGGCTCCCGTTT‑TTCATTCCCCCCAA‑CACCCGCGCTTTAA‑CCCAACGCC‑TTGTAAAAAAATCCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                      <  1:1209561/141‑1 (MQ=255)
                           cgaagcccTGTGTCAGGGGCTCCCGTGT‑TTCTTTCCCGCCAC‑CACCCGCGCTTTTA‑‑CCAACCCCAATGTGAAAAATTGCCTC‑CGGGGCGGCTTTCTGCTAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGc                                                                                                                     <  2:1139967/145‑1 (MQ=255)
                             cggccgtttgccggCGTCTACCGTAA‑TTCTTTCCCCCCAC‑CACCCGCGTTTTTA‑CCCACCCCC‑ATGAGAAAAAAGGCCTC‑CGGGGGGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGc                                                                                                                     <  1:922202/136‑1 (MQ=255)
                               gACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gg                                                                                                                        >  2:187327/1‑144 (MQ=255)
                                        cggcgcCTACCGTAA‑TTCTTTCCCCCCAA‑CACCCCCCCTTTTA‑TCCAACCCA‑ATGTGCAAAAATGCCTT‑CGGGGCGGCTGTCGGCCCGTTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg                                                                                                          <  2:515583/147‑1 (MQ=255)
                                            gcctccgttGA‑TTCTTTCCCCCCAA‑CACCCCCGCATTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGGTTTCCGGCCAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATa                                                                                                           <  1:115389/136‑1 (MQ=255)
                                            gcctcccGTGA‑TTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCGCC‑ATGTGAAAAAATGCCTC‑CGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCt                                                                                                        <  2:825477/142‑1 (MQ=255)
                                               tACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAGAttt                                                                                                      >  1:584735/1‑146 (MQ=255)
                                                cccGTGT‑TTCTTTCCCCCCAA‑CACCCGCCCTTTTA‑‑CCAACCCCAAGGTGAAAAAATCCCTC‑CGGGGGGGCTGTCGGCAAGTTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc                                                                                                <  1:139969/150‑1 (MQ=255)
                                                cccGTGT‑TTC‑TTTCCCCCAACCACCCGCGCTTTTA‑CCCACCGCC‑TGGTAAAAAAATGCCTC‑CGGGGCGGTTGTCGGCAAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGc                                                                                                  <  1:1061724/148‑1 (MQ=255)
                                                aCCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTaaaacc                                                                                                >  2:556077/1‑147 (MQ=25)
                                                 ccGTGT‑TTTTTTCCCCCCAA‑CACCCGCGCATTTA‑CCCAACGCC‑ATGTGCAAAAATGCCTC‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc                                                                                                <  2:159503/150‑1 (MQ=25)
                                                 ccGTGA‑TTCTTTCCCCCCAA‑CACCCGCGCATTTA‑CCCACCGCC‑TTGTGAAAAAATGCCTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTTCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGc                                                                                                  <  2:788332/148‑1 (MQ=255)
                                                    tGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGgaaaaagttt                                                                                                      >  2:592861/1‑134 (MQ=25)
                                                     gA‑TTCTTTCCCCCCAA‑CACCCGCCTATTTA‑CCCACCCCC‑AGGTCCAAAAAGCCCTC‑CGGGGGGCTTGTCCGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                              <  1:171552/148‑1 (MQ=255)
                                                        ttCTTTCCCGCAAA‑CACCCGCGATTTTA‑CCAAACGCC‑ATGTGAAAAATTGCCTC‑CGGGGGGGTTGTCGGCCACTTATTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                                          <  2:452240/150‑1 (MQ=255)
                                                              ccccccAA‑CACCCCCCCTTTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCTGCAAGCTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                   <  1:556026/147‑1 (MQ=255)
                                                                        acccGCGCATTTA‑CCCACCGCA‑ATGTGAAAAAATCCCTT‑CGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                             <  2:845587/146‑1 (MQ=255)
                                                                                 tttA‑TCCACCGCC‑ATGAAAAAAAGTCCCTG‑CGCGGGGGTTTCCGCCCAGTTTTTTTCCCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                    <  1:786247/148‑1 (MQ=255)
                                                                                       ccACCGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                              <  1:370226/149‑1 (MQ=255)
                                                                                       ccAACGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                               <  2:954857/148‑1 (MQ=255)
                                                                                           cGCC‑ATGTAAAAAAATGCCTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGACAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                            <  2:14748/147‑1 (MQ=255)
                                                                                             cc‑ATGTGAAAAAAGCCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                                       <  2:449158/150‑1 (MQ=255)
                                                                                                      aaaaaTGCTTG‑CGTGGCGGATGACTGCAACAAATTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                                 <  2:462029/148‑1 (MQ=255)
                                                                                                      aaaaaGGCCTG‑CGCGGCGGCTTTCTCCCATCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                                <  1:568480/149‑1 (MQ=255)
                                                                                                         aaaccctcccgcgCCGGTTGTCAGCCGTTTTTTTTCCCCCCCCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                            <  1:278901/142‑1 (MQ=255)
                                                                                                                   gggggggcTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                                      <  2:1066980/142‑1 (MQ=255)
                                                                                                                       cggcTTTCCGCCAGCTATTTCCCCGCCACAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                               <  1:812295/150‑1 (MQ=255)
                                                                                                                              cTGCCAGTTTTTTTCCCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                              <  1:1143432/144‑1 (MQ=255)
                                                                                                                               ggccggttATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                           <  1:1218756/139‑1 (MQ=255)
                                                                                                                                 ccAGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                          <  1:145731/145‑1 (MQ=255)
                                                                                                                                  ccactttttttCCCACCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                    <  1:1002127/145‑1 (MQ=255)
                                                                                                                                     tatttttcCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                     <  1:1230355/127‑1 (MQ=255)
                                                                                                                                      ttttttcCCGCCAAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                   <  1:1064137/145‑1 (MQ=255)
                                                                                                                                      ttttttcCCCCCGAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                  <  1:899604/146‑1 (MQ=255)
                                                                                                                                      tttttccccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                <  1:604086/143‑1 (MQ=255)
                                                                                                                                       tttttcCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                <  2:551856/148‑1 (MQ=255)
                                                                                                                                        ttttcCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa               <  2:25172/149‑1 (MQ=255)
                                                                                                                                        ttttcCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg              <  2:694655/150‑1 (MQ=255)
                                                                                                                                          ttcccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt            <  2:124324/147‑1 (MQ=255)
                                                                                                                                                ccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAAATTTAACCCGCCCTGTTTTGTGGGCATTGCTGCTGGTAACAATATGGGAAGGGGAGAGGGGGGGAACAGACCCAGGACCCCGCCCAATGAGttt           >  2:73755/1‑144 (MQ=255)
                                                                                                                                                   caacaaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAGGAAAAAAAATTTAAG‑CCGTCCTGTTGTGTGGGGAACGCTGCTGGGGACAAAATGGGAAAGGGAGAGGGGGGGATCAGCCCCCCAACCCCGCCAATTGAGTTTGCCCCCATg  >  2:984988/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357745‑358037

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑