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	43 [0]	[0] 43	[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)	40 (0.840)	33/264	0.5	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

AAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815803‑3816082                                                                                                                                           |                                                                                                                                              ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttacaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                                                           >  1:834247‑M2/139‑143 (MQ=255)  tgcttgctaatcgcaaattcaataaacttgcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          <  1:146940‑M2/6‑1 (MQ=255)                aaactcactaaactcgcgctgatatgcattcatgccgtcgctcctcatcttacttttctacatacaatacacaggccactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatga                                                                                                                            <  2:688767‑M2/20‑1 (MQ=255)                  attcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                          <  2:720357‑M2/22‑1 (MQ=255)                     caataaactggcgctgatatggtttcatccctgcgctcctcatcacacaagtctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                                                                       <  1:514622‑M2/25‑1 (MQ=255)                      aataaactggcgctgatatggtttcatgccttcgctcctcctcttacttttctacagacaaaaacaacgcgcctcctcactcgccttcaaaatccgtttgccagcgcccccttctgcCGTGCCCTGCACTTCAATCAT‑CGccc                                                                                                                     >  1:862044‑M2/118‑143 (MQ=255)                               gagctgatatggtttcatgccttcgctccgcaactcacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                                                             <  1:350568‑M2/35‑1 (MQ=255)                                            ttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                <  1:75178‑M2/48‑1 (MQ=255)                                            gtcctgccctcgcgcctcatcttacttttctacagacaaacaaaaggcgactcctcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                <  2:20838‑M2/48‑1 (MQ=255)                                             tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                               <  1:601800‑M2/49‑1 (MQ=255)                                                  ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGACATAACTACGTTCATGTc                                                                                          >  1:227736‑M2/90‑143 (MQ=255)                                                    ttcgctcctcatcttgcttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCCATCAAGCCCCCCTCTCCAGTCCTGACTACCATCATgtccc                                                                                        >  2:768816‑M2/88‑141 (MQ=255)                                                    ttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTCCATGTCgg                                                                                        >  1:766831‑M2/88‑143 (MQ=255)                                                                ctttcttttctactaacaaaaaaaagtctactcttcagtcgccttaaagatcagtttgccagcgccgcctgctgcCGTCCCCTGCACGCCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                            <  1:266558‑M2/68‑1 (MQ=255)                                                                  tacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATATCGGTCTCTGCGGCagag                                                                          >  2:342397‑M2/74‑143 (MQ=255)                                                                     ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttcgggCGACCCCTGCACTTCAATGCTGCGCACGTCTTCGGTCATCAATACGCTCATGACGGTCTCTGAGGCAGAGTCt                                                                       >  2:444564‑M2/71‑143 (MQ=255)                                                                      tttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGAACTCCGCCGCATAGTCtt                                                                      >  1:120837‑M2/70‑143 (MQ=255)                                                                       ttctacagacaaaaaaaaggcgactcatcagtcgccttagaaatcagtttgccagcgccgccttctgcAGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCCTCACCACGTTCATGTCGGTATCTGCGTGAGAGTCTTc                                                                     >  1:145224‑M2/69‑143 (MQ=255)                                                                       ttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                     >  1:805541‑M2/69‑143 (MQ=255)                                                                        tctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTCCa                                                                    >  2:262187‑M2/68‑143 (MQ=255)                                                                         ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                   <  1:333549‑M2/77‑1 (MQ=255)                                                                         ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTCCGGCAGAGTCTTCaa                                                                   <  2:771611‑M2/77‑1 (MQ=255)                                                                             agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTCAACGTa                                                               >  2:204618‑M2/63‑143 (MQ=255)                                                                             agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                               >  1:166536‑M2/63‑143 (MQ=255)                                                                              gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGCAGAGCCTTCAAcgtct                                                              >  2:445612‑M2/62‑141 (MQ=255)                                                                               acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTCAACGTAtt                                                             >  1:862234‑M2/61‑143 (MQ=255)                                                                               acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGACCAGGTCGGCCTCTGCGGCAGAGTCTTCAACGTAtt                                                             >  1:859553‑M2/61‑143 (MQ=255)                                                                                 aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTAAAGAGCGGTCTCTGCGGCAGAGCCTTCACCGCATccc                                                           >  2:641605‑M2/59‑143 (MQ=255)                                                                                  aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                          <  1:518752‑M2/86‑1 (MQ=255)                                                                                  aaaaaaaggcgactcatcagtagccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCGGTCATCACTACGCTCATGTCGGTCTCTGCGGCAGAGTATACAACGTATTCCa                                                          >  2:563000‑M2/58‑143 (MQ=255)                                                                                   aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGCGTTCAACGTGTGCCAg                                                         >  1:352799‑M2/57‑143 (MQ=255)                                                                                         ggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGCATTCCACGTCTTCCAGAGCGCa                                                   >  1:93944‑M2/51‑143 (MQ=255)                                                                                           cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCACCGTATTCCATATCGCaaa                                                 >  2:727206‑M2/49‑143 (MQ=255)                                                                                                 atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATCCCAGACCGCAAACcgccc                                           >  2:608181‑M2/43‑141 (MQ=255)                                                                                                   cagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCg                                         <  1:138586‑M2/103‑1 (MQ=255)                                                                                                   cagtcgccttaaaaaccagtttgccagcgccgccttccgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGACGAGTCTTCAATGTATTCCAGATCGCAAACCGCTTCg                                         <  1:473032‑M2/103‑1 (MQ=255)                                                                                                     gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCATGCACTTCAATGATGCGCCCGTCTTCCGTCATCACTCCGTTCATGCCGGTCTCTGCGGCAGAGTCTGCAACGTATTCCAGAACGCAAACCGCTTCGcc                                       >  2:31299‑M2/39‑143 (MQ=255)                                                                                                                aatcagttggccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                            <  1:9201‑M2/116‑1 (MQ=255)                                                                                                                    agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                        <  2:341148‑M2/120‑1 (MQ=255)                                                                                                                                ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGGCAGAAACTGCGGc            >  1:873635‑M2/12‑143 (MQ=255)                                                                                                                                  gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTa          <  2:166536‑M2/134‑1 (MQ=255)                                                                                                                                        tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATcc    <  1:400303‑M2/140‑1 (MQ=255)                                                                                                                                          cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt  <  1:60087‑M2/142‑1 (MQ=255)                                                                                                                                           |                                                                                                                                              AAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815803‑3816082

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