Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_vch_tpiA	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_vch_tpiA	1978503–1979270	1979278	9–776	43 [0]	[1] 44	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_tpiA	= 1978502	0 (0.000)	43 (1.190)	38/268	0.1	100%	noncoding (1978502/4641682 nt)	repeat_region	repeat region
	?	NC_000913_3_vch_tpiA	1979279 =	0 (0.000)					noncoding (1979279/4641682 nt)	repeat_region	repeat region

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_vch_tpiA/1978363‑1978641                                                                                                                                            |                                                                                                                                            gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGGCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                    >  2:343491‑M1/1‑140 (MQ=255)  tCTTTCCACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                   >  2:162922‑M1/1‑139 (MQ=255)     ttACACTTATCTAAGATTTTTTTTATTTTTTCCCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                >  2:45089‑M1/1‑136 (MQ=255)             aTCTAAGATTGTTCGTAAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTAATTAACATAAATGtataagtcatacttttgtttt                                                                                                                        <  2:317571‑M1/149‑22 (MQ=255)              tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                       >  1:450166‑M1/1‑127 (MQ=255)              tCTAAGATTTTTCCTAAATCGACGCAACTGTACCCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtaataattttgttttt                                                                                                                       >  2:153658‑M1/1‑127 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                     >  1:530185‑M1/1‑125 (MQ=255)                         tcataaaTAGACGAAACAGTACTCGTCACTAGACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                            <  2:303184‑M1/146‑34 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATATATGtataagtcttacttttgttttgggtgtctttcaat                                                                                                          >  2:464882‑M1/1‑114 (MQ=255)                              aaTCGACGCAACTGTACTCGTCACTTCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTGACATAAACGtataagtcatacttttgttatgggtgtatttcaatctg                                                                                                       >  1:513443‑M1/1‑111 (MQ=255)                                   aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcatacttttgttttcggtgtatttcaatctgtaaaa                                                                                                  >  1:74852‑M1/1‑106 (MQ=255)                                     gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                <  1:547494‑M1/149‑46 (MQ=255)                                         cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                            <  1:573876‑M1/149‑50 (MQ=255)                                          tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatgtgttaaaaagtttt                                                                                           >  1:6904‑M1/1‑99 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGCGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                      <  1:236095‑M1/149‑56 (MQ=255)                                                cGTCACTACACGCACATACAACGGAGGGGTTCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                     <  1:595515‑M1/149‑57 (MQ=255)                                                     cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                <  1:349645‑M1/149‑62 (MQ=255)                                                        cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                                                             >  1:439762‑M1/1‑85 (MQ=255)                                                                tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                     <  1:464882‑M1/149‑73 (MQ=255)                                                                tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                     >  1:365788‑M1/1‑77 (MQ=255)                                                                 aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttgcatcgtgtgaaaggtttttcgctacgctagcaagcaaaaat                                                                    <  1:162922‑M1/149‑74 (MQ=255)                                                                   aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgccagcaagcaaaaatga                                                                  >  2:520008‑M1/1‑74 (MQ=255)                                                                     cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                <  2:530185‑M1/149‑78 (MQ=255)                                                                      ggAGGGGGGCTCCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagctccgctagccagcaaaaatgaaat                                                               >  1:415795‑M1/1‑71 (MQ=255)                                                                       gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                              <  1:343227‑M1/149‑80 (MQ=255)                                                                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             <  1:77037‑M1/149‑81 (MQ=255)                                                                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             <  1:581167‑M1/149‑81 (MQ=255)                                                                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             <  1:45089‑M1/149‑81 (MQ=255)                                                                            gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                         >  2:321064‑M1/1‑65 (MQ=255)                                                                             ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                        <  1:343491‑M1/149‑86 (MQ=255)                                                                                   gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                  >  1:336989‑M1/1‑58 (MQ=255)                                                                                        tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                      <  2:309602‑M1/108‑56 (MQ=255)                                                                                        tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                      >  1:309602‑M1/1‑53 (MQ=255)                                                                                          aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                              <  1:589139‑M1/98‑48 (MQ=255)                                                                                          aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                              >  2:589139‑M1/1‑51 (MQ=255)                                                                                             aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                          >  2:630126‑M1/1‑48 (MQ=255)                                                                                             aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                          <  1:630126‑M1/83‑36 (MQ=255)                                                                                                    gATGCAGATCACACAAAACACGCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 <  2:660613‑M1/149‑109 (MQ=255)                                                                                                     aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttttttactttatcaagcaaaaatgaatcaggatgaatcgtaatgggatgttgcggatgg                                >  2:625170‑M1/1‑40 (MQ=255)                                                                                                     aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                <  1:470862‑M1/149‑110 (MQ=255)                                                                                                              acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                       >  2:16927‑M1/1‑31 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  >  2:73799‑M1/1‑10 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagaaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  >  2:73820‑M1/1‑10 (MQ=255)                                                                                                                                            |                                                                                                                                            GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_vch_tpiA/1978363‑1978641

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