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	37 [0]	[0] 37	[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)	36 (1.010)	32/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

CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAT  >  NC_000913/3815720‑3815991                                                                                                                                           |                                                                                                                                      ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAGAAAAAGGCGACTCATCAGTGGCCTTAAAAATCAGTTTGCTAGCGCCGCATTTTGCcgcc                                                                                                                                   >  1:82439‑M1/1‑139 (MQ=255)     cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                                                                                               <  2:443472‑M1/143‑9 (MQ=255)        gctAAGCTCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                                                            <  1:577164‑M1/143‑12 (MQ=255)              cgcAAATTCAATAGACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTAAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                      <  2:31547‑M1/143‑18 (MQ=255)                 aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCCCTAAAAAATCAGTTTGCCAGCGCCCCCTTCGGCcgtctcatgcccttctatga                                                                                                                   >  1:313413‑M1/1‑123 (MQ=255)                        atattctggcgCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTCCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                            <  1:517438‑M1/138‑28 (MQ=255)                        aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                            >  1:685094‑M1/1‑116 (MQ=255)                              tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgctgcgcccggcttc                                                                                                      >  1:296839‑M1/1‑110 (MQ=255)                                 cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                   <  2:180381‑M1/143‑37 (MQ=255)                                 cgcTGATATGGTGGCATGCCTTCGCTCCTCATCTTACTGTTCTACAGACAAAAGAAAGGCGACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                   <  2:337303‑M1/143‑37 (MQ=255)                                    tGATATGGTTTCATTCCTTCGCTCCTCATCTTACTTTTCTACAGACAAGAAAAAGGCGACTCATCAGGCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                <  1:279686‑M1/143‑40 (MQ=255)                                            tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCATAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcagtgatgcgcccgtcttcggtcatcactacgt                                                                                        >  2:58410‑M1/1‑96 (MQ=255)                                              tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgctcttcactgatgcgcccctcttcgttcttcactactttc                                                                                      >  1:100835‑M1/1‑94 (MQ=255)                                               cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgttca                                                                                     >  1:562646‑M1/1‑93 (MQ=255)                                                atcctttcgctcCTTATCCTCCTTTTCTACAGACAGAAAACTGGCGACTCACCAGTCGCCTTAAAATTCAGTTTGCCGGCGCCTCCTTCTGCcgtcccctgcccttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                    <  2:77427‑M1/138‑52 (MQ=255)                                                  gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccggtcttcggtcatcactacgttcatgt                                                                                  >  1:20839‑M1/1‑90 (MQ=255)                                                  gCCCTCTCTCCTCATCTTATTTGTCTACATACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                                                                  <  2:519678‑M1/143‑54 (MQ=255)                                                          tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctct                                                                          >  1:156104‑M1/1‑82 (MQ=255)                                                             tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaccttcatgtcggtctctgcg                                                                       >  2:125012‑M1/1‑79 (MQ=255)                                                                tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                    <  1:641128‑M1/143‑68 (MQ=255)                                                                      ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagaggat                                                              >  1:122631‑M1/1‑70 (MQ=255)                                                                      ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgctcatgtcggtctcagcggcagaggct                                                              >  2:558129‑M1/1‑70 (MQ=255)                                                                        ttCTCCAGACAAAAAAAAGGCGACTCGTCAGTCGCCTTAAAACGCAGTTTGCCAGCGCCGCCGTTGGCcgtcccctgcacttccatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                                                            <  1:623526‑M1/143‑76 (MQ=255)                                                                           tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTCAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagccttcaac                                                         >  1:372347‑M1/1‑65 (MQ=255)                                                                           tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                         <  1:566608‑M1/143‑79 (MQ=255)                                                                                aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttccgtcatcactacgttcatgtcgctctctgcgtccgcgacttctacgtctt                                                    >  2:675239‑M1/1‑60 (MQ=255)                                                                                  aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                  >  1:21045‑M1/1‑58 (MQ=255)                                                                                      ttaaggCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaccgtattccagcc                                              >  2:273148‑M1/3‑54 (MQ=255)                                                                                            cGACTCACCAGTCGCGTGAAAAACCAGTTTGCCAGCGCCGGCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaa                                        <  2:619203‑M1/143‑96 (MQ=255)                                                                                                         gCCTTAAAACTCTGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagttcgcaaaccgcttcgccgtt                           <  2:24602‑M1/143‑109 (MQ=255)                                                                                                           cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggccgagtcttcaacgtattccagatcgcaaaccgcttcgccgttca                         >  1:671645‑M1/1‑33 (MQ=255)                                                                                                           cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttca                         >  1:652100‑M1/1‑33 (MQ=255)                                                                                                            ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                        >  1:383298‑M1/1‑32 (MQ=255)                                                                                                                   tCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctccaacgtattccagatcgcaaaccgcttcgccggtcaccgttccg                 >  1:642535‑M1/1‑25 (MQ=255)                                                                                                                    cAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcctcggtcatcactacgttcatgtcggtgtctgcggcagagcctcccacgtctcccagctcgcacaccgcctcgcccctcaccatcccga                >  2:627995‑M1/1‑24 (MQ=255)                                                                                                                                gccgccTTCTGCcggccctctcactgcactggtgcgcccgtgtccggtcatcactacgtgcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg    <  2:535550‑M1/143‑132 (MQ=255)                                                                                                                                  cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcccaattccggcagcacctgcggct  >  1:73929‑M1/1‑10 (MQ=255)                                                                                                                                           |                                                                                                                                      CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAT  >  NC_000913/3815720‑3815991

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