Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,265,001	1086 bp→GAGG		[prfA]	[prfA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1265001	1266086	1086	41 [0]	[0] 31	[prfA]	[prfA]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1265000	0 (0.000)	31 (0.510) +GAGG	28/200	1.0	100%	intergenic (+30/‑12)	hemA/prfA	glutamyl tRNA reductase/peptide chain release factor RF‑1
	?	NC_000913	1266087 =	0 (0.000)					coding (1076/1083 nt)	prfA	peptide chain release factor RF‑1

CGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGACCAACTGGCGGCGTTGTCCGAGCAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTATTCTCGCCTTT  >  NC_000913/1265972‑1266219                                                                                                                    |                                                                                                                                     acgaaatcacttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTg                                                                                                               >  2:1263417‑M2/116‑139 (MQ=255)           ttcaacaggccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGaa                                                                                                     >  1:1097295‑M2/106‑139 (MQ=255)                   gccgcccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAg                                                                                             >  1:1006946‑M2/98‑139 (MQ=255)                       cccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGc                                                                                                   <  2:392506‑M2/36‑1 (MQ=255)                       cccgtgacggggataacgaacgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGc                                                                                                   >  1:392506‑M2/94‑129 (MQ=255)                                           cgcctgaatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGccc                                                                     >  1:88755‑M2/74‑139 (MQ=255)                                                 aatattctgcgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCgcggcg                                                               <  2:88755‑M2/72‑1 (MQ=255)                                                          cgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGt                                                              <  1:616445‑M2/73‑1 (MQ=255)                                                          cgcgacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGt                                                              >  2:616445‑M2/59‑131 (MQ=255)                                                             gacagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAAt                                                   <  1:1561483‑M2/84‑1 (MQ=255)                                                               cagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAaca                                                                                                                  >  2:831427‑M2/54‑74 (MQ=255)                                                               cagcctcgggctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAaca                                                                                                                  <  1:831427‑M2/21‑1 (MQ=255)                                                                        gctggagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTCCGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGGTGAAATCCTGCTGGAGc                                        <  2:1006946‑M2/95‑1 (MQ=255)                                                                            gagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGc                                                                                            <  2:1262838‑M2/43‑1 (MQ=255)                                                                            gagtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGc                                                                                            >  1:1262838‑M2/41‑83 (MQ=255)                                                                            gagtaacagtacatcattttctttttttacagggtggaggAGGAAT‑ATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGtt                                   >  1:730848‑M2/41‑139 (MQ=255)                                                                             agtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGcc                                                                                           >  1:673002‑M2/40‑83 (MQ=255)                                                                             agtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGcc                                                                                           <  2:673002‑M2/44‑1 (MQ=255)                                                                             agtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCctgc                                              <  1:744942‑M2/89‑1 (MQ=255)                                                                             agtaacagtacatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCctgc                                              >  2:744942‑M2/40‑128 (MQ=255)                                                                                    gtacatcattttctttctttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCAGCa                            <  2:1097295‑M2/107‑1 (MQ=255)                                                                                      acatcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCaga                          <  1:1263417‑M2/109‑1 (MQ=255)                                                                                       catcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCagag                         >  1:1031029‑M2/30‑139 (MQ=255)                                                                                        atcattttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAgg                        >  2:780614‑M2/29‑139 (MQ=255)                                                                                          cattttctttttttacagggtggaggAGGAATAATAGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGc                      <  2:1031029‑M2/113‑1 (MQ=255)                                                                                           attttctttttttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCg                     <  2:1302777‑M2/114‑1 (MQ=255)                                                                                                     tttacagggtggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTAtt           <  1:780614‑M2/124‑1 (MQ=255)                                                                                                      ttacagggtggaggAGGAAT‑ATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTATtct         <  2:730848‑M2/125‑1 (MQ=255)                                                                                                              tggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTATTCTCGCCttt  >  1:220517‑M2/7‑139 (MQ=255)                                                                                                              tggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTATTCTCGCCtt   >  2:484720‑M2/7‑138 (MQ=255)                                                                                                              tggaggAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTATTCTCGCCtt   <  1:484720‑M2/132‑1 (MQ=255)                                                                                                                    |                                                                                                                                     CGCATCAACCTGACGCTCTACCGCCTGGATGAAGTGATGGAAGGTAAGCTGGATATGCTGATTGAACCGATTATCCAGGAACATCAGGCCGACCAACTGGCGGCGTTGTCCGAGCAGGAATAATGGAATATCAACACTGGTTACGTGAAGCAATAAGCCAACTTCAGGCGAGCGAAAGCCCGCGGCGTGATGCTGAAATCCTGCTGGAGCATGTTACCGGCAGAGGGCGTACTTTTATTCTCGCCTTT  >  NC_000913/1265972‑1266219

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