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	36 [0]	[0] 38	[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)	35 (0.880)	30/264	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

ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815799‑3816083                                                                                                                                               |                                                                                                                                               tcacctgcttgctaagcgcaaccgcactaaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagaccaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                                <  1:45684‑M2/1‑1 (MQ=255)          tgctttgcgcaaatccaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGc                                                                                                                                       <  2:577924‑M2/10‑1 (MQ=255)            ctaagcgcaagttcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttccgcCGTCCCCcgccc                                                                                                                                     >  2:478864‑M2/132‑141 (MQ=255)                  gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgagtcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcTGTCTCCTGTACTTCAat                                                                                                                               >  1:76384‑M2/126‑143 (MQ=255)                     aattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                            <  1:240865‑M2/21‑1 (MQ=255)                        gcaagaacctggcgctgatatggtttcatgccgtcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcg                                                                                                                         <  2:145682‑M2/24‑1 (MQ=255)                                  ggcgctgctatggtttcatgtcttcgctcctcatcttactattctacagacaaaaaaaaggcgactgatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  2:372151‑M2/34‑1 (MQ=255)                                         atatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                                                                        >  1:385992‑M2/103‑143 (MQ=255)                                               tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                                  >  1:671244‑M2/97‑143 (MQ=255)                                                   atgtcttggctcctcatcttacttttctacagacaaaaaaaaggggactcatcagtcgccttaaaaatcagcttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                              <  2:78050‑M2/51‑1 (MQ=255)                                                   atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTAGAACACCACTACGTTCAt                                                                                              >  2:755599‑M2/93‑143 (MQ=255)                                                                   tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                              <  1:620380‑M2/67‑1 (MQ=255)                                                                       acttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                                                                          >  1:755355‑M2/73‑143 (MQ=255)                                                                        cttttctaccgacaaaagctaagtgcttcatcagtcgccttaaaactcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                         <  2:714588‑M2/72‑1 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                        >  1:482799‑M2/71‑143 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCCCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGCCCCAGTCt                                                                        >  2:337272‑M2/71‑143 (MQ=255)                                                                          tttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                       >  1:141494‑M2/70‑143 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTAATTTCATGTCGGTCTCTGCGGCAGAGTCTtcacc                                                                   >  1:26012‑M2/66‑141 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                   >  2:94088‑M2/66‑143 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgacagcaccgccttctacCGTCCTCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTAGGTCTCTGCAGCAGAGCCTTcaag                                                                   >  1:216089‑M2/66‑142 (MQ=255)                                                                                 agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTa                                                                >  1:137414‑M2/63‑143 (MQ=255)                                                                                     aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCCATGATTCGCCTGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCTTCTTCACCGCATTcc                                                            >  1:398794‑M2/59‑143 (MQ=255)                                                                                     aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCGTCACTACGTTCATGCCGGTGTCTGCGGCAGAGTCTTCAACGCATtca                                                            >  1:635696‑M2/59‑142 (MQ=255)                                                                                         aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTCTTCCAGAt                                                        >  1:13928‑M2/55‑143 (MQ=255)                                                                                                      tcaggcgccttaaaaatcagtttcccagcgcccccggctgcCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                           <  1:711871‑M2/102‑1 (MQ=255)                                                                                                            gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                     <  1:232020‑M2/108‑1 (MQ=255)                                                                                                             ccttaaaaatccgcttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                    <  2:663100‑M2/109‑1 (MQ=255)                                                                                                              cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTca                                   <  1:714695‑M2/110‑1 (MQ=255)                                                                                                              cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTca                                   <  1:723019‑M2/110‑1 (MQ=255)                                                                                                               ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                  >  2:520293‑M2/33‑143 (MQ=255)                                                                                                                 aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCACACCGCTTCGCCGGTCACaa                                >  2:109355‑M2/31‑143 (MQ=255)                                                                                                                     atcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                            <  1:369643‑M2/117‑1 (MQ=255)                                                                                                                       cagtttgccagccccgccttctgcCGTCCCTTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                          <  2:515886‑M2/119‑1 (MQ=255)                                                                                                                        agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                         <  2:141494‑M2/120‑1 (MQ=255)                                                                                                                         gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACa                        <  2:776491‑M2/121‑1 (MQ=255)                                                                                                                                              cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTCTTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAGCCAGCCTGCGGCTACCAtcccg   >  2:602215‑M2/2‑142 (MQ=255)                                                                                                                                               cGTCCCCTGCACTTCCATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt  <  1:230712/143‑1 (MQ=255)                                                                                                                                               cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt  <  2:385992/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815799‑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: ‑