Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,917,504	Δ28 bp	coding (735‑762/864 nt)	atpG ←	ATP synthase F1 complex subunit gamma

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3917504	3917531	28	42 [0]	[0] 42	atpG	ATP synthase F1 complex subunit gamma

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3917503	0 (0.000)	39 (1.130)	37/154	0.1	100%	coding (763/864 nt)	atpG	ATP synthase F1 complex subunit gamma
	?	NC_000913	3917532 =	0 (0.000)					coding (734/864 nt)	atpG	ATP synthase F1 complex subunit gamma

GGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCATTGTCGGTCGCGGCTTTCATCGCCACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGG  >  NC_000913/3917453‑3917612                                                                                |                                                                                 cgatctcggtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccac                                                                                 <  2:842312‑M2/2‑1 (MQ=255)  gatctcggtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccac                                                                                 >  2:842361‑M2/79‑80 (MQ=255)    tctcggtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccacca                                                                               >  2:842368‑M2/77‑80 (MQ=255)        ggtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATa                                                                             <  2:842313‑M2/6‑1 (MQ=255)        ggtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACgg                                                                          >  2:842358‑M2/73‑81 (MQ=255)         gtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACggg                                                                         >  2:842372‑M2/72‑81 (MQ=255)         gtgagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACggg                                                                         <  2:842314‑M2/10‑1 (MQ=255)           gagttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGggcg                                                                       >  2:842367‑M2/70‑81 (MQ=255)            agttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGggcgg                                                                      >  1:842356‑M2/69‑81 (MQ=255)             gttcctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGggcggc                                                                     <  1:842315‑M2/14‑1 (MQ=255)                 ctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGc                                                                 >  2:842373‑M2/64‑81 (MQ=255)                 ctgagtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGc                                                                 >  1:842365‑M2/64‑81 (MQ=255)                    agtaatgctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTc                                                               <  2:842316‑M2/20‑1 (MQ=255)                           ctggcctgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGcc                                                        <  1:842317‑M2/27‑1 (MQ=255)                                 tgacgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGtttt                                                 >  1:842357‑M2/48‑81 (MQ=255)                                   acgagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTc                                                >  1:842355‑M2/46‑80 (MQ=255)                                     gagctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAAc                                             >  1:842371‑M2/44‑81 (MQ=255)                                      agctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAAcc                                            <  2:842318‑M2/39‑1 (MQ=255)                                      agctttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAAc                                             >  2:842362‑M2/43‑80 (MQ=255)                                       gctttgttgtatatcaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCa                                           <  2:842319‑M2/40‑1 (MQ=255)                                          ttgttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGc                                        >  1:842390‑M2/39‑81 (MQ=255)                                            gttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGccc                                      >  2:842366‑M2/37‑81 (MQ=255)                                             ttgtataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCt                                     >  1:842381‑M2/36‑81 (MQ=255)                                                tataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGa                                   <  1:842328‑M2/48‑1 (MQ=255)                                                tataccaactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGAt                                  <  1:842320‑M2/49‑1 (MQ=255)                                                  taccaactgcagctctttaataaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATaa                                <  2:842323‑M2/51‑1 (MQ=255)                                                      aactgcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCt                            >  1:842379‑M2/27‑81 (MQ=255)                                                          gcagctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCt                            >  1:842376‑M2/23‑77 (MQ=255)                                                            agctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGAtt                      >  1:842370‑M2/21‑81 (MQ=255)                                                            agctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGAtt                      <  1:842321‑M2/61‑1 (MQ=255)                                                             gctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTc                     <  2:842322‑M2/62‑1 (MQ=255)                                                              ctctttaatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCg                    <  2:842326‑M2/63‑1 (MQ=255)                                                                   taatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACAt                <  2:842325‑M2/67‑1 (MQ=255)                                                                   taatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACAt                <  1:842327‑M2/67‑1 (MQ=255)                                                                    aatcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATaa              >  2:842386‑M2/13‑81 (MQ=255)                                                                     atcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAac             >  2:842394‑M2/12‑81 (MQ=255)                                                                     atcaggctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAac             <  2:842324‑M2/70‑1 (MQ=255)                                                                           ctgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAacgacg         >  1:842387‑M2/6‑79 (MQ=255)                                                                            tgccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACgcag      >  2:842375‑M2/5‑81 (MQ=255)                                                                             gccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACgcagc     >  1:842389‑M2/4‑81 (MQ=255)                                                                              ccACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACgcagca    <  2:842329‑M2/79‑1 (MQ=255)                                                                                accaTACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAgg  >  1:842392/1‑81 (MQ=255)                                                                                |                                                                                 GGCCTGACGAGCTTTGTTGTATACCAACTGCAGCTCTTTAATCAGGCTGCCGCCATTGTCGGTCGCGGCTTTCATCGCCACCATACGGGCGGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGG  >  NC_000913/3917453‑3917612

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