Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,042,457	C→A	G30G (GGC→GGA)	yodB →	cytochrome b561 homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,042,457	0	C	A	100.0%	96.0 / NA	28	G30G (GGC→GGA)	yodB	cytochrome b561 homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/18);  total (10/18)

AATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTG  >  NC_000913/2042314‑2042584                                                                                                                                                |                                                                                                                                aaTTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGtt                                                                                                                            <  5:136859/149‑1 (MQ=255)   tttttAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAATGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTc                                                                                                                          <  4:134646/149‑1 (MQ=255)     tttAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTcc                                                                                                                         <  3:29669/148‑1 (MQ=255)         aGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTaaa                                                                                                                     <  2:37519/148‑1 (MQ=255)         aGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTaaa                                                                                                                     <  3:85567/148‑1 (MQ=255)              cGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                <  4:85499/148‑1 (MQ=255)                 gggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtagtag                                                                                                             >  4:37009/1‑148 (MQ=255)                                   aaaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                <  4:75246/127‑1 (MQ=255)                                    aaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCgag                                                                                         <  5:69074/149‑1 (MQ=255)                                        acaGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGaaa                                                                                      <  2:105678/148‑1 (MQ=255)                                               cTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGCTTCCTaaa                                                                                                                     <  5:112772/110‑1 (MQ=255)                                               cTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGATTCCTaaa                                                                                                                     >  6:112772/1‑110 (MQ=255)                                                 aCAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACaa                                                                            >  5:34489/1‑149 (MQ=255)                                                    aaaaatcaCCGATTCTCAAAAACTCAAATTTATTTACCTTGGATAACGCTGCTTTTCGTTGCAATAACCAATGCCGCTATGGAACTCCGTGGAAGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGc                                                                         <  8:97017/142‑1 (MQ=255)                                                                   tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTg                                                          >  2:165948/1‑149 (MQ=255)                                                                   tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTg                                                          >  7:52465/1‑149 (MQ=255)                                                                     aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTggg                                                        >  1:23049/1‑149 (MQ=255)                                                                                    tttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATctgc                                                                                                  <  4:32122/92‑1 (MQ=255)                                                                                    tttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATctgc                                                                                                  >  3:32122/1‑92 (MQ=255)                                                                                    tttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCAc                                         <  8:52465/149‑1 (MQ=255)                                                                                    tttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGGTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCa                                          <  2:134219/148‑1 (MQ=255)                                                                                        cATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCt                                     <  2:23049/149‑1 (MQ=255)                                                                                            ggATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTgg                                                         >  2:94880/1‑125 (MQ=255)                                                                                              aTAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATtat                               <  1:165948/149‑1 (MQ=255)                                                                                                        ttttCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAtt                                                                                >  2:89600/1‑90 (MQ=255)                                                                                                        ttttCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAtt                                                                                <  1:89600/90‑1 (MQ=255)                                                                                                             gTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAg                <  3:37009/149‑1 (MQ=255)                                                                                                                           tGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGTGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATtg  >  5:102911/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                AATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTG  >  NC_000913/2042314‑2042584

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