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	45 [0]	[0] 44	[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)	41 (0.900)	36/262	0.2	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                                                                                                                                           |                                                                                                                                              ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcctcagtcgccgtaaaaatcagtttgccagcgccgccttctgcCCTcc                                                                                                                                           >  2:362298‑M2/139‑143 (MQ=255)      tgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGc                                                                                                                                      >  2:752505‑M2/134‑143 (MQ=255)         taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgcgcctcatcttacttttctacagacaaaaaaaaggcggcgcatcagtcgccgtaaaaatcagttgggcagcgcagcgttcggcCgtgcactgcacg                                                                                                                                   >  1:306597‑M2/131‑142 (MQ=255)          acgcgaaaattctataaactcgcgcttctatggtttagtgcctttgctcctcatcttccttttatacagagaagaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                                                                  <  2:459257‑M2/14‑1 (MQ=255)                 aattcaataagctggcgctgatatggtttcaggccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                           <  2:623175‑M2/21‑1 (MQ=255)                  attcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                          <  2:314861‑M2/22‑1 (MQ=255)                      aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                      <  1:437784‑M2/26‑1 (MQ=255)                        taaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagccgcattaaaaatcagtttgccagcgccgccctctgcCCTCCCCCGCACCTCAACGCTGCgcccc                                                                                                                    >  2:714454‑M2/116‑142 (MQ=255)                            ctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCtt                                                                                                                <  2:687139‑M2/32‑1 (MQ=255)                              gccgctgatacggtttcattccttcgctcctcatcttacttttctacagaccaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgtcagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  2:31996‑M2/34‑1 (MQ=255)                              cgcgctgatatggtttcatgccttcgctcctcatcttatgtttctacagacaaaaaaaacgcgactcataagtcgccttaaaaatccgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  1:121871‑M2/34‑1 (MQ=255)                                  ctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                                                          <  1:616004‑M2/38‑1 (MQ=255)                                     atatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                                                                       >  1:43910‑M2/103‑143 (MQ=255)                                       atggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                     >  1:733635‑M2/101‑143 (MQ=255)                                         ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttttggGGGCCCCAGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTAc                                                                                                   >  1:643529‑M2/99‑143 (MQ=255)                                             tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                               <  1:394244‑M2/49‑1 (MQ=255)                                               atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                             <  1:700725‑M2/51‑1 (MQ=255)                                                   cttcgctcctcatctgacttttctacagactaaaaaaaggcgcctcctcagtcgccttaacaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                         <  2:127272‑M2/55‑1 (MQ=255)                                                        ctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCCTCAATGATGCGCCCGTCTTCGCTCCTCACTACGTTCATGTCggtccc                                                                                    >  2:18857‑M2/84‑141 (MQ=255)                                                              ttcttgcttctctccaggcgtaacaaaggcgactcatcagtcgccttcaaaatcagtttgccagcgccgcctgctgcCGTCCTCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                                                              <  2:227079‑M2/66‑1 (MQ=255)                                                               tcttacttttctacagacaaaaaaaaggcgactcatcagccgccttacaaatcagtttgccagcgccgccctctgcCGTCCCCTGCACTTCAATGATGCGCCCGCCCTCGGTCATCACCACCTTCATGTCCGTCCCtgcgccc                                                                             >  2:887744‑M2/77‑140 (MQ=255)                                                                 ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                           >  1:581502‑M2/75‑143 (MQ=255)                                                                     ttttctacagacaaaaaaaaggcgcctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTTGGTCTCTTCGGCAGAGTCt                                                                       >  1:889380‑M2/71‑143 (MQ=255)                                                                     ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccggcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTCCGGCAGActct                                                                       >  2:793572‑M2/71‑143 (MQ=255)                                                                          tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                  >  1:75973‑M2/66‑143 (MQ=255)                                                                           acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATTACTACGTGCATGTCGGCCTCTGCGGCAGAGGCTTCAACg                                                                 >  2:616099‑M2/65‑143 (MQ=255)                                                                           acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCGGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                                                                 >  2:12272‑M2/65‑143 (MQ=255)                                                                                caaaagaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCGCGTCTTCGTTCACCACTTCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                            <  1:335874‑M2/84‑1 (MQ=255)                                                                                 aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTAAAAGTATTcc                                                           >  2:709579‑M2/59‑143 (MQ=255)                                                                                        aggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGc                                                    >  1:591752‑M2/52‑143 (MQ=255)                                                                                        aggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGACTTCAACGTATCCCATatcca                                                    >  1:232489‑M2/52‑141 (MQ=255)                                                                                               tcatcagccgccttaagactcgttttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCTGATCGCAAACCGc                                             <  2:408970‑M2/99‑1 (MQ=255)                                                                                                        gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                    <  1:563289‑M2/108‑1 (MQ=255)                                                                                                           ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTCCACCCTATTCCACACCACCAGCCGCTgcgccgccccc                                 >  2:812250‑M2/33‑138 (MQ=255)                                                                                                               aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTGTTCCAGATCGCAAACCGCTTCGCCGTTCACACtt                             >  2:65896‑M2/29‑143 (MQ=255)                                                                                                                 atcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                           <  1:50370‑M2/117‑1 (MQ=255)                                                                                                                 atcagttggccagcgccgacatctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                           <  1:581780‑M2/117‑1 (MQ=255)                                                                                                                    ggtgcgccagctccgcgctctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATGCCAGATCGCAGACCGCTTCGCCGTTCACAATTCCGAc                        <  2:89948‑M2/120‑1 (MQ=255)                                                                                                                          ccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGAATCTGCGACAGAGACATCAACGTATTCTAGTTTGCAAACCACACCGCCGTTCACAATTACGACagaacc                  >  1:161785‑M2/18‑141 (MQ=255)                                                                                                                                ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                   >  2:126802‑M2/12‑120 (MQ=255)                                                                                                                                ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                   <  1:126802‑M2/109‑1 (MQ=255)                                                                                                                                         gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATccc   <  1:137114‑M2/141‑1 (MQ=255)                                                                                                                                         gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGAAGAGTATTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATccc   >  2:133823‑M2/3‑143 (MQ=255)                                                                                                                                          cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAGTTCCGACAGCAACTGCGGCTACCATCCCt  >  1:799539‑M2/2‑143 (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: ‑