Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,062,410	G→T	P76Q (CCG→CAG)	yihQ ←	sulfoquinovosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,062,410	0	G	T	100.0%	136.2 / NA	34	P76Q (CCG→CAG)	yihQ	sulfoquinovosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/17);  total (17/17)

AAGGCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTA  >  NZ_CP009273/4062274‑4062541                                                                                                                                         |                                                                                                                                    aaGGCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGt                                                                                                                                   >  2:479059/1‑139 (MQ=255)  aGGCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGt                                                                                                                                   >  1:463997/1‑138 (MQ=255)          ccAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCt                                                                                                                           >  1:215104/1‑138 (MQ=255)              aTACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTACAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGa                                                                                                                      >  1:177523/1‑139 (MQ=255)                    ttGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGt                                                                                                                >  2:447700/1‑139 (MQ=255)                         gTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGt                                                                                                            >  2:198279/1‑138 (MQ=255)                            aaGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAgc                                                                                                        >  1:128294/1‑139 (MQ=255)                            aaGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAg                                                                                                         <  2:215104/138‑1 (MQ=255)                             aGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                       <  2:41155/139‑1 (MQ=255)                              ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                       >  1:412314/1‑138 (MQ=255)                                  gTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAAtt                                                                                                  >  1:524340/1‑139 (MQ=255)                                                                   aTCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAAt                                                                   <  2:558367/137‑1 (MQ=255)                                                                         ggCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgcg                                                            <  1:34387/138‑1 (MQ=255)                                                                         ggCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCgg                                                           <  1:4519/139‑1 (MQ=255)                                                                                     cAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGat                                                <  2:463997/138‑1 (MQ=255)                                                                                      aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                               <  1:343212/138‑1 (MQ=255)                                                                                      aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                               <  2:177523/138‑1 (MQ=255)                                                                                        cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgcg                                                            <  1:133890/123‑1 (MQ=255)                                                                                        cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgcg                                                            >  2:133890/1‑123 (MQ=255)                                                                                        cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                            >  2:304635/1‑139 (MQ=255)                                                                                           ggCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGc                                          >  1:229356/1‑138 (MQ=255)                                                                                            gcgcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCt                                         <  1:198786/138‑1 (MQ=255)                                                                                               cTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATAc                                     <  1:125843/139‑1 (MQ=255)                                                                                                         aaCCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAAt                           >  1:402100/1‑139 (MQ=255)                                                                                                            cACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                        <  2:93767/139‑1 (MQ=255)                                                                                                               ggCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                               >  2:25009/1‑113 (MQ=255)                                                                                                               ggCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                               <  1:25009/113‑1 (MQ=255)                                                                                                                gCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATaaa                    <  2:127297/139‑1 (MQ=255)                                                                                                                         aaTTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                               >  2:149539/1‑103 (MQ=255)                                                                                                                         aaTTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                               <  1:149539/103‑1 (MQ=255)                                                                                                                              accaaccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCttt      >  2:235008/1‑139 (MQ=255)                                                                                                                                  accaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTa  <  2:128294/139‑1 (MQ=255)                                                                                                                                   ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                       >  2:365407/1‑37 (MQ=255)                                                                                                                                   ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                       <  1:365407/37‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    AAGGCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTA  >  NZ_CP009273/4062274‑4062541

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