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	43 [0]	[0] 43	[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)	40 (0.840)	33/264	0.5	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

CTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815721‑3816000                                                                                                                                          |                                                                                                                                               cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                                                           >  1:834247‑M1/1‑138 (MQ=255)  tgcttgctAATCGCAAATTCAATAAACTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                          <  1:146940‑M1/143‑7 (MQ=255)                aaaCTCACTAAACTCGCGCTGATATGCATTCATGCCGTCGCTCCTCATCTTACTTTTCTACATACAATACACAGGCCACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                                                                                                            <  2:688767‑M1/143‑21 (MQ=255)                  aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                          <  2:720357‑M1/143‑23 (MQ=255)                     cAATAAACTGGCGCTGATATGGTTTCATCCCTGCGCTCCTCATCACACAAGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                       <  1:514622‑M1/143‑26 (MQ=255)                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGACAAAAACAACGCGCCTCCTCACTCGCCTTCAAAATCCGTTTGCCAGCGCCCCCTTCTGCcgtgccctgcacttcaatcatcgccc                                                                                                                      >  1:862044‑M1/1‑117 (MQ=255)                               gagctgATATGGTTTCATGCCTTCGCTCCGCAACTCACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                                                             <  1:350568‑M1/141‑36 (MQ=255)                                            ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                                <  1:75178‑M1/143‑49 (MQ=255)                                            gtcctgcccTCGCGCCTCATCTTACTTTTCTACAGACAAACAAAAGGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                                <  2:20838‑M1/139‑49 (MQ=255)                                             tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                               <  1:601800‑M1/143‑50 (MQ=255)                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggacataactacgttcatgtc                                                                                          >  1:227736‑M1/1‑89 (MQ=255)                                                    ttCGCTCCTCATCTTGCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttccatcaagcccccctctccagtcctgactaccatcatgtccc                                                                                        >  2:768816‑M1/1‑87 (MQ=255)                                                    ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtccatgtcgg                                                                                        >  1:766831‑M1/1‑87 (MQ=255)                                                                ctttctttTCTACTAACAAAAAAAAGTCTACTCTTCAGTCGCCTTAAAGATCAGTTTGCCAGCGCCGCCTGCTGCcgtcccctgcacgccaatgatgcgcccgtcttcgggcatcactacgttcatgtcggtctctgcggcag                                                                            <  1:266558‑M1/143‑69 (MQ=255)                                                                  tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatatcggtctctgcggcagag                                                                          >  2:342397‑M1/1‑73 (MQ=255)                                                                     ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGGcgacccctgcacttcaatgctgcgcacgtcttcggtcatcaatacgctcatgacggtctctgaggcagagtct                                                                       >  2:444564‑M1/1‑70 (MQ=255)                                                                      tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgaactccgccgcatagtctt                                                                      >  1:120837‑M1/1‑69 (MQ=255)                                                                       ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAGAAATCAGTTTGCCAGCGCCGCCTTCTGCagtcccctgcacttcaatgatgcgcccggcttcggtcctcaccacgttcatgtcggtatctgcgtgagagtcttc                                                                     >  1:145224‑M1/1‑68 (MQ=255)                                                                       ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                                                                     >  1:805541‑M1/1‑68 (MQ=255)                                                                        tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctcca                                                                    >  2:262187‑M1/1‑67 (MQ=255)                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaa                                                                   <  1:333549‑M1/143‑78 (MQ=255)                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctccggcagagtcttcaa                                                                   <  2:771611‑M1/143‑78 (MQ=255)                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcatgtcggtctctgcggcagagccttcaacgta                                                               >  2:204618‑M1/1‑62 (MQ=255)                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgta                                                               >  1:166536‑M1/1‑62 (MQ=255)                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgccggtctctgcggcagagccttcaacgtct                                                              >  2:445612‑M1/1‑61 (MQ=255)                                                                               aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagccttcaacgtatt                                                             >  1:862234‑M1/1‑60 (MQ=255)                                                                               aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgaccaggtcggcctctgcggcagagtcttcaacgtatt                                                             >  1:859553‑M1/1‑60 (MQ=255)                                                                                 aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtaaagagcggtctctgcggcagagccttcaccgcatccc                                                           >  2:641605‑M1/1‑58 (MQ=255)                                                                                  aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcca                                                          <  1:518752‑M1/143‑87 (MQ=255)                                                                                  aaaaaaaGGCGACTCATCAGTAGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcctcggtcatcactacgctcatgtcggtctctgcggcagagtatacaacgtattcca                                                          >  2:563000‑M1/1‑57 (MQ=255)                                                                                   aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgcgttcaacgtgtgccag                                                         >  1:352799‑M1/1‑56 (MQ=255)                                                                                         ggCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgcggcagagcattccacgtcttccagagcgca                                                   >  1:93944‑M1/1‑50 (MQ=255)                                                                                           cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcttcaccgtattccatatcgcaaa                                                 >  2:727206‑M1/1‑48 (MQ=255)                                                                                                 atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatcccagaccgcaaaccgccc                                           >  2:608181‑M1/1‑42 (MQ=255)                                                                                                   caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcg                                         <  1:138586‑M1/143‑104 (MQ=255)                                                                                                   caGTCGCCTTAAAAACCAGTTTGCCAGCGCCGCCTTCCGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggacgagtcttcaatgtattccagatcgcaaaccgcttcg                                         <  1:473032‑M1/143‑104 (MQ=255)                                                                                                     gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccatgcacttcaatgatgcgcccgtcttccgtcatcactccgttcatgccggtctctgcggcagagtctgcaacgtattccagaacgcaaaccgcttcgcc                                       >  2:31299‑M1/1‑38 (MQ=255)                                                                                                                aaTCAGTTGGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtcgctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattc                            <  1:9201‑M1/143‑117 (MQ=255)                                                                                                                    aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac                        <  2:341148‑M1/143‑121 (MQ=255)                                                                                                                                ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccggcagaaactgcggc            >  1:873635‑M1/1‑11 (MQ=255)                                                                                                                                  gccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggcta          <  2:166536‑M1/143‑135 (MQ=255)                                                                                                                                        tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatcc    <  1:400303‑M1/143‑141 (MQ=255)                                                                                                                                          ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccct  <  1:60087‑M1/143‑143 (MQ=255)                                                                                                                                          |                                                                                                                                               CTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815721‑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: ‑