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)	42 (0.950)	41/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

ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGT  >  NC_000913/3815719‑3815997                                                                                                                                            |                                                                                                                                            acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCTCCTTCTGCcgt                                                                                                                                          >  1:17111‑M1/1‑140 (MQ=255)   cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTACATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAGAGGCGACTCATCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTACTGCcgtcc                                                                                                                                        >  2:691810‑M1/1‑138 (MQ=255)    cgcttgctaAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                       <  1:448060‑M1/142‑7 (MQ=255)           tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcact                                                                                                                                <  1:601349‑M1/143‑14 (MQ=255)              gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcc                                                                                                                             >  1:4547‑M1/1‑127 (MQ=255)                    aTTCAATAAACTGGCGCTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                       <  2:126316‑M1/143‑23 (MQ=255)                       cAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                    <  2:754720‑M1/143‑26 (MQ=255)                        aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                   >  1:705748‑M1/1‑117 (MQ=255)                         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcacttatgcgcct                                                                                                                  >  1:205608‑M1/1‑116 (MQ=255)                                ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTCATACAGACAAAAAAAAGGCGACTCGTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                           <  1:539043‑M1/143‑35 (MQ=255)                                     ttatatggTTTCATGGCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                      <  2:33949‑M1/141‑40 (MQ=255)                                       atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttccgtcatca                                                                                                    >  2:79274‑M1/1‑102 (MQ=255)                                              ttCATCCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCCACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                             <  2:332005‑M1/143‑49 (MQ=255)                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                                                                        >  2:167230‑M1/1‑90 (MQ=255)                                                    ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGTCTCATCAGGCGCTTTAAAAGTCCGTGTGCCAGCACCGCCTTCGGCcggccccggcacttcaatgtggcgcccgtcttcggtcatcactacgttcatgtc                                                                                       <  1:839139‑M1/143‑55 (MQ=255)                                                     cTTCTCTCCCCATCTTGCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                                                                      <  1:391272‑M1/143‑56 (MQ=255)                                                         gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                  >  1:298830‑M1/1‑84 (MQ=255)                                                           tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctct                                                                                >  1:569903‑M1/1‑82 (MQ=255)                                                            cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactgcgttcatgtcggtctctg                                                                               >  1:160055‑M1/1‑81 (MQ=255)                                                            cctcATCTTACTTCTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctg                                                                               <  2:759421‑M1/143‑63 (MQ=255)                                                             ctcATCCGACTCTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgttgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                                                              <  1:517563‑M1/143‑64 (MQ=255)                                                              tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcg                                                                             <  2:69855‑M1/143‑65 (MQ=255)                                                                aTCTTACTTTTGTACAGACACAAACAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccggtcttcggtcatcactacgttcatgtcggtgtctgcggc                                                                           <  2:312145‑M1/143‑67 (MQ=255)                                                                            tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                               >  1:406125‑M1/1‑65 (MQ=255)                                                                            tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                               <  1:852604‑M1/143‑79 (MQ=255)                                                                               aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccggcttcggtcatcactacgttcatgtcgggctctgcggccgagtcttcaacgta                                                            >  1:654458‑M1/1‑62 (MQ=255)                                                                                gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtat                                                           >  1:284907‑M1/1‑61 (MQ=255)                                                                                    aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccggcgtcggtcatcactacgttcatgtcggtctctgcggcagagtcttccacgtattcca                                                       >  1:197136‑M1/1‑57 (MQ=255)                                                                                     aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                                                      >  1:691560‑M1/1‑56 (MQ=255)                                                                                      aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctggacttcaatgatgcgcccgtcttcggtcatcactgagttcatgtcggtctctaaggcagagcctgcaacgtcttccaga                                                     >  2:274207‑M1/1‑55 (MQ=255)                                                                                                 tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcggagtcttcaacgtattccagatcgccaaccgc                                          >  1:783970‑M1/1‑44 (MQ=255)                                                                                                    tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcgaagcgcttc                                       >  2:154881‑M1/1‑41 (MQ=255)                                                                                                     caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcg                                      >  2:155868‑M1/1‑40 (MQ=255)                                                                                                      aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgc                                     >  1:288741‑M1/1‑39 (MQ=255)                                                                                                        tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcacttggttcacccccctctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccg                                   >  1:397817‑M1/1‑37 (MQ=255)                                                                                                        tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctgcggtcatcactacgttcatgtcggtctctgcgacagagtgttccacgtattccagatagcaacccgcgtcgccg                                   >  2:460530‑M1/1‑37 (MQ=255)                                                                                                          gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgccttcaacgcattccagctcgcaaaccgctgcgccgtc                                 >  2:638307‑M1/1‑35 (MQ=255)                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcaagccgacccgttcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaatt                          >  1:629545‑M1/1‑28 (MQ=255)                                                                                                                       gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttccacgtattccatatcgcaaaccgcttcgccgctcacaattccgaca                    >  2:389533‑M1/1‑22 (MQ=255)                                                                                                                       gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgaca                    <  1:185496‑M1/143‑122 (MQ=255)                                                                                                                              aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcacaaaccgcttagccgttaacaattccgacagaaactg             >  1:522520‑M1/1‑15 (MQ=255)                                                                                                                                cgccgcCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactccg           >  1:693107‑M1/1‑13 (MQ=255)                                                                                                                                         cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatc  <  1:603589‑M1/143‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGT  >  NC_000913/3815719‑3815997

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