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	49 [0]	[0] 49	[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)	47 (0.870)	40/262	0.3	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

AAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCC  >  NC_000913/3815804‑3816081                                                                                                                                          |                                                                                                                                             tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatccgtcgccttaaaaatcagtttgccagcgccgccttctgcctgccc                                                                                                                                         >  2:410376‑M2/138‑138 (MQ=255)  gcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCgtcctc                                                                                                                                        >  1:632984‑M2/137‑141 (MQ=255)  gcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                                                                        >  1:911721‑M2/137‑143 (MQ=255)   ctggctcagcgcaaagtcaataaactgggcgtgatatggttgcatgccttcgctcctcatctgacgtgcctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                                                                       <  2:1011065‑M2/8‑1 (MQ=255)    ttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTg                                                                                                                                      <  2:201795‑M2/9‑1 (MQ=255)        taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttatacaaagccaaaataggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACt                                                                                                                                  >  1:610085‑M2/131‑143 (MQ=255)            cccaaattcgatccaccggcgctgttctggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatctgtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                              <  1:491916‑M2/17‑1 (MQ=255)                aattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaacaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgcctcctgcCGTGCCCACCacctcactagt                                                                                                                          >  1:534538‑M2/123‑138 (MQ=255)                  ttcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCCTCCCCTGCACTTCAATGATgc                                                                                                                        >  2:598414‑M2/121‑143 (MQ=255)                     aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttacaactctgtgttgtagctccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                     >  1:199352‑M2/118‑143 (MQ=255)                     aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                     >  1:565597‑M2/118‑143 (MQ=255)                          actggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaacaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCAGCACTTGCATGATCGGCCCGTCt                                                                                                                >  2:198581‑M2/113‑143 (MQ=255)                             ggcgccgatacggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                             <  1:939179‑M2/34‑1 (MQ=255)                                gctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                          <  2:331466‑M2/37‑1 (MQ=255)                                    atatggtctcatgccttcgcctctcatcttacgtttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                                                                      <  2:599220‑M2/41‑1 (MQ=255)                                     tatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                     >  1:590234‑M2/102‑143 (MQ=255)                                      atggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                    >  1:135694‑M2/101‑143 (MQ=255)                                        ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTAGGTCCTCcctgc                                                                                                  >  2:627366‑M2/99‑141 (MQ=255)                                            tcatgccttcgctcctcatcttacttttctacagacagaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                              <  2:479393‑M2/49‑1 (MQ=255)                                              atgccttcgctcctccgcttacttttctacagaaaataaaaacgcgccccatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                            <  2:221909‑M2/51‑1 (MQ=255)                                                gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCACCTGCACTTCAATGATGCGCCAGTCTTCGGCCATCACCGCGTTCATGt                                                                                          >  2:968131‑M2/91‑143 (MQ=255)                                                       atcctcatcttacgtttctcacaacaaaaaaacgccgactcatcagtcgccttaaaaatcagttggccagcgccgcctcctgcCGTCCCCTGCACGTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                                                                   <  2:40282‑M2/60‑1 (MQ=255)                                                         cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTGATCACGACGGTAATGTCGGGGTCTg                                                                                 >  1:189703‑M2/82‑143 (MQ=255)                                                         cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                                 <  2:911721‑M2/62‑1 (MQ=255)                                                          ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGc                                                                                >  1:366153‑M2/81‑143 (MQ=255)                                                                    ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCGGAGTCt                                                                      >  1:560882‑M2/71‑143 (MQ=255)                                                                     tttctacgcacaaaaaaaaggcgatgcagcagtcgccttaaaaatcaggttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                     <  2:409683‑M2/74‑1 (MQ=255)                                                                      ttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                    <  1:298500‑M2/75‑1 (MQ=255)                                                                            agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGGGTCTTCAACGTa                                                              >  1:61023‑M2/63‑143 (MQ=255)                                                                             gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctccCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCCTGTCGGCCTCCGCGGCAGAGCCTtcaacccat                                                             >  2:107956‑M2/62‑139 (MQ=255)                                                                               caaaaaaaaggcgactcatcagtcgccttaaaaagcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                           <  1:261553‑M2/84‑1 (MQ=255)                                                                                   tttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGACTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                       >  1:1026653‑M2/56‑143 (MQ=255)                                                                                   aaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa                                                       <  1:589896‑M2/88‑1 (MQ=255)                                                                                    aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCATCGCCTTCCAGAt                                                      >  2:925062‑M2/55‑143 (MQ=255)                                                                                        ggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCa                                                  >  1:833917‑M2/51‑143 (MQ=255)                                                                                        ggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCAGCACTACGTTCATGTCGATCTCTGCGACGGAGTCTTCAACgtagtccatcaggac                                                  >  2:745495‑M2/51‑136 (MQ=255)                                                                                         gcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaa                                                 >  1:984545‑M2/50‑143 (MQ=255)                                                                                               cagcggtcgccttaaaaatcagtctgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt                                           <  2:364290‑M2/100‑1 (MQ=255)                                                                                                atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTAGATCACTGCGGAAGAGTCTTCAACGTATTCCAGAACGCAAAACACtt                                          >  1:1029883‑M2/43‑143 (MQ=255)                                                                                                        ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                  >  1:1007226‑M2/35‑143 (MQ=255)                                                                                                        ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                  <  1:302849‑M2/109‑1 (MQ=255)                                                                                                        ccttaaaaatcagtctgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGTCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                  <  2:842069‑M2/109‑1 (MQ=255)                                                                                                        ccttaaaaatcacctcgccagcgccgccttctgcCCTCCCCTGCACTTCCATGCTGCGCCCGTCTTCGGCCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                  <  2:923763‑M2/109‑1 (MQ=255)                                                                                                          ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCGTCAACGTATTCCAGATCGCACACCGCTTCGCCGTTcac                                >  1:91366‑M2/33‑143 (MQ=255)                                                                                                                atcagtttgccagcgctgccttctgcCGTCACAGGCAATACAAGGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                          <  1:469656‑M2/117‑1 (MQ=255)                                                                                                                     tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCCACGTATTCCAGATCGCAAACCGCTTCGCCGGTCACAATTCCGACAg                     >  1:112257‑M2/22‑143 (MQ=255)                                                                                                                        gccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATCCCGCCAGgaa                  >  2:78991‑M2/19‑143 (MQ=255)                                                                                                                                  ccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATACCCACCGTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTACTCCAGATCGCAAACCGCTTCGGCGTTCACCATTCCGACAGCACCTGCGGCCAc        >  2:652700‑M2/9‑143 (MQ=255)                                                                                                                                        gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACGATTCCGACAGAAACTGCGGCTACCATccc  >  1:617183‑M2/3‑143 (MQ=255)                                                                                                                                          |                                                                                                                                             AAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCC  >  NC_000913/3815804‑3816081

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