Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	3,815,859	Δ82 bp	100%		[rph]–[rph]	[rph], [rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3815859	3815940	82	77 [0]	[0] 77	[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)	74 (1.550)	53/394	NT	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

AGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTC  >  NC_000913/3815632‑3816085                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    aGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                                                                                                                    <  2:366701‑M1/228‑2 (MQ=255)       cgcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCCCCTAATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTTCCCGCGCCGCCTTCTGCcgccccc                                                                                                                                                                                                                              >  2:468640‑M1/1‑221 (MQ=255)       cgcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAAGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                                                                                                                                                              >  1:178332‑M1/1‑221 (MQ=255)          ccGGTATTAAACAGCCCGGCTTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                                                                                                                                            <  1:614890‑M1/227‑10 (MQ=255)             gTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCCTCcgtcccccgcac                                                                                                                                                                                                                         >  2:615616‑M1/1‑215 (MQ=255)                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTTAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                                                                                                                                                                                           >  2:118333‑M1/1‑204 (MQ=255)                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCCACTCATCAGTCGCCTTAAAAATCATTTTGACAGCGCCGCCTTCTTCcgtcccctgcacttcaatgaatga                                                                                                                                                                                                             >  2:230302‑M1/1‑204 (MQ=255)                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCCGCcgccccctgcacttcaatgatgc                                                                                                                                                                                                              >  2:168766‑M1/1‑204 (MQ=255)                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCACCGCCGCCTTCTCCcgtaccctgcacttcaataatggg                                                                                                                                                                                                             >  2:638711‑M1/1‑204 (MQ=255)                        cccGGCGTTGAAGAAATAGGGGCTTTTGAGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCACTCATCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATTAGTTCGCCATCGCCCACTTCTGCggtcccctgcacttcaatcatacg                                                                                                                                                                                                             >  1:53695‑M1/1‑204 (MQ=255)                         ccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTTAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                                                                                                                                                          <  1:118333‑M1/214‑12 (MQ=255)                         ccGGCGTTAAAAAAAAAGGGGCGTTTGCGCACGGATCTCCGCGGAAACTCGCCAAGCTTTAACACCTGCTTGCTAAGCGCAAATTCAAGAAACTGGCGCTGATAGGTTTTCATGCCTTCGCTCCTCATCGTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                                                                                                            <  2:611350‑M1/228‑26 (MQ=255)                          cGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                                                                                                            <  1:468640‑M1/227‑26 (MQ=255)                          cGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCACCCTTCTGCcgccccctgcacgtcaatgatgcgcc                                                                                                                                                                                                           >  1:651876‑M1/1‑202 (MQ=255)                                 gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                                                                                                                                                >  1:633126‑M1/1‑195 (MQ=255)                                  aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                                                                                                               <  2:633126‑M1/200‑7 (MQ=255)                                    gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                                                                                                                   <  1:63506‑M1/226‑35 (MQ=255)                                      aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                                                                                                                                                                                               <  1:586070‑M1/228‑39 (MQ=255)                                       aTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTCCCATCGCCGCCTTCTCCcgtcccctgcacttcaatgatgcgcccctcttcggtcat                                                                                                                                                                                              >  1:590319‑M1/1‑189 (MQ=255)                                            ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatccacccgtcttcgatcatcacta                                                                                                                                                                                         >  1:170348‑M1/1‑184 (MQ=255)                                            ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                                                                                                                                                         >  1:324248‑M1/1‑184 (MQ=255)                                              cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCCAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCATCGCCCCCTTCTGCcgtcccctgcacttcaatgatgcgcccgccttcggccatcactac                                                                                                                                                                                        >  2:78816‑M1/1‑182 (MQ=255)                                                tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAACTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTTATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgt                                                                                                                                                                                      >  1:460706‑M1/1‑180 (MQ=255)                                                tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATTAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCACCGCCGCCTTCTGCcgtcccctgcacttcactgatgcgcccgccttcggccatctctccgt                                                                                                                                                                                      >  1:52415‑M1/1‑180 (MQ=255)                                                     gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCATCTGCcctaccctgcacttcaattatgcgcccgtcttcggtcacacctacgttcatgc                                                                                                                                                                                >  2:579963‑M1/1‑175 (MQ=255)                                                          gATTTCAGTGTAAACCTGCCAACATTAAACCCCGTCTTGCTAAGCGCAAAATCAATGACCTGGCGCTGATAGGGGTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAGGGCGACGCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt                                                                                                                                                                            <  2:33808‑M1/227‑58 (MQ=255)                                                           aTTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                                                                                                          <  1:169850‑M1/228‑60 (MQ=255)                                                           aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                                                                                                          <  1:544780‑M1/228‑60 (MQ=255)                                                           aTTTCAGCGTAAACTCGCCAAACTTGAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                                                                                                                                          <  1:617502‑M1/228‑60 (MQ=255)                                                                   gTAAACTCGCCAAACTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                                                                                                                                                  <  1:242901‑M1/228‑68 (MQ=255)                                                                   gTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactccgttcatttcggttcctgcggc                                                                                                                                                                   >  1:549475‑M1/1‑161 (MQ=255)                                                                     taacTCCCCAAACTTTAACCCCTGCTTGCTAAGCGCAAATTAAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                                                                                                                <  2:651876‑M1/227‑70 (MQ=255)                                                                     aaaCTCTCCAACATTTAACCCCTGTTTGCTAAGCGCAAATTCAATAAACTGGCCCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                                                                                                                <  2:558627‑M1/228‑70 (MQ=255)                                                                      aaCCCGCCAAACTTTAGCGCCTTGCTGCGAAGGGCAAATTCAAAAGACTGGCGCTGTTGTTGTTTCCTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                                                                                                                <  2:52415‑M1/227‑70 (MQ=255)                                                                          ctccaacacttaacccctgcTTGCTAAGCGGAAATTCAAAAAACTGGCGATGATATGGTTTCATGCCTTCGCCCCTCATCTTACTTTTCTACAGGCAAAAAAAAGGCGACTCATCAGTCGACTTAAAAATCAGTTGGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                                                                                                           <  2:368368‑M1/219‑75 (MQ=255)                                                                            ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgcaagagtcttca                                                                                                                                                         >  1:412619‑M1/1‑152 (MQ=255)                                                                                            gcttgctAAGCGCACATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgttcatgtcggtctctgcgtcagagtcttcaacggattccaagtcc                                                                                                                                          >  1:244086‑M1/1‑136 (MQ=255)                                                                                            gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcc                                                                                                                                          >  1:24441‑M1/1‑136 (MQ=255)                                                                                            gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgccggtccctgcggcagagtcttcacccgattccaaaccc                                                                                                                                          >  2:262301‑M1/1‑136 (MQ=255)                                                                                                    aGCGCAAATTCAATAAACTGGCGCTGATATGGATTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAACAAAGGAGACTCATCAGTCGCCTTCAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccccttcacttcaatgatgcgcccgtctccgttcatcactacgatcctgtcattctctgcggccgagtctttcacgtattccaaaacgcaaaccgct                                                                                                                                 >  2:189950‑M1/1‑128 (MQ=255)                                                                                                      cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                                                                                                               >  1:551571‑M1/1‑126 (MQ=255)                                                                                                              taatAAACTGGCGTTTATATGGTTTCATGCCTTCGCTCCTCATCTTCATTTTGATCAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGGGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgccttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttca                                                                                                                       <  2:590319‑M1/227‑111 (MQ=255)                                                                                                                      tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttccgaattccg                                                                                                               >  1:431449‑M1/1‑110 (MQ=255)                                                                                                                       ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaccgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccg                                                                                                              >  1:393113‑M1/1‑109 (MQ=255)                                                                                                                       ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                                                                                                               >  2:623559‑M1/1‑108 (MQ=255)                                                                                                                        gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                                                                                                                                                              <  1:623559‑M1/114‑8 (MQ=255)                                                                                                                                       cacttcttcgaCCCACATCTTGCATTTCCGCAAGCAAAAAAAAGGCGCCTCACTGGACTCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCcgtcccctgcccgtcaatgatgcgcccgtctgcggtcatcactaggttcatgtcgttctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggcta                                                                                               <  2:53695‑M1/222‑135 (MQ=255)                                                                                                                                       cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCCACTCATCAGTCGCTGTAAAAATCAGTTTGCCAGCGCCTCCTTCTGCcgtccccttcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtctgtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc                                                                                                                        <  1:351002‑M1/202‑110 (MQ=255)                                                                                                                                        aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctcctccgttcacaattacgccagaaactgcgcgtga                                                                                              >  1:502982‑M1/1‑92 (MQ=255)                                                                                                                                        aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaagttcctgtcggtctctgcggcagagtattcaccggattccagatcccaaaacgcttccccgttcaccattccgacagaaactgcggctacc                                                                                             >  2:425922‑M1/1‑92 (MQ=255)                                                                                                                                         tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca                                                                                            >  1:555396‑M1/1‑91 (MQ=255)                                                                                                                                         tGCCTTCGCTCCTCATCCTACGTTTCCACCGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca                                                                                            <  1:246398‑M1/228‑138 (MQ=255)                                                                                                                                          tccctccctcctcatttttcTTTTACACAGACAAAAAAAAGGCTGCTCATCAGGCGCCTTAAAACTCAGTTTGCCAGCGCCGACTTCTGCcttcccctgcacttcaatgatgcgcccgtcttcggtcatcacgacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca                                                                                            <  2:502982‑M1/220‑138 (MQ=255)                                                                                                                                             ttCGCTCCTCATCTTACTTTTTCTCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAGCAGTTTGCCCGCGCCGCCTTCTGCcggcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccc                                                                                        <  2:324248‑M1/228‑142 (MQ=255)                                                                                                                                             ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgcattccagatcgcaaaccgcttcgccgttaccaattccgacaaaaaatgcggccaccacccc                                                                                        >  2:395929‑M1/1‑87 (MQ=255)                                                                                                                                                       tgtgtacttTTC‑ACAGACAAAAAAAAGGCTACTCATCAGGCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctggacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtctgtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggat                                                                             <  1:160966‑M1/224‑153 (MQ=255)                                                                                                                                                        tCTTATTTTTTCACAGACAAAAAAAAGGCGACCCAACAGTCGTCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCcgtcccctggacttcaatgatgcgcccgtcttcggtcatcacttcgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggat                                                                             <  2:587801‑M1/228‑153 (MQ=255)                                                                                                                                                                 tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcccaaaccgcttcgccgttcacaattccgacagaaactgcggctacactccctttaatcggattggttttc                                                                     >  1:230171‑M1/1‑67 (MQ=255)                                                                                                                                                                  cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccct                                                                                       >  2:426496‑M1/1‑66 (MQ=255)                                                                                                                                                                   taaagacAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctt                                                                                      <  1:426496‑M1/205‑144 (MQ=255)                                                                                                                                                                    aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgccagaaactgaggctaccatccctttcatcggattggttttcacct                                                                 >  1:11352‑M1/1‑64 (MQ=255)                                                                                                                                                                     cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagccttcaacgtattccagatcgcaaaccgcttccccgttcccaattcccaccgaaaccgcggccaccatacccttcaccggattgggtttcgcct                                                                 >  2:238486‑M1/1‑63 (MQ=255)                                                                                                                                                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcccaattccgacagaacctgcggctaccatccctttcaccggcttggttttcaccttgc                                                              >  1:241304‑M1/1‑60 (MQ=255)                                                                                                                                                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgcc                                                             >  2:520960‑M1/1‑60 (MQ=255)                                                                                                                                                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgc                                                              >  1:5574‑M1/1‑60 (MQ=255)                                                                                                                                                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatctccgtctctgcggcagagtcttcaccgtattccagatcgcaaaccgcttccccgttccccattcctcccgaaactgcgcctcccatccctttcatcggcttggttttcaaattcc                                                              >  2:57151‑M1/1‑60 (MQ=255)                                                                                                                                                                                      aCTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcacagtcttcaacgtattccagatcgcaaaccgcttccccgttaccagttccgcccgacactgcggccaccatccctttcaccggatagggtttagacttgccgttttccccaccc                                                >  2:112475‑M1/1‑46 (MQ=255)                                                                                                                                                                                           tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgt                                          <  1:520960‑M1/228‑188 (MQ=255)                                                                                                                                                                                            tcgtcTCCTTAAAAATTAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacgtcaatgatgcgcccttcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgta                                         <  1:189950‑M1/226‑189 (MQ=255)                                                                                                                                                                                                  ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgtagcgca                                    >  1:389298‑M1/1‑34 (MQ=255)                                                                                                                                                                                                  ccTTAAAAAGCAGGTTGCCCGCGCCCCCTTTGGCcgtcccctgcacttcaatgatgcgcccggcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggtttttagcttgccgttttccaccagcttctgtagcgca                                    <  1:78816‑M1/227‑194 (MQ=255)                                                                                                                                                                                                    ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatctcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggaatggttttcagcttgccgttttccacccgcttctgaagcgccaca                                 >  2:231483‑M1/1‑32 (MQ=255)                                                                                                                                                                                                         aaTCAGTTTGCCAGCGCCCGCTTCGGCcgtcccctgcccttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgtagcgcatctaccag                            <  1:504770‑M1/228‑202 (MQ=255)                                                                                                                                                                                                           tCATTTTGCCAGCGCCCCCTTCGTCcgtcccctgctcttcaatgagtcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgtagcgcatctaccagcg                          <  2:555396‑M1/228‑204 (MQ=255)                                                                                                                                                                                                             aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgccagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgtagcgcatctcccagcgc                         >  1:535872‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                 tGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactccggctaccatccctttcatcggattgggtttaaccttgccgttttccaccatcttctgtagcccatctaccacagcc                        >  1:327511‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                   ccgtcccctgcccttcaatgttgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcccaattccgacagaaactgcggctaccatccctttcatcggattggttttcagcttgccgttttccaccagcttctgtagcgcatctaccagcgccacgcaggcacccgtaatcgacg  <  1:262553‑M1/228‑228 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    AGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTC  >  NC_000913/3815632‑3816085

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑