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

AACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815716‑3816000                                                                                                                                               |                                                                                                                                               aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAACAACAGATTGACAGCGCCGCCTTCTGC                                                                                                                                                >  2:345715‑M1/1‑143 (MQ=255)      cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCCTCAGTCGCCGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcctcc                                                                                                                                           >  2:362298‑M1/1‑138 (MQ=255)           tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                                                                                                      >  2:752505‑M1/1‑133 (MQ=255)              tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCGCATCAGTCGCCGTAAAAATCAGTTGGGCAGCGCAGCGTTCGGCcgtgcactgcacg                                                                                                                                   >  1:306597‑M1/1‑130 (MQ=255)               acgcgaaaATTCTATAAACTCGCGCTTCTATGGTTTAGTGCCTTTGCTCCTCATCTTCCTTTTATACAGAGAAGAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                                                                                  <  2:459257‑M1/141‑15 (MQ=255)                      aaTTCAATAAGCTGGCGCTGATATGGTTTCAGGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                           <  2:623175‑M1/143‑22 (MQ=255)                       aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                          <  2:314861‑M1/143‑23 (MQ=255)                           aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                      <  1:437784‑M1/143‑27 (MQ=255)                             tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGCCGCATTAAAAATCAGTTTGCCAGCGCCGCCCTCTGCcctcccccgcacctcaacgctgcgcccc                                                                                                                    >  2:714454‑M1/1‑115 (MQ=255)                                 cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                                                                                                                <  2:687139‑M1/143‑33 (MQ=255)                                   gccgctgATACGGTTTCATTCCTTCGCTCCTCATCTTACTTTTCTACAGACCAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGTCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  2:31996‑M1/141‑35 (MQ=255)                                   cgcgctGATATGGTTTCATGCCTTCGCTCCTCATCTTATGTTTCTACAGACAAAAAAAACGCGACTCATAAGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  1:121871‑M1/142‑35 (MQ=255)                                       cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                                                                                                          <  1:616004‑M1/143‑39 (MQ=255)                                          atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                                                       >  1:43910‑M1/1‑102 (MQ=255)                                            atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                                                                                     >  1:733635‑M1/1‑100 (MQ=255)                                              ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTTTGGgggccccagcacttccatgatgcgcccgtcttcggtcatcactac                                                                                                   >  1:643529‑M1/1‑98 (MQ=255)                                                  tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                               <  1:394244‑M1/143‑50 (MQ=255)                                                    aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                             <  1:700725‑M1/143‑52 (MQ=255)                                                        cTTCGCTCCTCATCTGACTTTTCTACAGACTAAAAAAAGGCGCCTCCTCAGTCGCCTTAACAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                                                                         <  2:127272‑M1/143‑56 (MQ=255)                                                             ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccctcaatgatgcgcccgtcttcgctcctcactacgttcatgtcggtccc                                                                                    >  2:18857‑M1/1‑83 (MQ=255)                                                                   ttcTTGCTTCTCTCCAGGCGTAACAAAGGCGACTCATCAGTCGCCTTCAAAATCAGTTTGCCAGCGCCGCCTGCTGCcgtcctctgcacttcaatgatgcgcccggcttcggtcatcactacgttcatgtcggtctctgcggc                                                                              <  2:227079‑M1/142‑67 (MQ=255)                                                                    tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGCCGCCTTACAAATCAGTTTGCCAGCGCCGCCCTCTGCcgtcccctgcacttcaatgatgcgcccgccctcggtcatcaccaccttcatgtccgtccctgcgccc                                                                             >  2:887744‑M1/1‑76 (MQ=255)                                                                      ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                           >  1:581502‑M1/1‑74 (MQ=255)                                                                          ttttCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgttggtctcttcggcagagtct                                                                       >  1:889380‑M1/1‑70 (MQ=255)                                                                          ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCGGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcatgtcggtctctccggcagactct                                                                       >  2:793572‑M1/1‑70 (MQ=255)                                                                               tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                                  >  1:75973‑M1/1‑65 (MQ=255)                                                                                aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcattactacgtgcatgtcggcctctgcggcagaggcttcaacg                                                                 >  2:616099‑M1/1‑64 (MQ=255)                                                                                aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccggcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                                                 >  2:12272‑M1/1‑64 (MQ=255)                                                                                     cAAAAGAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcgcgtcttcgttcaccacttcgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                            <  1:335874‑M1/143‑85 (MQ=255)                                                                                      aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttaaaagtattcc                                                           >  2:709579‑M1/1‑58 (MQ=255)                                                                                             aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgc                                                    >  1:591752‑M1/1‑51 (MQ=255)                                                                                             aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagacttcaacgtatcccatatcca                                                    >  1:232489‑M1/1‑51 (MQ=255)                                                                                                    tcatcaGCCGCCTTAAGACTCGTTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcctgatcgcaaaccgc                                             <  2:408970‑M1/143‑100 (MQ=255)                                                                                                             gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgtt                                    <  1:563289‑M1/143‑109 (MQ=255)                                                                                                                ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtctccaccctattccacaccaccagccgctgcgccgccccc                                 >  2:812250‑M1/1‑32 (MQ=255)                                                                                                                    aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtgttccagatcgcaaaccgcttcgccgttcacactt                             >  2:65896‑M1/1‑28 (MQ=255)                                                                                                                      aTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                           <  1:50370‑M1/143‑118 (MQ=255)                                                                                                                      aTCAGTTGGCCAGCGCCGACATCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                           <  1:581780‑M1/143‑118 (MQ=255)                                                                                                                         ggtgcgccagcTCCGCGCTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatgccagatcgcagaccgcttcgccgttcacaattccgac                        <  2:89948‑M1/138‑121 (MQ=255)                                                                                                                               ccAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgttcatgtcggaatctgcgacagagacatcaacgtattctagtttgcaaaccacaccgccgttcacaattacgacagaacc                  >  1:161785‑M1/1‑17 (MQ=255)                                                                                                                                     ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                   >  2:126802‑M1/1‑11 (MQ=255)                                                                                                                                     ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                   <  1:126802‑M1/120‑110 (MQ=255)                                                                                                                                              gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccc   <  1:137114‑M1/143‑142 (MQ=255)                                                                                                                                              gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggaagagtattcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccc   >  2:133823‑M1/1‑2 (MQ=255)                                                                                                                                               ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacagttccgacagcaactgcggctaccatccct  >  1:799539‑M1/1‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               AACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815716‑3816000

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