Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,818 A→G K168K (AAA→AAG)  nadR → DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,8180AG91.4% 93.7 / ‑1.7 36K168K (AAA→AAG) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (17/15);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATC  >  NC_000913/4627587‑4628048
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
tGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGGTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGgacctg                                                                                                                                                                                                                        <  1:629942/250‑1 (MQ=255)
                   aTGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCTATATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGATCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGTGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGaagaag                                                                                                                                                                                                   <  1:51231/250‑1 (MQ=255)
                               gATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGAATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGCAGAAGCCTATGCGCcac                                                                                                                                                                                       >  2:202391/1‑250 (MQ=255)
                                           cgcgcaTTGTTTGACGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTCTTCCAAACTTTTACATATCAGCAACATATTCGCATTCATGCTTTCAACGAAGAGGGAATGGAGCCGTATCCGAACGGCTGGGTTGTGTGGAACAACGGCATCACAAAGTTTATGGCTGCAAAAGGGATTCGGCCGGGCCTG‑AT‑AAAACACTCGGAAGAAGCCGATGACCCCCAGTATATggcac                                                                                                                                                                           >  2:426006/1‑247 (MQ=255)
                                                                         cagcagCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCGAATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                            >  2:362649/1‑251 (MQ=255)
                                                                                          ggATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGa                                                                                                                           <  1:188639/251‑1 (MQ=255)
                                                                                             tcgtctgctctcgttagtacaaccTTTTAAATTTCAGACAAATATTCGCTATCATGCTTTAAACGAAGAGGGAATGGAGCATTATTCGGACGGCTGGGATGGGTGGAGCAACTGCATCAAAAAGTTTATGTCTGTAAAGGGTATTGAGCCGGACGTT‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCCCAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGt                                                                                                                         <  2:45637/227‑1 (MQ=255)
                                                                                                   gcgctGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCg                                                                                                                  >  1:512062/5‑251 (MQ=255)
                                                                                                        ggTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcg                                                                                                             <  1:480301/251‑1 (MQ=255)
                                                                                                           tATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                                                                                                            >  1:548739/1‑249 (MQ=255)
                                                                                                            attacAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAg                                                                                                          <  1:572529/246‑1 (MQ=255)
                                                                                                            attacAAACTTTTAAAGATCAGAAAAATAGTCGCTTTCATGATTTCAACGAAGAGGGCATGGAGCCGTTTCCGGTCTGCGGGGATGTGGGGAGCAACGGCAGCAAAAAGTTTATGGTTGAAAAGGATATGCAGCCCGCCCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCGGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAg                                                                                                          <  2:672892/246‑1 (MQ=255)
                                                                                                                cAAACTTTTAAATATCAGAAAAATACTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGAGGTGCGCAGATCcg                                                                                                     >  2:667311/1‑251 (MQ=255)
                                                                                                                cAAACTTTTAAAAATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGCTTATGGATGAAAAGGGGATTCAGACGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTAGCTCCGAAACGTACCTTTATGAGTATAAGCGGTGCGCAGATccc                                                                                                     >  2:312347/1‑250 (MQ=255)
                                                                                                                         aaaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGAACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGGCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCCGATCCGCGAAAAccc                                                                                            >  1:76345/1‑251 (MQ=255)
                                                                                                                           atatCAGAAAAATATTCGCATTCAGGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGACTGGGATGTGTGGCGCAACGGCATCAAAAAGTCTATGGCTGACAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGGGACACAGTATATGGAACATATGGGGATCGCTACGGTGCTGGTCGAACGTAAACGTACCTTTATGAGTCTAAGGGGTGCGCAGATCCGCGCAAACCCGt                                                                                          >  2:527331/1‑251 (MQ=255)
                                                                                                                               cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCg                                                                                      <  1:295456/251‑1 (MQ=255)
                                                                                                                                  aaaaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCtc                                                                                   >  1:561247/1‑250 (MQ=255)
                                                                                                                                  aaaaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAggg                                                                                                                                                                                                                                       >  2:274119/1‑105 (MQ=255)
                                                                                                                                  aaaaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAggg                                                                                                                                                                                                                                       <  1:274119/105‑1 (MQ=255)
                                                                                                                                   aaaaTAATCGCAGTCATGCTTTCAACGAAGAGGGCATGGATCCGTAGCCGCACGGCTGGGATGTGTGGAGCAACGGCACCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTAc                                                                                  <  1:618578/251‑1 (MQ=255)
                                                                                                                                     aaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg                                                                                >  1:475808/1‑251 (MQ=255)
                                                                                                                                      atatTCGCATTCATGCTTTCAACGAAGAGGGGATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAGAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACTC‑CTCGGAAGAATCCGATGCGCCACAGTATATTGAACATCTGGGGATCGTTACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTAGCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTgg                                                                               <  1:384599/251‑1 (MQ=255)
                                                                                                                                       tatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTcc                                                                                       <  1:100543/242‑1 (MQ=255)
                                                                                                                                       tatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTggg                                                                              >  1:175525/1‑251 (MQ=255)
                                                                                                                                       tatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAAAATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCCCAGTTCCGCGAAAACCCGTTcc                                                                                       >  2:100543/1‑242 (MQ=255)
                                                                                                                                            gCATTCATGCTTACAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCTTTCCGCTACTGGGAatat                                                                         <  1:426006/251‑1 (MQ=255)
                                                                                                                                                   tGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCAGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGAACCTTTATGATTATCAGCCGTGCACAGATCCGCGAAAACCCGGTcctattatggggatataatcct                                                                   >  2:355115/1‑230 (MQ=255)
                                                                                                                                                        tCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTAcgga                                                              <  1:544492/250‑4 (MQ=255)
                                                                                                                                                         caaagaggagGGCATGTAGCTGTTTCCGCACGGGTGGGATGTGGTGTGCATCGTCCTCAAAAATGTTATGGCTGAAAATGGGATTCAGACGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAGCCGCGAAAACCCGTTCCGCGACTGGGAATGTATTccgacggaa                                                             <  2:92499/243‑8 (MQ=255)
                                                                                                                                                                         gAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTtgt                                            >  1:205443/1‑251 (MQ=255)
                                                                                                                                                                                            gggATGTCTGGATCAACGGCATCAAAAAGGTTATCACAGTAATGGGGATTCATCCGGA‑CTGGATCTACGC‑CTCGGAATAGTCCGATTCGCCACAGGATATTGAACATGTGGGTACCGATACGGTGCTGGTCTATCCGAAACGGATCTCTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGGTCCGCTACTGGGAATATTATCCTACGGAAGTGAAGCCGTGCCTTGTACGTACCGCGTCGATccgt                         <  2:99347/251‑3 (MQ=255)
                                                                                                                                                                                                   gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTggtgtcg                  >  1:111290/1‑246 (MQ=255)
                                                                                                                                                                                                             ggCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCgg        <  2:205443/251‑1 (MQ=255)
                                                                                                                                                                                                              gCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGt       >  1:140523/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaaaaGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATc  >  2:210179/1‑251 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
TGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTG‑ATCTACAC‑CTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATC  >  NC_000913/4627587‑4628048

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