Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_002947	1,385,534	Δ69 bp	coding (861‑929/1266 nt)	oprD ←	basic amino acid specific porin OprD

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_002947	1385534	1385602	69	34 [0]	[0] 33	oprD	basic amino acid specific porin OprD

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 1385533	0 (0.000)	31 (1.290)	24/258	0.2	100%	coding (930/1266 nt)	oprD	basic amino acid specific porin OprD
	?	NC_002947	1385603 =	0 (0.000)					coding (860/1266 nt)	oprD	basic amino acid specific porin OprD

CCGGCGACGTTGGCGTCGGAACCGCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCA  >  NC_002947/1385384‑1385671                                                                                                                                                      |                                                                                                                                           ccGGCGACGTTGGCGTCGGAACCGCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTAGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGG                                                                                                                                            <  1:188267‑M1/150‑1 (MQ=255)             gCGTCGGAACCGCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtga                                                                                                                                <  1:256790‑M1/150‑13 (MQ=255)                tCGGAACCGCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcg                                                                                                                             >  1:42558‑M1/1‑135 (MQ=255)                               tAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaa                                                                                                                                       >  1:81100‑M1/1‑120 (MQ=255)                               tAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaa                                                                                                                                       <  2:81100‑M1/125‑6 (MQ=255)                                aaCGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtga                                                                                                                                <  2:348391‑M1/131‑13 (MQ=255)                                aaCGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtga                                                                                                                                >  1:348391‑M1/1‑119 (MQ=255)                                   gggTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccg                                                                                                          >  2:278592‑M1/1‑116 (MQ=255)                                      tCATGAAGGTCAGGCCTGGAATGCCGTATTCGACCGAGTGCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgaca                                                                                                       >  1:147609‑M1/1‑113 (MQ=255)                                                 aGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcg                                                                                            >  1:218837‑M1/1‑102 (MQ=255)                                                       ggAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttg                                                                                      <  2:618493‑M1/150‑55 (MQ=255)                                                         aaTGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACAGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtc                                                                                    <  2:264063‑M1/150‑57 (MQ=255)                                                           tGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtcca                                                                                  >  1:123100‑M1/1‑92 (MQ=255)                                                              cGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagct                                                                               >  1:469326‑M1/1‑89 (MQ=255)                                                              cGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtc                                                                                    >  2:606933‑M1/1‑89 (MQ=255)                                                              cGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtc                                                                                    <  1:606933‑M1/145‑57 (MQ=255)                                                                 aTTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgt                                                                            <  1:278592‑M1/150‑65 (MQ=255)                                                                      gCGAAGTTCAGGTAGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgcca                                                                       >  2:308599‑M1/1‑81 (MQ=255)                                                                           gTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaa                                                                  <  2:42558‑M1/150‑75 (MQ=255)                                                                           gTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgccgacgtcgaa                                                                  <  2:310653‑M1/150‑75 (MQ=255)                                                                                  tCGTAGCGAGCCCGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcgg                                                           >  2:125435‑M1/1‑69 (MQ=255)                                                                                           gCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggca                                                  <  2:149801‑M1/150‑91 (MQ=255)                                                                                           gCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggca                                                  <  1:308599‑M1/150‑91 (MQ=255)                                                                                              tGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagag                                               <  1:27468‑M1/150‑94 (MQ=255)                                                                                                      tttCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgct                                       >  2:407109‑M1/1‑49 (MQ=255)                                                                                                        tctcGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgctct                                     <  1:407109‑M1/150‑104 (MQ=255)                                                                                                         ctcGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgctctt                                    <  2:469326‑M1/150‑105 (MQ=255)                                                                                                          tcGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgctcttg                                   >  1:202528‑M1/1‑45 (MQ=255)                                                                                                          tcGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgctcttg                                   >  1:293884‑M1/1‑45 (MQ=255)                                                                                                                  ggCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgctcttggtgtcata                           <  2:293884‑M1/150‑114 (MQ=255)                                                                                                                                       cACGGAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgcggtactcggcagagccttcgctcttggtgtcatagatgttgaagtcgaacaccag      <  1:125435‑M1/150‑135 (MQ=255)                                                                                                                                           gAGTTGGCCAGGgtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcacggtactcggcagagccttcgctcttggtgtcatagatgttgaagtcgaacaccaggctt  <  2:123100‑M1/150‑139 (MQ=255)                                                                                                                                                      ggtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgc                                                             >  1:610722‑M1/1‑1 (MQ=255)                                                                                                                                                      ggtgaaggtgtgagcgccgatgttgtaggcgcccgacaggctgaacgcgcggttgtccagcttgtcgccatcgaaggcgc                                                             <  2:610722‑M1/80‑80 (MQ=255)                                                                                                                                                      |                                                                                                                                           CCGGCGACGTTGGCGTCGGAACCGCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGAAGATGGTGCCGCCACCGTCGATACCATAACCATAGTCGCCGTCGCCGGTGACTTTCTGGTAGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCA  >  NC_002947/1385384‑1385671

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