Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	448,092	Δ1 bp	coding (569/615 nt)	cyoC ←	cytochrome bo3 ubiquinol oxidase subunit 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	448,090	0	A	.	100.0%	189.2 / NA	41	coding (571/615 nt)	cyoC	cytochrome bo3 ubiquinol oxidase subunit 3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (24/17);  total (24/17)

ATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGC  >  NC_000913/448011‑448161                                                                                |                                                                        aTGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCaa                                                                        >  1:86768/1‑81 (MQ=255)  tGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCaac                                                                       <  1:86723/81‑2 (MQ=255)     gACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCaacca                                                                     >  1:86758/1‑77 (MQ=255)      aCGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACcaca                                                                  <  1:86724/81‑1 (MQ=255)       cgcgCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACAt                                                                 >  1:86787/1‑81 (MQ=255)       cgcgCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACAt                                                                 >  1:86771/1‑81 (MQ=255)        gcgcCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATc                                                                >  1:86775/1‑81 (MQ=255)          gccgcTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATc                                                                >  2:86785/1‑79 (MQ=255)          gccgcTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATc                                                                >  2:86781/1‑79 (MQ=255)          gccgcTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATc                                                                >  2:86774/1‑79 (MQ=255)          gccgcTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATc                                                                >  2:86766/1‑79 (MQ=255)           ccgcTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAg                                                             <  1:86725/81‑1 (MQ=255)                gtgATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGt                                                        >  2:86757/1‑81 (MQ=255)                 tgATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGt                                                        <  1:86726/80‑1 (MQ=255)                  gATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGc                                                      >  1:86779/1‑81 (MQ=255)                   aTCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGcc                                                     <  2:86727/81‑1 (MQ=255)                   aTCGGTAGAATGACTCATTACATCG‑CCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGcc                                                     <  1:86728/80‑1 (MQ=255)                      ggTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGa                                                  <  2:86730/81‑1 (MQ=255)                      ggTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGa                                                  <  1:86729/81‑1 (MQ=255)                       gTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGaa                                                 >  1:86788/1‑81 (MQ=255)                       gTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGaa                                                 <  2:86731/81‑1 (MQ=255)                         aGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGaa                                                 >  2:86773/1‑79 (MQ=255)                          gAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAg                                              >  2:86764/1‑81 (MQ=255)                           aaTGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAgg                                             >  1:86760/1‑81 (MQ=255)                           aaTGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAgg                                             >  1:86798/1‑81 (MQ=255)                            aTGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGc                                            <  2:86732/81‑1 (MQ=255)                              gACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCt                                           >  2:86776/1‑80 (MQ=255)                                  cATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAgg                                       >  2:86793/1‑80 (MQ=255)                                   aTTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGc                                      <  1:86733/80‑1 (MQ=255)                                      aCATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACat                                  >  1:86777/1‑81 (MQ=255)                                       cATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACatg                                 <  1:86734/81‑1 (MQ=255)                                         tCGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACatgat                               <  2:86735/81‑1 (MQ=255)                                          tcgccccATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACatgatg                              <  2:86736/78‑1 (MQ=255)                                          cGCCCCCATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACatgatg                              >  2:86782/1‑81 (MQ=255)                                             ccccATCAGATAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCgg                           >  2:86783/1‑81 (MQ=255)                                                      aTAAACAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTa                  <  2:86737/81‑1 (MQ=255)                                                        aaaCAACAGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGt                >  2:86799/1‑81 (MQ=255)                                                          acaacaGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGc              >  2:86786/1‑81 (MQ=255)                                                           caacaGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCt             <  1:86738/81‑1 (MQ=255)                                                            aacaGTGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTg            >  2:86778/1‑81 (MQ=255)                                                                 tGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAg       >  2:86792/1‑81 (MQ=255)                                                                 tGAACACACAGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAg       <  1:86739/81‑1 (MQ=255)                                                                     cacacaGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCcg   >  2:86784/1‑81 (MQ=255)                                                                      acacaGATCC‑AACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCcgc  <  1:86740/81‑1 (MQ=255)                                                                                |                                                                        ATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGC  >  NC_000913/448011‑448161

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