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	29 [0]	[0] 29	[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)	28 (0.740)	24/264	0.6	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

TGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCT  >  NC_000913/3815722‑3815999                                                                                                                                         |                                                                                                                                              tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                         <  2:549038‑M1/143‑7 (MQ=255)   cttgcgAAGCGCAAATGCAATAAGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                                                                                                       <  1:213683‑M1/143‑9 (MQ=255)           gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGAGCCGCCTTCTGGcggccactgcccttca                                                                                                                               >  1:115288‑M1/1‑127 (MQ=255)                aaTTCAACAATCCGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                          <  2:53450‑M1/143‑22 (MQ=255)                aaTGCAATCAGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                          <  2:78621‑M1/143‑22 (MQ=255)                  ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgatgc                                                                                                                        >  1:390710‑M1/1‑120 (MQ=255)                    cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                      <  1:375182‑M1/143‑26 (MQ=255)                    cAAGAAACTGGCGAGAATATGGTTTCATGCCTTCGGTCCTCATCATACTTTACTAAAGACAAAAAAAAGGCGACTCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                      <  2:137997‑M1/143‑26 (MQ=255)                     aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTGGTGTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                     >  1:172457‑M1/1‑117 (MQ=255)                      aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGACGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacctcaatgatgcgcac                                                                                                                    >  1:138187‑M1/1‑116 (MQ=255)                                   gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgatgcgcccgtcttcggtcatc                                                                                                       >  1:67321‑M1/1‑103 (MQ=255)                                              aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                            <  1:628626‑M1/143‑52 (MQ=255)                                                gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTCAAAATCAGTTTGCCAGAGCCGCCTCCCGCcgttcccggcccttcaatgatgcgcccgtcttcggtcatcactacgttcttgt                                                                                          >  1:270888‑M1/1‑90 (MQ=255)                                                  ctgcgctTCTCATCTGACTTGACGACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAGATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                                                                        <  2:448391‑M1/140‑56 (MQ=255)                                                      gCTCCTCATCTGACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                    <  2:15845‑M1/143‑60 (MQ=255)                                                          gtcaggttactgttCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccttcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                                                                <  1:402735‑M1/137‑64 (MQ=255)                                                                ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                          <  1:10474‑M1/143‑70 (MQ=255)                                                                    ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtcgctgcggcagagtct                                                                      >  1:97806‑M1/1‑70 (MQ=255)                                                                      ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactgcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcacagtcttc                                                                    >  2:482518‑M1/1‑68 (MQ=255)                                                                         tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                                 >  1:267956‑M1/1‑65 (MQ=255)                                                                              aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatt                                                            >  1:363226‑M1/1‑60 (MQ=255)                                                                                atttttaaggcgacTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                          >  2:594293‑M1/7‑58 (MQ=255)                                                                                    aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                                      <  1:53096‑M1/143‑90 (MQ=255)                                                                                              tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcaagatcgcaaaccgc                                            >  1:555423‑M1/1‑44 (MQ=255)                                                                                              tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaaagtcttcaacgtattccagatcgaaaacggc                                            >  1:426489‑M1/1‑44 (MQ=255)                                                                                                          ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttccc                                >  1:89470‑M1/1‑32 (MQ=255)                                                                                                            aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcctcggccatcactacctccacgtcggtctcttcggcagagccttcaccgtattccagctctcaaaccgcttcgcccttcccaa                              >  2:712222‑M1/1‑30 (MQ=255)                                                                                                              aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcctgccgcgcacggcggcagagtcgtcaacgtatgccacatcgcaaaccgctccgccgtccacaatt                            >  1:336266‑M1/1‑28 (MQ=255)                                                                                                                                        gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgatctctgcggcagagtcttcaacgtattcccaatcgcacaccgcttcgccgttcccaattccgacagaaactgcggctagcagccc  >  2:406248‑M1/1‑2 (MQ=255)                                                                                                                                         |                                                                                                                                              TGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCT  >  NC_000913/3815722‑3815999

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