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)	47 (0.870)	40/262	0.3	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

TGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCT  >  NC_000913/3815722‑3815999                                                                                                                                          |                                                                                                                                              tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCctgccc                                                                                                                                         >  2:410376‑M1/1‑137 (MQ=255)  gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcctc                                                                                                                                        >  1:632984‑M1/1‑136 (MQ=255)  gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                                                                        >  1:911721‑M1/1‑136 (MQ=255)   ctggctcagCGCAAAGTCAATAAACTGGGCGTGATATGGTTGCATGCCTTCGCTCCTCATCTGACGTGCCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                                                                                                       <  2:1011065‑M1/140‑9 (MQ=255)    ttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                                                      <  2:201795‑M1/143‑10 (MQ=255)        tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTAT‑ACAAAGCCAAAATAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcact                                                                                                                                  >  1:610085‑M1/1‑130 (MQ=255)            cccaaATTCGATCCACCGGCGCTGTTCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                              <  1:491916‑M1/141‑18 (MQ=255)                aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTCCTGCcgtgcccaccacctcactagt                                                                                                                          >  1:534538‑M1/1‑122 (MQ=255)                  ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcctcccctgcacttcaatgatgc                                                                                                                        >  2:598414‑M1/1‑120 (MQ=255)                     aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTACAACTCTGTGTTGTAGCTCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                     >  1:199352‑M1/1‑117 (MQ=255)                     aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                     >  1:565597‑M1/1‑117 (MQ=255)                          aCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAACAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccagcacttgcatgatcggcccgtct                                                                                                                >  2:198581‑M1/1‑112 (MQ=255)                             ggCGCCGATACGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                             <  1:939179‑M1/143‑35 (MQ=255)                                gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                                                          <  2:331466‑M1/143‑38 (MQ=255)                                    atatGGTCTCATGCCTTCGCCTCTCATCTTACGTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                                                      <  2:599220‑M1/143‑42 (MQ=255)                                     tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                                                                                                     >  1:590234‑M1/1‑101 (MQ=255)                                      atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                                                                                    >  1:135694‑M1/1‑100 (MQ=255)                                        ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttaggtcctccctgc                                                                                                  >  2:627366‑M1/1‑98 (MQ=255)                                            tCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAGAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                              <  2:479393‑M1/143‑50 (MQ=255)                                              aTGCCTTCGCTCCTCCGCTTACTTTTCT‑ACAGAAAATAAAAACGCGCCCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                            <  2:221909‑M1/143‑52 (MQ=255)                                                gCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcacctgcacttcaatgatgcgccagtcttcggccatcaccgcgttcatgt                                                                                          >  2:968131‑M1/1‑90 (MQ=255)                                                       atcctcaTCTTACGTTTCTCACA‑ACAAAAAAACGCCGACTCATCAGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTCCTGCcgtcccctgcacgtcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                                                   <  2:40282‑M1/142‑61 (MQ=255)                                                         cctcATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtgatcacgacggtaatgtcggggtctg                                                                                 >  1:189703‑M1/1‑81 (MQ=255)                                                         cctcATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctg                                                                                 <  2:911721‑M1/143‑63 (MQ=255)                                                          ctcATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcctgtcggtctctgc                                                                                >  1:366153‑M1/1‑80 (MQ=255)                                                                    ttttCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcggagtct                                                                      >  1:560882‑M1/1‑70 (MQ=255)                                                                     tttCT‑ACGCACAAAAAAAAGGCGATGCAGCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgaggcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                     <  2:409683‑M1/143‑75 (MQ=255)                                                                      ttCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                                                                    <  1:298500‑M1/143‑76 (MQ=255)                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagggtcttcaacgta                                                              >  1:61023‑M1/1‑62 (MQ=255)                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcctgtcggcctccgcggcagagccttcaacccat                                                             >  2:107956‑M1/1‑61 (MQ=255)                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                           <  1:261553‑M1/143‑85 (MQ=255)                                                                                    tttaaggcGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggactctgcggcagagtcttcaacgtattccaga                                                       >  1:1026653‑M1/4‑55 (MQ=255)                                                                                    aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccaga                                                       <  1:589896‑M1/143‑89 (MQ=255)                                                                                     aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcatcgccttccagat                                                      >  2:925062‑M1/1‑54 (MQ=255)                                                                                         ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgca                                                  >  1:833917‑M1/1‑50 (MQ=255)                                                                                         ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacttcaatgatgcgcccgtcttcggtcagcactacgttcatgtcgatctctgcgacggagtcttcaacgtagtccatcaggac                                                  >  2:745495‑M1/1‑50 (MQ=255)                                                                                          gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaa                                                 >  1:984545‑M1/1‑49 (MQ=255)                                                                                                cagcggtcgccTTAAAAATCAGTCTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgct                                           <  2:364290‑M1/138‑101 (MQ=255)                                                                                                 atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtagatcactgcggaagagtcttcaacgtattccagaacgcaaaacactt                                          >  1:1029883‑M1/1‑42 (MQ=255)                                                                                                         ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                  >  1:1007226‑M1/1‑34 (MQ=255)                                                                                                         ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                  <  1:302849‑M1/143‑110 (MQ=255)                                                                                                         ccTTAAAAATCAGTCTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgtccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                  <  2:842069‑M1/143‑110 (MQ=255)                                                                                                         ccTTAAAAATCACCTCGCCAGCGCCGCCTTCTGCcctcccctgcacttccatgctgcgcccgtcttcggccctcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                  <  2:923763‑M1/143‑110 (MQ=255)                                                                                                           ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcgtcaacgtattccagatcgcacaccgcttcgccgttcac                                >  1:91366‑M1/1‑32 (MQ=255)                                                                                                                 aTCAGTTTGCCAGCGCTGCCTTCTGCcgtcacaggcaatacaaggatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                          <  1:469656‑M1/143‑118 (MQ=255)                                                                                                                      tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttccacgtattccagatcgcaaaccgcttcgccggtcacaattccgacag                     >  1:112257‑M1/1‑21 (MQ=255)                                                                                                                         gCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaatcccgccaggaa                  >  2:78991‑M1/1‑18 (MQ=255)                                                                                                                                   ccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatacccaccgtcatgtcggtctctgcggccgagtcttcaacgtactccagatcgcaaaccgcttcggcgttcaccattccgacagcacctgcggccac        >  2:652700‑M1/1‑8 (MQ=255)                                                                                                                                         gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacgattccgacagaaactgcggctaccatccc  >  1:617183‑M1/1‑2 (MQ=255)                                                                                                                                          |                                                                                                                                              TGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCT‑ACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCT  >  NC_000913/3815722‑3815999

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