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	55 [0]	[2] 57	[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)	51 (0.970)	45/264	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

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTC  >  NC_000913/3815717‑3815998                                                                                                                                              |                                                                                                                                             acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                             <  2:151065‑M1/143‑2 (MQ=255) acacCTGCTCGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                             <  2:786690‑M1/143‑2 (MQ=255)  cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                                                            >  1:306966‑M1/1‑141 (MQ=255)       gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTACTGCcgtcccc                                                                                                                                       >  2:95882‑M1/1‑136 (MQ=255)         gtgctaAGCGCAAATCCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                                                     <  1:597497‑M1/142‑10 (MQ=255)                 cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                             >  1:25833‑M1/1‑126 (MQ=255)                 cgcAAATGCAATGANCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                             <  1:832988‑M1/143‑18 (MQ=255)                     aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                         >  2:512338‑M1/1‑122 (MQ=255)                          aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                    <  1:745129‑M1/143‑27 (MQ=255)                           aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAAAGACGACTCATCAGCCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCcgtcccccgacactcaatgctgggccc                                                                                                                   >  2:337068‑M1/1‑116 (MQ=255)                             aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCGTTCTGCcgtccctcgcacttcaacgatgcgccccc                                                                                                                 >  2:72920‑M1/1‑114 (MQ=255)                             aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGAGCAGACGTCAGCcgtcccctgcacttcactgctgcgcccgt                                                                                                                 >  1:686683‑M1/1‑114 (MQ=255)                              aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttccatgctgcgcccgcc                                                                                                                >  2:499731‑M1/1‑113 (MQ=255)                               aCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTACTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtgt                                                                                                               >  1:192156‑M1/1‑112 (MQ=255)                                  ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                            <  1:156873‑M1/143‑35 (MQ=255)                                   gcgcTGCTATCGTCTCATGCCTTCGCTCCTCACCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTCGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                                                           <  2:69376‑M1/143‑36 (MQ=255)                                   gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAACGCCGCCTGCTCCcgccccctgcacatcactcatgcgcccgcctccgg                                                                                                           >  2:587275‑M1/1‑108 (MQ=255)                                     gcTGATATGGTTGCATGCCTCCGCTCCTCGTCTTACTTTCCTACAGACAGAAGAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                                                         <  2:594364‑M1/143‑38 (MQ=255)                                       tGAGATGGTTTCATGCCGTCGCTCCTCATCTTACTTGTCTACAGACAAAAAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                       <  2:70251‑M1/143‑40 (MQ=255)                                             ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgccttcggtcatcactcc                                                                                                 >  2:846217‑M1/1‑98 (MQ=255)                                                ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                              >  1:804854‑M1/1‑95 (MQ=255)                                                 tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                             >  1:601582‑M1/1‑94 (MQ=255)                                                 ccacgccttCGCTCCTCATCTTACTTGCCTACCGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgaggcgcccgtcttcggtcatcactacgttc                                                                                             <  1:116153‑M1/139‑50 (MQ=255)                                                   aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacggtcct                                                                                           >  1:860790‑M1/1‑92 (MQ=255)                                                    tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                                                          >  1:306600‑M1/1‑91 (MQ=255)                                                          cGCTTCTCATCTGACTCTTCTACCGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGTCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                                                    <  2:991278‑M1/143‑59 (MQ=255)                                                           gCTCCTCATCTTACTTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                   <  1:335697‑M1/143‑60 (MQ=255)                                                               ctcACCTTACCTTTCTACCCACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                                                               <  1:351721‑M1/143‑64 (MQ=255)                                                                         ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                                                     >  1:155134‑M1/1‑70 (MQ=255)                                                                              tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcggcttcaac                                                                >  1:381336‑M1/1‑65 (MQ=255)                                                                               aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgggcgagtcttcaacg                                                               >  2:512784‑M1/1‑64 (MQ=255)                                                                                   aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttccacggatt                                                           >  1:621690‑M1/1‑60 (MQ=255)                                                                                    cAAAAGAAAGGCGCGTCTTCGGTCTCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcgtcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                          <  2:45038‑M1/143‑85 (MQ=255)                                                                                    cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTCGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                          <  1:4672‑M1/143‑85 (MQ=255)                                                                                     aaaaaaaaGGCGACTCATCAGTCGCATTAAATATCAGTTTTCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                         <  1:956373‑M1/143‑86 (MQ=255)                                                                                      aaaaaaaGGCGTCTCAGCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGCcgtcccctgcacttccatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcatcgtatttta                                                        <  2:175240‑M1/143‑87 (MQ=255)                                                                                      aaaaaaaGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcca                                                        <  1:607243‑M1/143‑87 (MQ=255)                                                                                       aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcgtcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtagtgcag                                                       >  1:560419‑M1/1‑56 (MQ=255)                                                                                        ataaaggCGACTCATCAGTCGCCTCAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga                                                      <  1:566725‑M1/141‑89 (MQ=255)                                                                                        aaaaaGGCGACTCATCAGTCGCTTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga                                                      <  2:428142‑M1/143‑89 (MQ=255)                                                                                               cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaa                                               <  1:390911‑M1/143‑96 (MQ=255)                                                                                                      tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaataatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                        >  1:663911‑M1/1‑41 (MQ=255)                                                                                                            gCCTTAAAAATCGGTGCGCCAGCGCCGCCTTCTCCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgtt                                  <  2:80588‑M1/143‑109 (MQ=255)                                                                                                              cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcctcgccgttca                                >  2:227276‑M1/1‑33 (MQ=255)                                                                                                              cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcagtcatcactacgctcatgtcggtctctgcggcacagccttcaccgtattccacctcgcacaccccaccgcccttca                                >  2:586448‑M1/1‑33 (MQ=255)                                                                                                                   aaaGCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggccgagccttccacgtattccagatcgcaaaccgcttcgccgtccacactt                           >  2:677874‑M1/1‑28 (MQ=255)                                                                                                                    aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagatcgcaaaccgcttcgccgttcacaattg                          >  1:352576‑M1/1‑27 (MQ=255)                                                                                                                       catgtcgccagcgCCGCCTACTGCcgtcccctgcacatcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcagcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccga                       <  2:544710‑M1/137‑120 (MQ=255)                                                                                                                              ccAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgccagaaac                >  1:254811‑M1/1‑17 (MQ=255)                                                                                                                                aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgaccgaaactg              >  1:647377‑M1/1‑15 (MQ=255)                                                                                                                                    ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacttattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggc          <  1:95882‑M1/143‑133 (MQ=255)                                                                                                                                     cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaacagcggct         >  2:810591‑M1/1‑10 (MQ=255)                                                                                                                                            tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatcc  <  1:512338‑M1/143‑141 (MQ=255)                                                                                                                                            tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatcc  <  1:803530‑M1/143‑141 (MQ=255)                                                                                                                                            tGCcgtccccgtacattcaatgctgcgcctgtctacgggcatcactacgttcttgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatcc  <  2:294932‑M1/143‑141 (MQ=255)                                                                                                                                              |                                                                                                                                             ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTC  >  NC_000913/3815717‑3815998

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