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	63 [0]	[0] 64	[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)	58 (0.940)	53/264	0.1	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

CAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815800‑3816083                                                                                                                                              |                                                                                                                                               cacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                                                               >  1:1024417‑M2/142‑143 (MQ=255) cacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                                                               >  1:12856‑M2/142‑143 (MQ=255)   cctgctggctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacacacaaaaaaaaggccactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                             <  1:364850‑M2/4‑1 (MQ=255)    ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                                                            >  1:709988‑M2/139‑143 (MQ=255)          gctaagcgcaaattcaataaactggcgctgatatgggttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCAGCa                                                                                                                                      >  1:419744‑M2/133‑143 (MQ=255)            taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGGCCCCCGCCCt                                                                                                                                    >  2:181804‑M2/131‑143 (MQ=255)               gcgcaaattcaacgagcggccgctgatatggtttcatgccttcgctcctcatcttacttttcgacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCa                                                                                                                                 <  1:724566‑M2/16‑1 (MQ=255)                 gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                               <  2:579940‑M2/18‑1 (MQ=255)                    aattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTcaatgct                                                                                                                            >  2:727387‑M2/123‑141 (MQ=255)                        caataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                                                                        <  1:505199‑M2/25‑1 (MQ=255)                         aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatccgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                       >  2:133252‑M2/118‑143 (MQ=255)                          ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttacaaacccgtttcccggccccgccttctgcCCTCCCCTGCCCTTCACTGATCCGccc                                                                                                                      >  1:808063‑M2/117‑143 (MQ=255)                          ataaactggcgctgagatggtttcatgccttcgctcctcatcttacttctctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccggcctctgcCGGCCCCTCCActtcaatgttggctcc                                                                                                                      >  2:284702‑M2/117‑138 (MQ=255)                                 ggccctcatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  2:1180519‑M2/34‑1 (MQ=255)                                  gcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                                                              <  1:11105‑M2/35‑1 (MQ=255)                                    gctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                            >  1:305792‑M2/107‑143 (MQ=255)                                     ctgagagggtgtcgtgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtatgccagcgccgccttctgcCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGtca                                                                                                           <  2:203309‑M2/38‑1 (MQ=255)                                       gatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTCCAATGATGCCCCCGTCTTCggtcacc                                                                                                         >  2:26431‑M2/104‑141 (MQ=255)                                       gatatggtttcatgcctccgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                         <  1:1100457‑M2/40‑1 (MQ=255)                                             gtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCACTGCACTTCAATGATGCGCCCGTCCTAGGCCATGACCacc                                                                                                   >  2:603880‑M2/98‑142 (MQ=255)                                              tttcaggcctgcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                                  <  2:571335‑M2/47‑1 (MQ=255)                                              tagaatgacttcggtcctcatcttactttgcgacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                                  <  2:1091904‑M2/47‑1 (MQ=255)                                               ttcagtccttggcccctcatcttaccgttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                 <  2:406296‑M2/48‑1 (MQ=255)                                                   tgccttcgctcctcatcttacttttctacagacaaaaaaaaggagactcatcagtcgccttaaaaattagtttggcagcgccgccttttgcCGTCCCCTGGCCTTCAATGATGCCGCCGTCTTCCGTCACCACTCCGGTCATg                                                                                             >  1:117914‑M2/92‑143 (MQ=255)                                                    gccttcgctcctcatcgtacttgtctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                            <  1:438772‑M2/53‑1 (MQ=255)                                                     ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAATGATGCGCCCGTCCTCGGTCCTCACTACGTCCAtctc                                                                                           >  2:169732‑M2/90‑143 (MQ=255)                                                     ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                           <  1:231382‑M2/54‑1 (MQ=255)                                                     ccgtccgtcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTACAATGATGCTCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                           <  1:822034‑M2/54‑1 (MQ=255)                                                      cttcgctccgcatcttacttgtctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccatctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                          <  2:1020136‑M2/55‑1 (MQ=255)                                                            tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTACACGGCGgtctcg                                                                                    >  2:331014‑M2/83‑142 (MQ=255)                                                                  tcttacttttcttcagacaaaaaaacggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                              <  2:709988‑M2/67‑1 (MQ=255)                                                                    ttacctttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                            <  1:751121‑M2/69‑1 (MQ=255)                                                                    cttacttttctacagacaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCTGTCTCTGCGGCaga                                                                            >  1:661515‑M2/75‑143 (MQ=255)                                                                        ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                        >  1:943336‑M2/71‑143 (MQ=255)                                                                          tttctacagacaaagaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                      <  2:827953‑M2/75‑1 (MQ=255)                                                                             taccgtcaaaaagaaggcgactcagcagtggccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAGTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                   <  2:184782‑M2/78‑1 (MQ=255)                                                                                 gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGACAGAGTCTTCAACGtac                                                               >  2:692173‑M2/62‑142 (MQ=255)                                                                                  acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTCATGTCGGTCTCTGCGTCAGAGGCTTCAACGTAtt                                                              >  1:311840‑M2/61‑143 (MQ=255)                                                                                   caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTCAACGTATTc                                                             >  1:255351‑M2/60‑143 (MQ=255)                                                                                    aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTCCAACGTCTTcc                                                            >  2:758141‑M2/59‑143 (MQ=255)                                                                                    aaaaaaaaggcgactcatcagtcgccttaaaaatcactttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                            <  2:1081819‑M2/85‑1 (MQ=255)                                                                                      aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCACCGTCTTCAACGTAGTCCAg                                                          >  1:566165‑M2/57‑143 (MQ=255)                                                                                      aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTAGGTCTCTGCGGCCGAGTCTTCCACGTATTCCAg                                                          >  2:162367‑M2/57‑143 (MQ=255)                                                                                        aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCCCTACGTTCATGTCGGTCCCTGCGGCCGAGTCGTCAACGTATTCCAgac                                                        >  2:197070‑M2/55‑142 (MQ=255)                                                                                         aagggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCACGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                       <  1:85582‑M2/90‑1 (MQ=255)                                                                                                  tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGc                                              <  1:369370‑M2/99‑1 (MQ=255)                                                                                                    atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                            >  1:854897‑M2/43‑143 (MQ=255)                                                                                                     ttagtcgccttaaaaatcggtttgtcagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                           <  2:305792‑M2/102‑1 (MQ=255)                                                                                                          cgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGt                                      <  2:1069193‑M2/107‑1 (MQ=255)                                                                                                            ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATCATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGATCTCTGCGGCAGAGTCTTCAACGTATTCCCGATGGCAACACGCATCGGcgtca                                    >  2:535236‑M2/35‑141 (MQ=255)                                                                                                              ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAGACTTCAACGCATTGCAGATCCCAAAccgcaccgccgcacac                                  >  2:1027981‑M2/33‑138 (MQ=255)                                                                                                                   aatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                             <  2:12856‑M2/116‑1 (MQ=255)                                                                                                                         tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTGCACCATTCCGACAg                       >  1:296238‑M2/22‑143 (MQ=255)                                                                                                                              cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCCTCGCCGTGCACAATTCCGCCAGAAacg                  >  2:233571‑M2/17‑142 (MQ=255)                                                                                                                               agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGCCAGAGTCTTCAACGTCTTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCCACATAAACTg                 >  1:645793‑M2/16‑143 (MQ=255)                                                                                                                                tcgccgtcttgtgcCGCTCCCTGCACTTCAGTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCTTCGTTCACAATTCCGACAGAAACTGc                <  2:692961‑M2/129‑1 (MQ=255)                                                                                                                                gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCATCATCGGTCATCACTACGATCATGTCGGTCTCGGCGGCAGACTATTCAACGTATTCCAGAGCGCAAACACCTTCGCAGATCACAAAGCCGACAGAAACAGc                >  2:92051‑M2/15‑143 (MQ=255)                                                                                                                                 cgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg               <  1:577744‑M2/130‑1 (MQ=255)                                                                                                                                     gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCCGTTCACAATTCTCCCACACACTGCGGctc           >  1:207850‑M2/10‑142 (MQ=255)                                                                                                                                       cttctgcCGTCCCCCGCCCTCCACCGATGCGCCCGTCTTCGGTCATCACTACGTCCATGTCGGTCTCTGCGGCAGAGTCATCATCGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGACAGAAACTGCGGCTAcc         <  1:169732‑M2/136‑1 (MQ=255)                                                                                                                                         ccggcCGTCCCCTGCATCGCAATGTTGCGCCCGTCTTCGGTCATCACTACGTGCATGGCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt       <  2:375643‑M2/138‑1 (MQ=255)                                                                                                                                           tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGGTCACAAGTCCGGCCGAAACTGCGGCTACGATcc     >  1:9986‑M2/4‑143 (MQ=255)                                                                                                                                             cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCGTCGCCGTTCACAATTCCGACAGAAACTGCGGGTACCATCCCt   >  1:537177‑M2/2‑143 (MQ=255)                                                                                                                                              cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGTAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt  <  1:727387/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               CAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815800‑3816083

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: ‑