Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_002947	3,788,572	Δ23 bp	100%	coding (41‑63/1392 nt)	PP_3350 ←	hypothetical protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_002947	3788572	3788594	23	53 [0]	[0] 53	PP_3350	hypothetical protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3788571	0 (0.000)	49 (0.950)	40/264	NT	100%	coding (64/1392 nt)	PP_3350	hypothetical protein
	?	NC_002947	3788595 =	0 (0.000)					coding (40/1392 nt)	PP_3350	hypothetical protein

GGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTG  >  NC_002947/3788448‑3788739                                                                                                                                                    |                                                                                                                                                 ccaggccgactctgacagccgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtacc                                                                                                                                                 <  6:288008‑M2/2‑1 (MQ=255)   aggccgactctgacagctgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCAt                                                                                                                                               <  4:99183‑M2/4‑1 (MQ=255)     gccgactctgacagccgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGa                                                                                                                                             <  3:35732‑M2/6‑1 (MQ=255)      ccgactctgacagccgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAg                                                                                                                                            <  5:72149‑M2/7‑1 (MQ=255)           tctgacagccgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAgggg                                                                                                                                         >  7:199924‑M2/138‑147 (MQ=255)               acatcagattgtagctcaggccatcgaaaggggcagtgtccttggcggtgtcgtcgagaaagcggtagtcctcgatccagcgagtggcccctggcctccccggccttctctcttcaagcgggctggccgcgtaCCATGAGGggcag                                                                                                                                      <  1:205460‑M2/13‑1 (MQ=255)               acagccgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGggcag                                                                                                                                      <  8:199924‑M2/13‑1 (MQ=255)                cagccgatggtagcgcaggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagacaacgatccagcgaggggcccctggcctccccggccggctcactacaagcgggctggccgcgtaCCATGAGGggaaggca                                                                                                                                   >  4:16977‑M2/133‑148 (MQ=255)                                 ggccatcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCtt                                                                                                                 <  2:234213‑M2/34‑1 (MQ=255)                                    catcgaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGat                                                                                                              >  1:8307‑M2/113‑149 (MQ=255)                                        gaaaggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGatcatc                                                                                                          <  7:122034‑M2/41‑1 (MQ=255)                                            ggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGaa                                                                                                      >  6:131668‑M2/105‑149 (MQ=255)                                            ggggcagtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGaa                                                                                                      >  3:184777‑M2/105‑149 (MQ=255)                                                 agtggccttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTcc                                                                                                 <  2:79637‑M2/50‑1 (MQ=255)                                                        ttggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATAAGAATTTCACTacaat                                                                                          >  8:200591‑M2/93‑147 (MQ=255)                                                          ggcggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTg                                                                                        >  3:119505‑M2/91‑149 (MQ=255)                                                            cggggtcgtcgagaaagcggtagtcctcgatccagcgagggtcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGtt                                                                                      <  4:77971‑M2/61‑1 (MQ=255)                                                            cggggtcgtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGtt                                                                                      <  4:184777‑M2/61‑1 (MQ=255)                                                                   gtcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATaa                                                                               >  8:121932‑M2/82‑149 (MQ=255)                                                                    tcgagaaagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAAt                                                                              >  5:325013‑M2/81‑149 (MQ=255)                                                                          aagcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCt                                                                        <  7:121932‑M2/75‑1 (MQ=255)                                                                           agcggtagtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCt                                                                        >  4:226967‑M2/74‑148 (MQ=255)                                                                                  gtcctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACt                                                                <  3:309345‑M2/83‑1 (MQ=255)                                                                                     ctcgctccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTtat                                                                                 <  7:142880‑M2/66‑1 (MQ=255)                                                                                     ctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTtat                                                                                 >  8:142880‑M2/64‑129 (MQ=255)                                                                                     ctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTgg                                                              <  1:41841‑M2/85‑1 (MQ=255)                                                                                     ctcgatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGa                                                             >  1:162864‑M2/64‑149 (MQ=255)                                                                                        gatccagcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGAt                                                            <  2:8307‑M2/87‑1 (MQ=255)                                                                                              gcgaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGt                                                    <  2:56039‑M2/95‑1 (MQ=255)                                                                                                gaggggcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGtt                                                   <  3:16977‑M2/96‑1 (MQ=255)                                                                                                     gcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTtat                                                                                 >  8:207302‑M2/48‑113 (MQ=255)                                                                                                     gcccctggcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGat                                             >  5:171336‑M2/48‑149 (MQ=255)                                                                                                            gcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTAttttt                                      <  5:131668‑M2/109‑1 (MQ=255)                                                                                                            gcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTAttttt                                      <  4:57189‑M2/109‑1 (MQ=255)                                                                                                            gcctccccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAAGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTAttttt                                      <  3:164062‑M2/109‑1 (MQ=255)                                                                                                                  ccggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAAc                                <  2:133533‑M2/115‑1 (MQ=255)                                                                                                                   cggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTAc                                                                     >  3:209352‑M2/34‑111 (MQ=255)                                                                                                                   cggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTAc                                                                     <  4:209352‑M2/78‑1 (MQ=255)                                                                                                                   cggccggctctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAAcc                               <  6:325013‑M2/116‑1 (MQ=255)                                                                                                                           tctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATgcgc                       >  5:288688‑M2/26‑149 (MQ=255)                                                                                                                           tctcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATgcgc                       >  2:80272‑M2/26‑149 (MQ=255)                                                                                                                             tcttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCg                     >  3:65122‑M2/24‑149 (MQ=255)                                                                                                                               ttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGt                    >  6:4317‑M2/22‑148 (MQ=255)                                                                                                                               ttcaagcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGt                    >  7:45346‑M2/22‑148 (MQ=255)                                                                                                                                   agcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCaa                 >  5:100086‑M2/18‑147 (MQ=255)                                                                                                                                   agcgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATa               <  5:4317‑M2/132‑1 (MQ=255)                                                                                                                                     cgggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTAGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGc             <  1:80272‑M2/134‑1 (MQ=255)                                                                                                                                      gggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCg            <  7:200591‑M2/135‑1 (MQ=255)                                                                                                                                      gggctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCg            >  8:151294‑M2/15‑149 (MQ=255)                                                                                                                                        gctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTc                                                                                             <  3:245123‑M2/54‑1 (MQ=255)                                                                                                                                        gctggccgcgtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTc                                                                                             >  4:245123‑M2/13‑66 (MQ=255)                                                                                                                                                 gtaCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCgtg  >  8:278578‑M2/4‑148 (MQ=255)                                                                                                                                                   aCCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGAt                                                            <  6:153176‑M2/87‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTG  >  NC_002947/3788448‑3788739

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