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	26 [0]	[0] 25	[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)	24 (0.640)	22/264	0.8	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

AACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGT  >  NC_000913/3815716‑3815993                                                                                                                                               |                                                                                                                                        aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                                                         <  1:521591‑M1/143‑1 (MQ=255)          ttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACGTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                                                <  2:10439‑M1/143‑10 (MQ=255)                    cAAATTCAATACACCGGCGCTGATATGGTTTCATGCCTTCGCCCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                                                      <  1:522987‑M1/143‑20 (MQ=255)                      aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                    <  1:651028‑M1/143‑22 (MQ=255)                         tCAACACACTGGTGCTGATATGGTGTCATGCCTTCGCTCCTCATCTTACTTCTCTACAGGCAAAAAACAGGCGACTCATCAGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcg                                                                                                                 <  2:547561‑M1/143‑25 (MQ=255)                           aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                               <  1:649586‑M1/143‑27 (MQ=255)                            aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcttcccctgcacttccatgatgcaccc                                                                                                              >  1:46021‑M1/1‑116 (MQ=255)                            aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCAGCagtccccagaacttcaatgatgcgccc                                                                                                              >  2:616091‑M1/1‑116 (MQ=255)                                   ggCGCTGATATGGTTTCAGGCCTGCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                       <  1:138890‑M1/143‑35 (MQ=255)                                                 ttCATGCCTTCGCCCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                         <  2:262617‑M1/143‑49 (MQ=255)                                                      gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                                                                    >  1:464436‑M1/1‑90 (MQ=255)                                                                          ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcctcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtccctgcggccgagtct                                                                >  2:444869‑M1/1‑70 (MQ=255)                                                                               tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                           >  1:536633‑M1/1‑65 (MQ=255)                                                                                acagacagAAAACAGGCGACTCATCAGTCGCCCTAAAAACCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                                          <  1:11719‑M1/143‑80 (MQ=255)                                                                                  aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgta                                                        <  1:149731‑M1/143‑82 (MQ=255)                                                                                   gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtactcaacgtat                                                       >  2:92765‑M1/1‑61 (MQ=255)                                                                                      aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagaggcgccaacgtcgccg                                                    >  2:488758‑M1/1‑58 (MQ=255)                                                                                          aaagGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                                <  2:462722‑M1/143‑90 (MQ=255)                                                                                          aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                                <  1:391556‑M1/143‑90 (MQ=255)                                                                                                 gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagatcgcaacc                                         >  1:513323‑M1/1‑47 (MQ=255)                                                                                                  aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaacc                                        >  1:250703‑M1/1‑46 (MQ=255)                                                                                                    tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCccgccgctgcacttcaatgatgcgcccgtcttcggtcatcactgcgttcaggtcggtctctgcggcagagcctgcaacgtattccggatcgcgaaccgc                                      >  2:398463‑M1/1‑44 (MQ=255)                                                                                                               cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcacaccgcttcgccgttca                           >  1:97973‑M1/1‑33 (MQ=255)                                                                                                                  aaaaaTCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaa                        <  1:492306‑M1/143‑114 (MQ=255)                                                                                                                     agtccgtttgcCAGCGCCGCCTGCTGCcgtcccctgcacttcaatgacgcgcccgtcttcggtcatcactccgctcatgtcggtctctgcggcagcgtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattc                     <  2:645305‑M1/138‑117 (MQ=255)                                                                                                                                        ccTTCTGCcgtcccctgcactgcattgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctac  <  1:450859‑M1/143‑136 (MQ=255)                                                                                                                                               |                                                                                                                                        AACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGT  >  NC_000913/3815716‑3815993

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑