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	36 [0]	[0] 38	[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)	35 (0.880)	30/264	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

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815717‑3816000                                                                                                                                              |                                                                                                                                               tcaccTGCTTGCTAAGCGCAACCGCACTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACCAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                               <  1:45684‑M1/142‑2 (MQ=255)          tgctttgcgcaaatccAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                                                                                                      <  2:577924‑M1/137‑11 (MQ=255)            ctAAGCGCAAGTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCCGCcgtcccccgccc                                                                                                                                    >  2:478864‑M1/1‑131 (MQ=255)                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGAGTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCtgtctcctgtacttcaat                                                                                                                              >  1:76384‑M1/1‑125 (MQ=255)                     aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                           <  1:240865‑M1/143‑22 (MQ=255)                        gcaAGAACCTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcg                                                                                                                        <  2:145682‑M1/142‑25 (MQ=255)                                  ggCGCTGCTATGGTTTCATGTCTTCGCTCCTCATCTTACTATTCTACAGACAAAAAAAAGGCGACTGATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  2:372151‑M1/143‑35 (MQ=255)                                         atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                                                       >  1:385992‑M1/1‑102 (MQ=255)                                               tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                                                                 >  1:671244‑M1/1‑96 (MQ=255)                                                   aTGTCTTGGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGGGACTCATCAGTCGCCTTAAAAATCAGCTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                             <  2:78050‑M1/143‑52 (MQ=255)                                                   aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttagaacaccactacgttcat                                                                                             >  2:755599‑M1/1‑92 (MQ=255)                                                                   tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                             <  1:620380‑M1/143‑68 (MQ=255)                                                                       aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt                                                                         >  1:755355‑M1/1‑72 (MQ=255)                                                                        cTTTTCTACCGACAAAAGCTAAGTGCTTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                                                                        <  2:714588‑M1/143‑73 (MQ=255)                                                                         ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                                                       >  1:482799‑M1/1‑70 (MQ=255)                                                                         ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcccccgtcttcggtcatcactccgttcatgtcggtctctgcgccccagtct                                                                       >  2:337272‑M1/1‑70 (MQ=255)                                                                          tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                      >  1:141494‑M1/1‑69 (MQ=255)                                                                              tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactaatttcatgtcggtctctgcggcagagtcttcacc                                                                  >  1:26012‑M1/1‑65 (MQ=255)                                                                              tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                                  >  2:94088‑M1/1‑65 (MQ=255)                                                                              tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCACCGCCTTCTACcgtcctctgcacttcaatgatgcgcccgtcttcgggcatcactacgttcatgtaggtctctgcagcagagccttcaag                                                                  >  1:216089‑M1/1‑65 (MQ=255)                                                                                 aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaccgta                                                               >  1:137414‑M1/1‑62 (MQ=255)                                                                                     aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttccatgattcgcctgtcttcggtcatcactacgttcatgtcggtctctgcggcagcttcttcaccgcattcc                                                           >  1:398794‑M1/1‑58 (MQ=255)                                                                                     aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcgtcactacgttcatgccggtgtctgcggcagagtcttcaacgcattca                                                           >  1:635696‑M1/1‑58 (MQ=255)                                                                                         aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagat                                                       >  1:13928‑M1/1‑54 (MQ=255)                                                                                                      tcaGGCGCCTTAAAAATCAGTTTCCCAGCGCCCCCGGCTGCcgtcccctgcacttcaatgatgcgcccggcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                          <  1:711871‑M1/143‑103 (MQ=255)                                                                                                            gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgtt                                    <  1:232020‑M1/143‑109 (MQ=255)                                                                                                             ccTTAAAAATCCGCTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                   <  2:663100‑M1/143‑110 (MQ=255)                                                                                                              cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttca                                  <  1:714695‑M1/143‑111 (MQ=255)                                                                                                              cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttca                                  <  1:723019‑M1/143‑111 (MQ=255)                                                                                                               ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                                 >  2:520293‑M1/1‑32 (MQ=255)                                                                                                                 aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcacaccgcttcgccggtcacaa                               >  2:109355‑M1/1‑30 (MQ=255)                                                                                                                     aTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcc                           <  1:369643‑M1/143‑118 (MQ=255)                                                                                                                       cAGTTTGCCAGCCCCGCCTTCTGCcgtcccttgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccga                         <  2:515886‑M1/143‑120 (MQ=255)                                                                                                                        aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac                        <  2:141494‑M1/143‑121 (MQ=255)                                                                                                                         gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgaca                       <  2:776491‑M1/143‑122 (MQ=255)                                                                                                                                              ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttccagatcgcaaaccgcttcgccgttcacaattccgagccagcctgcggctaccatcccg  >  2:602215‑M1/1‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815717‑3816000

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