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	42 [0]	[0] 42	[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)	41 (1.190)	35/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

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGG  >  NC_000913/3815717‑3815996                                                                                                                                              |                                                                                                                                           ccaccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCACGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                           <  1:530746‑M1/142‑2 (MQ=255) acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                           <  1:141101‑M1/143‑2 (MQ=255)    ccTGCTTGCTAAGCGCAAATTCAAGAAACTGGCGCTGATATGGTTTCACGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                                                                        <  1:544325‑M1/143‑5 (MQ=255)      tgcttgctAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCTGCcgtccc                                                                                                                                      >  1:590363‑M1/1‑137 (MQ=255)      tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCACTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                      >  1:448629‑M1/1‑137 (MQ=255)                  gacaattcAACAGACTGGCGCTGATATGGTGTCATGCCTTCGCTGCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                          <  2:170078‑M1/140‑19 (MQ=255)                    aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCCGCcgtcacctgcactccaatga                                                                                                                        >  2:501882‑M1/1‑123 (MQ=255)                        tCAATAAAGTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTGCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcg                                                                                                                    <  2:183305‑M1/143‑25 (MQ=255)                             caacTCGCCCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACACACAAATAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                                                                                               <  2:340727‑M1/142‑30 (MQ=255)                             aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCACCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                                                                                               <  2:373699‑M1/143‑30 (MQ=255)                                 tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGCCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctccacttcaatgatgcgcccgtcttc                                                                                                           >  1:162584‑M1/1‑110 (MQ=255)                                        gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCTGCcgtcccccgcccctccatgctgcacccgtcgtcggcactg                                                                                                    >  1:512126‑M1/1‑103 (MQ=255)                                        gATATGATTCCACGTCTTCGCTCCTCACCTTACTTTTCTACAGACACAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCACTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                                                                    <  1:245195‑M1/143‑41 (MQ=255)                                           atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgacgcgcccgcccccgctcagcact                                                                                                 >  2:540741‑M1/1‑100 (MQ=255)                                            tGGTTTCATGCCTTCGCTCCTCATCTTACGCGTCTACAGACAAAAAAGACGCGACTCATCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                                                                <  2:264628‑M1/143‑45 (MQ=255)                                               tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                                                             >  1:630228‑M1/1‑96 (MQ=255)                                                        ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatatcgg                                                                                    >  2:827‑M1/1‑87 (MQ=255)                                                                   tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcagtgatgcccccgcctccggccctcactacgttcatgtggcgccctgcggca                                                                         >  2:580961‑M1/1‑76 (MQ=255)                                                                     ttACTGTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                       <  2:571818‑M1/143‑70 (MQ=255)                                                                             cTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgttccctgcacttcaatgatgcgcccgtcgtcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaa                                                               <  1:488570‑M1/143‑78 (MQ=255)                                                                              ttaagacAGAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                              <  2:67181‑M1/140‑79 (MQ=255)                                                                              tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                              <  2:179268‑M1/143‑79 (MQ=255)                                                                               aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactccgtccgtgtcggtctctgcggcagactcttcaacg                                                             >  1:294045‑M1/1‑64 (MQ=255)                                                                                 aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttccctgatgcgcccgtcttcggtcttccctacgttcatgtcggtcactacagcagaatattcaccgta                                                           >  1:520714‑M1/1‑62 (MQ=255)                                                                                 aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgta                                                           <  1:368458‑M1/143‑82 (MQ=255)                                                                                   aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgaacttcaatgatgcgcccgtcttcggtcatccctacgttcatgtcggcctacccgcctgaggcttcaacgactt                                                         >  2:261947‑M1/1‑60 (MQ=255)                                                                                     aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGATTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcgtcactacgttcatgtcgctctctgcggcagagacatcaacgtattcc                                                       >  2:321331‑M1/1‑58 (MQ=255)                                                                                        aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcatcccctgcacttcaatgatgcgcccgtcttcggtcatcactaccatcatgacgttgtctgcggcagagtcttcaacgtattccaga                                                    <  1:641466‑M1/143‑89 (MQ=255)                                                                                           aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCCGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcg                                                 >  1:28522‑M1/1‑52 (MQ=255)                                                                                                  cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtagtccagatcgcaaaccg                                          >  2:124630‑M1/1‑45 (MQ=255)                                                                                                     atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                       >  1:361203‑M1/1‑42 (MQ=255)                                                                                                     atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatagcaaaccgctg                                       >  1:361209‑M1/1‑42 (MQ=255)                                                                                                     atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccatgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                       >  2:348953‑M1/1‑42 (MQ=255)                                                                                                      tcaGTCGTGTTAAAAATCAGTGTGCCAGCGGCGCCTTCTGCcgtcccctgcacttcagtgatgcgcccgtcttcggtcatcattacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                      <  2:176655‑M1/143‑103 (MQ=255)                                                                                                          tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcctcaacgtcttccagatcgcagaccgcttcgccc                                  >  1:44258‑M1/1‑37 (MQ=255)                                                                                                              cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttccacgtcttccagctcgcacaccgcttcgcctttca                              >  1:542207‑M1/1‑33 (MQ=255)                                                                                                               ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                             >  1:477569‑M1/1‑32 (MQ=255)                                                                                                                     aTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                       >  1:481144‑M1/1‑26 (MQ=255)                                                                                                                          tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacacttcctccag                  >  1:578193‑M1/1‑21 (MQ=255)                                                                                                                           ttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcctgtcggtctctgcggcagcgtcttcaacgtattcccgatcgcgaaccgcttcgccgttcccacttccgccaca                 >  1:165433‑M1/1‑20 (MQ=255)                                                                                                                             gCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaga               >  1:269685‑M1/1‑18 (MQ=255)                                                                                                                                          tCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcaccattccgacagaaactgcggctaccat  >  1:400691‑M1/1‑5 (MQ=255)                                                                                                                                              |                                                                                                                                           ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGG  >  NC_000913/3815717‑3815996

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