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	37 [0]	[0] 37	[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)	36 (1.010)	32/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

AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCT  >  NC_000913/3815802‑3816073                                                                                                                                            |                                                                                                                                     cctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaagaaaaaggcgactcatcagtggccttaaaaatcagtttgctagcgccgcattttgccgcc                                                                                                                                   >  1:82439‑M2/140‑141 (MQ=255)     cttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttatacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                                                               <  2:443472‑M2/8‑1 (MQ=255)        gctaagctcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                            <  1:577164‑M2/11‑1 (MQ=255)              cgcaaattcaatagactggcgctgatatggtttcatgccttcgctcctcatcttacttaaaaaaaaaaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                      <  2:31547‑M2/17‑1 (MQ=255)                 aaattcaataaactggcgctgatatggtttcatgccttcgcacctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcccctaaaaaatcagtttgccagcgcccccttcggcCGTCTCATGCCCTTCTatga                                                                                                                   >  1:313413‑M2/124‑143 (MQ=255)                        atattctggcgctgatatggtttcatgccttcgctcctcatcttactttcctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                            <  1:517438‑M2/27‑1 (MQ=255)                        ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                                                            >  1:685094‑M2/117‑143 (MQ=255)                              tggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGCTGCGCCCGGCTTc                                                                                                      >  1:296839‑M2/111‑143 (MQ=255)                                 cgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                   <  2:180381‑M2/36‑1 (MQ=255)                                 cgctgatatggtggcatgccttcgctcctcatcttactgttctacagacaaaagaaaggcgacgcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                   <  2:337303‑M2/36‑1 (MQ=255)                                    tgatatggtttcattccttcgctcctcatcttacttttctacagacaagaaaaaggcgactcatcaggcgccttaaaaatctgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                <  1:279686‑M2/39‑1 (MQ=255)                                            tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccataaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAGTGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                        >  2:58410‑M2/97‑143 (MQ=255)                                              tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaactcagtttgccagcgccgccttctgcCGTCCCCTGCTCTTCACTGATGCGCCCCTCTTCGTTCTTCACTACTTTc                                                                                      >  1:100835‑M2/95‑143 (MQ=255)                                               catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCa                                                                                     >  1:562646‑M2/94‑143 (MQ=255)                                                atcctttcgctccttatcctccttttctacagacagaaaactggcgactcaccagtcgccttaaaattcagtttgccggcgcctccttctgcCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                    <  2:77427‑M2/51‑1 (MQ=255)                                                  gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCGGTCTTCGGTCATCACTACGTTCATGt                                                                                  >  1:20839‑M2/91‑143 (MQ=255)                                                  gccctctctcctcatcttatttgtctacatacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                  <  2:519678‑M2/53‑1 (MQ=255)                                                          tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                          >  1:156104‑M2/83‑143 (MQ=255)                                                             tcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACCTTCATGTCGGTCTCTGCg                                                                       >  2:125012‑M2/80‑143 (MQ=255)                                                                tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                    <  1:641128‑M2/67‑1 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCggcagaggat                                                              >  1:122631‑M2/71‑140 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCAGCGGCAGAGGCt                                                              >  2:558129‑M2/71‑143 (MQ=255)                                                                        ttctccagacaaaaaaaaggcgactcgtcagtcgccttaaaacgcagtttgccagcgccgccgttggcCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                            <  1:623526‑M2/75‑1 (MQ=255)                                                                           tacagacaaaaaaaaggcgactcatcagtcgccttcaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTCAAc                                                         >  1:372347‑M2/66‑143 (MQ=255)                                                                           tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                         <  1:566608‑M2/78‑1 (MQ=255)                                                                                acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCCGTCATCACTACGTTCATGTCGCTCTCTGCGTCCGCGACTTCTACGTCtt                                                    >  2:675239‑M2/61‑143 (MQ=255)                                                                                  aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                  >  1:21045‑M2/59‑143 (MQ=255)                                                                                      ttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCcagcc                                              >  2:273148‑M2/55‑141 (MQ=255)                                                                                            cgactcaccagtcgcgtgaaaaaccagtttgccagcgccggcttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                        <  2:619203‑M2/95‑1 (MQ=255)                                                                                                         gccttaaaactctgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGTTCGCAAACCGCTTCGCCGtt                           <  2:24602‑M2/108‑1 (MQ=255)                                                                                                           cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTca                         >  1:671645‑M2/34‑143 (MQ=255)                                                                                                           cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTca                         >  1:652100‑M2/34‑143 (MQ=255)                                                                                                            ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                        >  1:383298‑M2/33‑143 (MQ=255)                                                                                                                   tcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTCCAACGTATTCCAGATCGCAAACCGCTTCGCCGGTCACCGTTCCg                 >  1:642535‑M2/26‑143 (MQ=255)                                                                                                                    cagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCGGTCATCACTACGTTCATGTCGGTGTCTGCGGCAGAGCCTCCCACGTCTCCCAGCTCGCACACCGCCTCGCCCCTCACCATCCCGa                >  2:627995‑M2/25‑143 (MQ=255)                                                                                                                                gccgccttctgcCGGCCCTCTCACTGCACTGGTGCGCCCGTGTCCGGTCATCACTACGTGCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg    <  2:535550‑M2/131‑1 (MQ=255)                                                                                                                                  cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGGCAGCACCTGCGGCt  >  1:73929‑M2/11‑143 (MQ=255)                                                                                                                                            |                                                                                                                                     AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCT  >  NC_000913/3815802‑3816073

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