Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,815,859	Δ82 bp		[rph]–[rph]	[rph], [rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3815859	3815940	82	62 [0]	[0] 63	[rph]–[rph]	[rph],[rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	0 (0.000)	56 (0.850)	45/262	0.4	100%	pseudogene (24/48 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
	?	NC_000913	3815941 =	0 (0.000)					pseudogene (609/669 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGT  >  NC_000913/3815719‑3815997                                                                                                                                            |                                                                                                                                            acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATAGGGCTCCTTAAACATCAGTCGGCCAGCGCCGCCTTCTGCcgt                                                                                                                                          >  1:474788‑M1/1‑140 (MQ=255) acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCCCCTCATCCTACCTTCCTACAGACAAAACAAACCCGACCCACCAGACGCCTTAAAAATCAGTTTACCAGAGCCGCCTTCTGCcgt                                                                                                                                          >  1:275290‑M1/1‑140 (MQ=255)     gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                                                                      >  1:817926‑M1/1‑136 (MQ=255)      cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                                                                                                     <  1:1255191‑M1/143‑9 (MQ=255)         cttaaGCGCAAATTCAATATGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                                                                  <  1:379504‑M1/141‑12 (MQ=255)             cgcgcaAATTCAATAAACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                                                                                                              <  2:396691‑M1/142‑16 (MQ=255)               cgcAAATTCTATAAACTGGCGCAGAGATGGAGTCATGCCTACGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                            <  1:120074‑M1/143‑18 (MQ=255)                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcaat                                                                                                                           >  1:126101‑M1/1‑125 (MQ=255)                     ttCAATAAACTGGCGCTGATATGGTTACATGCCTTCGCTCCTCATCTTACCTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgc                                                                                                                      <  1:396453‑M1/143‑24 (MQ=255)                       cAATAACCGGGCGCTGATATGCTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                    <  2:362337‑M1/143‑26 (MQ=255)                       cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                    <  1:837886‑M1/143‑26 (MQ=255)                        aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcactgatgcgcc                                                                                                                   >  1:1121170‑M1/1‑117 (MQ=255)                                gtggctgaGAGGTTTTCATGCCTTCGCTCGTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                           <  2:43283‑M1/140‑35 (MQ=255)                                ggCGCTGATATGGTTTCAGGCCGTCGCGCGTCATCTTACTTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                           <  2:1205100‑M1/143‑35 (MQ=255)                                   gcTGATATGGTTTCATGCCTTCACTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                                                        >  1:189611‑M1/1‑106 (MQ=255)                                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTTCAGACGAAAAAAAGCCCACTCATCTGCCGCTTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacctcaatgatgcgcccgtcttcggtcat                                                                                                      <  1:48218‑M1/143‑40 (MQ=255)                                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                      >  1:1004857‑M1/1‑104 (MQ=255)                                       acatggTTTCATGCCTTCGCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                                                    <  1:747320‑M1/141‑42 (MQ=255)                                         atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAATTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                                                                                  >  1:1017631‑M1/1‑100 (MQ=255)                                            gTTTCATGCCTTCGCTCCTCATCCTCCTATCCTCTACACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                                                               <  1:651262‑M1/143‑47 (MQ=255)                                              ttCATGCCTTCGCTCCTCATCTTATTTTTCTACAGTCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                             <  1:341456‑M1/143‑49 (MQ=255)                                                 ctgCCTTCGCTCCTCATCTGGCGTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                                                         <  1:929160‑M1/142‑53 (MQ=255)                                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                          >  2:403455‑M1/1‑92 (MQ=255)                                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttact                                                                                          >  2:651856‑M1/1‑92 (MQ=255)                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcaatgatgcgcccgccttcggtcatcactacgttcatgt                                                                                        >  1:626988‑M1/1‑90 (MQ=255)                                                     cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcgtcggtcatcactccgttcatgccg                                                                                      >  1:416319‑M1/1‑88 (MQ=255)                                                     cTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                                                                      <  1:432512‑M1/143‑56 (MQ=255)                                                      ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                                                                     >  1:915060‑M1/1‑87 (MQ=255)                                                         gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                  >  2:886857‑M1/1‑84 (MQ=255)                                                             ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccagcaatgaaatgattcgcgcgtatgagcgcagcactacgttaatgtcggtctctgc                                                                              >  2:967887‑M1/1‑80 (MQ=255)                                                                 tCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCGTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                          <  2:160831‑M1/143‑68 (MQ=255)                                                                 tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                          >  1:1193400‑M1/1‑76 (MQ=255)                                                                   ttACTCCTCTACAGACAAAAAAAAGGCGCCTCATCCGTAGCCTTAAAAATCAGTCTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                        <  1:652844‑M1/143‑70 (MQ=255)                                                                       ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccttcttcggtcatcaatacgttcatgtaaataactgcagcagagtgt                                                                    >  1:563269‑M1/1‑70 (MQ=255)                                                                       ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                                                    >  1:267617‑M1/1‑70 (MQ=255)                                                                       ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTACTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaagtccatatcggtctctgctgccgagtcc                                                                    >  1:1167337‑M1/1‑70 (MQ=255)                                                                        tctccacagacAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                   <  1:254733‑M1/138‑75 (MQ=255)                                                                             aCAGACAAAAAAAGGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaaggatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                                              <  2:650641‑M1/143‑80 (MQ=255)                                                                                gACAAAAAAAAGGCGACTACTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgccttcggtcatctcctcgttcttgtctgtctctgcggcagagtcttcaacgtat                                                           <  1:54674‑M1/143‑83 (MQ=255)                                                                                 aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctccggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacctagt                                                          >  2:980333‑M1/1‑60 (MQ=255)                                                                                 aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccagaacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgacgccgcgtcttcgacgtact                                                          >  2:246336‑M1/1‑60 (MQ=255)                                                                                        aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                                                   <  2:189611‑M1/143‑91 (MQ=255)                                                                                        aaaGGCGACTCATCAGTCGCCTTAAAAAGCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                                                   <  1:403422‑M1/143‑91 (MQ=255)                                                                                           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgca                                                <  2:1274562‑M1/143‑94 (MQ=255)                                                                                              gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgttcctatccatctctgcggcagagtcttccaccaattccagatcgcacac                                             >  2:293196‑M1/1‑47 (MQ=255)                                                                                              gACTCATCAGTCGCCTGAAACTTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaac                                             <  2:979185‑M1/143‑97 (MQ=255)                                                                                                cTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccg                                           <  2:758092‑M1/143‑99 (MQ=255)                                                                                                   atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                        >  2:61368‑M1/1‑42 (MQ=255)                                                                                                        tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgatgcgcccgtcttcgatcatccctccgtccatgtccgtctctgcggccgagtcttccacgccttccagatcgcaacctccttcgccg                                   >  1:316743‑M1/1‑37 (MQ=255)                                                                                                        tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcatcgtcttccagctcgcaacccgcttcgccg                                   >  1:272438‑M1/1‑37 (MQ=255)                                                                                                           ccTTAAAAAGCGGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactgccatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                <  1:530474‑M1/143‑110 (MQ=255)                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcacgtcttcggccatccctccgctcatctcggtctctgcggcagactcttcaacgtataccagatcgcaacccccctcgccgttcacaatt                          >  2:162188‑M1/1‑28 (MQ=255)                                                                                                                         tcgacagctCCACCTTCTGCcttcctctgcacttcaatgatgcgcctgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacaga                  <  2:1003407‑M1/139‑124 (MQ=255)                                                                                                                               gcgcCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcaccattccgccagcaactgc            >  1:896905‑M1/1‑14 (MQ=255)                                                                                                                                 gccgccTTCTGCcgtcccctgcacttcactgttgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagagcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg          <  2:633314‑M1/143‑132 (MQ=255)                                                                                                                                 gccgccTTCTGCcgtcccctgcactgcaatgatgcgcccgtcttcggtcatcactacgttcatgtctgtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg          <  2:180265‑M1/143‑132 (MQ=255)                                                                                                                                   cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccactgcgccggtcacaattccgccagagactgcggcc        >  1:1237489‑M1/1‑10 (MQ=255)                                                                                                                                     ctttctGCtgttccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatgccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctac      <  1:157359‑M1/141‑136 (MQ=255)                                                                                                                                     ccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcagcggattggacatcccaaaccgcttcggcgtgcacaattccgacagaaactgcggccac      >  2:1120319‑M1/1‑8 (MQ=255)                                                                                                                                       ttCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccaacagtaaatgcggctaaaa    >  1:156614‑M1/1‑6 (MQ=255)                                                                                                                                         cTGCcgtcccctgcacttcaattatgtctccgtcttctgtcatcacaacgttcatgtcggtctcggcggcagaggcgtcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatc  <  1:651856‑M1/143‑140 (MQ=255)                                                                                                                                         cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatc  <  1:403455‑M1/143‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGT  >  NC_000913/3815719‑3815997

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