Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,476,629	A→G	intergenic (+119/+65)	torI → / ← dsdC	prophage CPS‑53 recombination directionality factor and response regulator inhibitor/DNA‑binding transcriptional dual regulator DsdC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,476,629	0	A	G	100.0%	138.8 / NA	36	intergenic (+119/+65)	torI/dsdC	prophage CPS‑53 recombination directionality factor and response regulator inhibitor/DNA‑binding transcriptional dual regulator DsdC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/21);  total (15/21)

TCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCC  >  NC_000913/2476386‑2476845                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         tCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGta                                                                                                                                                                                                                    <  1:13614/250‑1 (MQ=255)         gACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGtaaattaaat                                                                                                                                                                                                            <  1:336850/250‑1 (MQ=255)                        ggCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCg                                                                                                                                                                                            <  1:103835/251‑1 (MQ=255)                           gACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAAt                                                                                                                                                                                         <  1:156803/251‑1 (MQ=255)                                gCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGc                                                                                                                                                                                    >  1:116182/1‑251 (MQ=255)                                                             cAAGATGGTTATATCGTGACCATTGTGAATTCAAAACTAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAtt                                                                                                                                                        >  2:360790/1‑250 (MQ=255)                                                             cAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCAtt                                                                                                                                                        >  2:281974/1‑250 (MQ=255)                                                              aaGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTAt                                                                                                                                                      <  2:220228/251‑1 (MQ=255)                                                                gATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGa                                                                                                                                                    <  2:252096/251‑1 (MQ=255)                                                                         atCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGttt                                                                                                                                            <  1:85854/250‑1 (MQ=255)                                                                                       gaatgcgaagaTAAGCTCTTAAGCCGAGCCAATGGGTAAAATAGCGGGTAAAATATTTCCCACATCTAAAAAACACCATTACAATCAATCCACTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGATTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTAATTATGAAGTGGTTTTTACTTATTCTCttt                                                                                                                             >  2:303209/11‑250 (MQ=255)                                                                                             aaaaaTAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGcc                                                                                                                          >  2:387446/1‑248 (MQ=255)                                                                                                   aaGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGataata                                                                                                                  <  2:165884/250‑1 (MQ=255)                                                                                                     gCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCAAATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAAc                                                                                                               >  1:358626/1‑251 (MQ=255)                                                                                                           caGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTATATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTc                                                                                                          <  1:317494/249‑1 (MQ=255)                                                                                                                       tGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCAc                                                                                             >  1:435065/1‑251 (MQ=255)                                                                                                                           tAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCCTTCAAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCActgcc                                                                                         >  2:387489/1‑251 (MQ=255)                                                                                                                            aaaaTAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCActgcct                                                                                        <  1:66317/251‑1 (MQ=255)                                                                                                                            aaaaTAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCActgcc                                                                                         <  2:223644/250‑1 (MQ=255)                                                                                                                                             atTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACAACATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATgtggtg                                                                        >  1:167336/1‑250 (MQ=255)                                                                                                                                               ttCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATg                                                                     >  1:99660/1‑251 (MQ=255)                                                                                                                                                  tcACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGta                                                                   >  1:341488/1‑250 (MQ=255)                                                                                                                                                   cacaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaat                                                                 <  2:204149/251‑1 (MQ=255)                                                                                                                                                   cacaTCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGtaat                                                                 >  1:260010/1‑251 (MQ=255)                                                                                                                                                        cTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGc                                                            >  2:455449/1‑251 (MQ=255)                                                                                                                                                          aaaaaaCACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCt                                                           <  1:364693/250‑1 (MQ=255)                                                                                                                                                                 cccATTCCAATCAATCCCCTGCTGCGTCAAGTAGATGTCTGCAGGGGACACCAGCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCAt                                                   <  2:351079/250‑1 (MQ=255)                                                                                                                                                                                ccccTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATc                                     <  1:281974/250‑1 (MQ=255)                                                                                                                                                                                        gctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACg                             <  1:387489/250‑1 (MQ=255)                                                                                                                                                                                         ctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgc                           >  1:203670/1‑251 (MQ=255)                                                                                                                                                                                           tCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgcg                          <  2:134172/250‑1 (MQ=255)                                                                                                                                                                                            cAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGcgcg                        <  2:260010/251‑1 (MQ=255)                                                                                                                                                                                              aGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGcgac                      <  1:319973/251‑1 (MQ=255)                                                                                                                                                                                                               gccacCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATACTGCTGATGGCATTTCACCGTCATATCGCCAAACGGCCCGACGAGCTCACCACTGGCa      >  2:248973/3‑250 (MQ=255)                                                                                                                                                                                                                  cccATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTAATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGcc  <  2:247459/250‑1 (MQ=255)                                                                                                                                                                                                                  acCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGcc  <  1:454250/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         TCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCC  >  NC_000913/2476386‑2476845

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