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	33 [0]	[0] 33	[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)	33 (0.880)	32/264	0.2	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

GCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA  >  NC_000913/3815732‑3815990                                                                                                                               |                                                                                                                                     gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGCCGCCTCAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                                                                                                      >  1:724207‑M1/1‑127 (MQ=255)       tttcAATAACCCGGCGCTGATATGGTTTCATGCCTTCGCTCCGCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                <  1:643036‑M1/142‑23 (MQ=255)        ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCATGATCTTACTTTTGTAGAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgc                                                                                                               >  1:124354‑M1/1‑120 (MQ=255)         tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatcatgcg                                                                                                              >  1:493297‑M1/1‑119 (MQ=255)            aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                           <  2:264459‑M1/143‑28 (MQ=255)              aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctggacttcagtggtgcccccgt                                                                                                         >  1:391160‑M1/1‑114 (MQ=255)                                  tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                                                     <  1:590644‑M1/143‑50 (MQ=255)                                      gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgaagttctacttttggcccgtcttcggtcatcactacgttcatgt                                                                                 >  1:495704‑M1/1‑90 (MQ=255)                                         ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCTACTCATCAGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTACTACcgtccccggcacttcaatgatgcgcccgccttcggtcatcactacgttcatgtcgg                                                                              >  2:484870‑M1/1‑87 (MQ=255)                                                 tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAATTTGCCACCGCCGCCTTCTGCcgtcccctgcacttccatgatgcgcccgtcttcggccatccctacgttcacgccggtgcctgcg                                                                      >  2:234165‑M1/1‑79 (MQ=255)                                                  cATCTTACTGTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgg                                                                     <  2:395361‑M1/143‑66 (MQ=255)                                                        aCTTTTCTACAGACAAAAAAAAGGCGACGCATCAGTCGCCTTAAAACTCATTTTGCCTGCGGCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt                                                               <  2:352320‑M1/143‑72 (MQ=255)                                                         cttatgaacagagaaAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttctatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                                                              <  1:232792‑M1/136‑73 (MQ=255)                                                          ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                                             >  1:396083‑M1/1‑70 (MQ=255)                                                            ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatggcggtctctgcggcagagtcttc                                                           >  1:220643‑M1/1‑68 (MQ=255)                                                                aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                                                       >  1:417489‑M1/1‑64 (MQ=255)                                                                  aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTGGCCAGCACCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgta                                                     >  1:543122‑M1/1‑62 (MQ=255)                                                                   gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaccgtat                                                    >  1:226382‑M1/1‑61 (MQ=255)                                                                      cagaaaaaggcGACTCATCAGTCGCCTTAAAAGTCAGTTTGCCAGCGCCGCCTTCTGCcttgccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                 <  2:634293‑M1/140‑86 (MQ=255)                                                                         aaaaaTGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcccga                                              >  2:560173‑M1/1‑55 (MQ=255)                                                                          aaaaGGCTACTCATCAGTAGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                             <  1:492388‑M1/143‑90 (MQ=255)                                                                           agaggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccttcacttctatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                                            <  2:634682‑M1/141‑91 (MQ=255)                                                                           aaaGGCGCCTCCTCAGTCGCCTTAAAGATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatggttcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                                            <  2:101056‑M1/143‑91 (MQ=255)                                                                                    tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtatgtgcggctgagtcttcaacgtattccagatcgccaaccgt                                   >  1:495683‑M1/1‑44 (MQ=255)                                                                                      atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccaatcaacagacatcgctacgcccacctcgttctttgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                 >  1:200573‑M1/1‑42 (MQ=255)                                                                                       tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCAGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                <  1:54499‑M1/143‑103 (MQ=255)                                                                                           tCGCCTTAAACATCCGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccg                            <  2:561969‑M1/143‑107 (MQ=255)                                                                                                   aagaccagtttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaat                    <  1:437739‑M1/138‑115 (MQ=255)                                                                                                     aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtcttcccgatcgcaaaccgctccgcccttcccaattc                  >  2:166095‑M1/1‑27 (MQ=255)                                                                                                       ccaGTTCGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaacgatgcgcccgtcttcggtcatcactacgtccatgtcggtctctgcggcacagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccg                <  2:462052‑M1/142‑119 (MQ=255)                                                                                                         aGTTTGCCTGCGCCGCCTTCTGCcgtcccctgcacgtcaatggtgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac              <  2:727005‑M1/143‑121 (MQ=255)                                                                                                                 aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgccagaaactg      >  1:506522‑M1/1‑15 (MQ=255)                                                                                                                     ccgccTTCTGCcctcccctgcacttcattgatgcgcccgtcttcggtcatcactacgttcatgtcggtcgctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggc  <  2:194192‑M1/143‑133 (MQ=255)                                                                                                                               |                                                                                                                                     GCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA  >  NC_000913/3815732‑3815990

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