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	49 [0]	[0] 49	[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)	48 (1.080)	42/262	0.1	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

ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815719‑3816000                                                                                                                                            |                                                                                                                                               gcctGCTTGCGAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                                                             <  1:589286‑M1/142‑4 (MQ=255)       ttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTACATGCCTTAGCTCATCATCTGACTTTTCTACAGACAAAAAAAAGGCCACTCATCAGTCGCCTTACACATCAGTTTGCCAGCACCGCCTCCTGAcgtccccgg                                                                                                                                       >  2:387304‑M1/1‑134 (MQ=255)            aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                                                                                  <  1:527732‑M1/143‑15 (MQ=255)               cgcAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                               <  1:69793‑M1/143‑18 (MQ=255)                  aaaTTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                                                                                                            <  1:51972‑M1/143‑21 (MQ=255)                  aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatta                                                                                                                            >  1:585206‑M1/1‑123 (MQ=255)                   aaTTCCACAAACTGGCGCTGATATGGTTTCATGCCCTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                           <  1:137608‑M1/143‑22 (MQ=255)                   aaTTCAATAGACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgttccctgcacttcaatgat                                                                                                                           <  2:808015‑M1/143‑22 (MQ=255)                    aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                          <  1:689954‑M1/143‑23 (MQ=255)                           aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAcgtcccctgcacttcaatgatgcgcccgc                                                                                                                   >  2:598093‑M1/1‑114 (MQ=255)                                ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  1:270847‑M1/143‑35 (MQ=255)                                ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCTCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  2:715442‑M1/143‑35 (MQ=255)                                ggCGCTGATATGGTTTCATGCCTTCACTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              >  1:370431‑M1/1‑109 (MQ=255)                                   gcTGATAGGGTTACATGCCGTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                                                           <  1:514936‑M1/143‑38 (MQ=255)                                     tGATATGGTTTCATGCCTTCGCTCCTCGTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                         <  2:127542‑M1/143‑40 (MQ=255)                                       atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                                                       >  1:679926‑M1/1‑102 (MQ=255)                                         atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcact                                                                                                     >  1:650861‑M1/1‑100 (MQ=255)                                         atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcccttcaatgctgcccccgccttcggtcctccct                                                                                                     >  1:391699‑M1/1‑100 (MQ=255)                                            gTTGCATGACGTTGCGCCTCATCTTGCGTTTCTGCAGGCAAGAAGAAGGCGACTCTTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                                                                  <  2:490392‑M1/143‑47 (MQ=255)                                              ttCATGCCATCGCTCCTCATCTTACTGTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                                <  2:626148‑M1/143‑49 (MQ=255)                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcatgt                                                                                           >  1:226942‑M1/1‑90 (MQ=255)                                                     cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTAGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactcacatgatgcgcccgccctcggtaataactacgctcacgacg                                                                                         >  2:673705‑M1/1‑88 (MQ=255)                                                       tCGCTCCTCATCTGAACTTGCTCCACACAAAAAACAGGCGACTCATCAGTCCCGCTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgattcccccgtcttcggtcatcactacgttcatgtcggt                                                                                       <  2:511470‑M1/143‑58 (MQ=255)                                                        cGCTCCTCATCTTACTTTTCTACATACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                                                      <  1:633640‑M1/143‑59 (MQ=255)                                                        cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                                                      >  1:262310‑M1/1‑85 (MQ=255)                                                          ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                                                    >  1:207210‑M1/1‑83 (MQ=255)                                                             ctcATCTTTCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                                                                 <  1:519783‑M1/143‑64 (MQ=255)                                                             ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgtcccttttcgttcatcactacgttcatgtcgttctcctc                                                                                 >  1:395346‑M1/1‑80 (MQ=255)                                                                 tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                             <  1:60485‑M1/143‑68 (MQ=255)                                                                            tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGATTGCCAGCGCCGCCTGCTGCcgtcccctggacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcgatctctgcggcacgggcttcaac                                                                  >  2:60598‑M1/1‑65 (MQ=255)                                                                              cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgt                                                                >  2:271395‑M1/1‑63 (MQ=255)                                                                                   atttttaaggcgacTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                           >  1:493944‑M1/7‑58 (MQ=255)                                                                                   aacaaaaaggcGACTCATCAGTCGCCTTAAAGATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagactcttcaacgtattcc                                                           <  2:743262‑M1/140‑86 (MQ=255)                                                                                   aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                           >  1:527069‑M1/1‑58 (MQ=255)                                                                                    aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcaattcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcca                                                          >  1:539240‑M1/1‑57 (MQ=255)                                                                                     aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcctcccctgcgcgacaatgatgcgcccgtcttcggtcctcactacctccatgtcggtctctgcggcagagccgtccgcgcattccac                                                         >  2:424955‑M1/1‑56 (MQ=255)                                                                                      tttaaggcGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcgtgtctgtctctgcggcacagtcttcaacgtattccaga                                                        >  1:41994‑M1/4‑55 (MQ=255)                                                                                                 tcatcaGTCGCCTTAACAATCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgc                                             <  1:799104‑M1/143‑100 (MQ=255)                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcctcggtcatcactacgttcatgtcggtctctgcggcacagtcttcaacgtattccagagcgcaaaccccttcgccgttcaaaact                             >  1:854432‑M1/1‑28 (MQ=255)                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatcatgcgcccgtcttcggtcatcactacgttcatgtcggcctctgcggcagcgtcttcaccgtattccagaccgcaacccgcctcgcccctcacactt                             >  2:832097‑M1/1‑28 (MQ=255)                                                                                                                     cAGTCTTCCAGCGCCGCCTTCTGCcgtcccctgcacttcaagcatgcgcccgtcttcggtcatcactacgttcatgtcggtccctccggcacagtcttcaacgtactccagctcgcaaaccgcttcgccgttcacaattccga                         <  2:239727‑M1/143‑120 (MQ=255)                                                                                                                       gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagatcgcaaacagcttcgacgttcacactaccgact                       >  1:73288‑M1/1‑22 (MQ=255)                                                                                                                         ttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagatcccccaccccttcgccgttcacaattccgaaaga                     >  1:316820‑M1/1‑20 (MQ=255)                                                                                                                              aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgttcatgtcggtctctgcggcagagtctccaacgtattcaagatagcgacccgcttcgccgctcaccactccgccacaacctg                >  1:595947‑M1/1‑15 (MQ=255)                                                                                                                               gcgcCGCCTTCTGCcgtccccagcactacaacgacgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgc               <  1:43061‑M1/143‑130 (MQ=255)                                                                                                                                cgccgcCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcg              >  1:64798‑M1/1‑13 (MQ=255)                                                                                                                                  ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggc            >  1:882143‑M1/1‑11 (MQ=255)                                                                                                                                  ccgccTCCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcacccctacgttcatgtcggtctctgcggcagaggctgcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggc            <  1:13666‑M1/143‑133 (MQ=255)                                                                                                                                            ccgtccgccgcagttcactgatgcgcccgtcttcggtcatcactccgttcatgtcgttctctgcggcagagtcttctacttattctagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccct  <  2:542853‑M1/143‑143 (MQ=255)                                                                                                                                            |                                                                                                                                               ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815719‑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: ‑