Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	773,267	C→T	P76S (CCG→TCG)	cydB →	cytochrome d terminal oxidase, subunit II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	773,267	0	C	T	100.0%	124.5 / NA	33	P76S (CCG→TCG)	cydB	cytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/16);  total (17/16)

GACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGA  >  NC_000913/773126‑773400                                                                                                                                              |                                                                                                                                      gACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCt                                                                                                                              >  2:268733/1‑151 (MQ=255)  aCATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTa                                                                                                                             <  1:180166/151‑1 (MQ=255)      gggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGcc                                                                                                                         >  2:596981/1‑151 (MQ=255)      gggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGcc                                                                                                                         >  2:132659/1‑151 (MQ=255)                aTGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTtct                                                                                                               >  2:474840/1‑151 (MQ=255)                  gCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTtctc                                                                                                              <  2:158122/150‑1 (MQ=255)                       cccGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCtt                                                                                                        >  1:13934/1‑151 (MQ=255)                              ccTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTAtgtg                                                                                                 >  1:136858/1‑151 (MQ=255)                              ccTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTAtgtg                                                                                                 >  2:574852/1‑151 (MQ=255)                              ccTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTAtgtg                                                                                                 <  2:378506/151‑1 (MQ=255)                                     cGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCgatga                                                                                          <  2:113796/151‑1 (MQ=255)                                           gccacCGAGCGTCGAATTATGATTAACTCCATTGCTCCACACTGGGACGTTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCCCTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCg                                                                                    <  1:379181/149‑1 (MQ=255)                                                    cGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGt                                                                           >  2:62869/1‑151 (MQ=255)                                                     gTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTc                                                                          >  1:455897/1‑151 (MQ=255)                                                             aTGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGttct                                                                  >  2:407580/1‑151 (MQ=255)                                                               gATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGttctt                                                                 <  2:619061/150‑1 (MQ=255)                                                                aTTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTtcc                                                               >  2:607269/1‑151 (MQ=255)                                                                    aCTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTtccgtcc                                                           <  1:371798/151‑1 (MQ=255)                                                                     cTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTtccgtccg                                                          >  1:137302/1‑151 (MQ=255)                                                                       ccATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTtccgtccg                                                          <  2:297717/149‑1 (MQ=255)                                                                       ccATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCcggt                                                        <  1:389815/151‑1 (MQ=255)                                                                        cATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCcggtc                                                       <  2:443292/151‑1 (MQ=255)                                                                              accacACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGtttt                                                 >  2:460826/1‑151 (MQ=255)                                                                                 acacTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTACGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGAc                                              >  1:666865/1‑151 (MQ=255)                                                                                              aaCCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGa                                 <  1:106529/151‑1 (MQ=255)                                                                                                 cAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTg                              >  2:176241/1‑151 (MQ=255)                                                                                                       tGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTgaagaa                         <  1:598716/150‑1 (MQ=255)                                                                                                         gCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAAcc                      <  2:484680/151‑1 (MQ=255)                                                                                                                 ccGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCg              <  1:168066/151‑1 (MQ=255)                                                                                                                 ccGCGGGCGGCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCg              >  2:444437/1‑151 (MQ=255)                                                                                                                     gggcggCGCACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAAc          >  2:50699/1‑151 (MQ=255)                                                                                                                            gcACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTggg   <  1:293986/151‑1 (MQ=255)                                                                                                                             cACTCTTTGCTGCCTGGTCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGa  <  1:261754/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                      GACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGA  >  NC_000913/773126‑773400

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