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%	209.7 / NA	51	P76Q (CCG→CAG)	yihQ	sulfoquinovosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/28);  total (23/28)

TACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAT  >  NZ_CP009273/4062288‑4062542                                                                                                                           |                                                                                                                                     tACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGAt                                                                                                                      >  2:855606/1‑139 (MQ=255)      ttGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGt                                                                                                                 >  1:160309/1‑139 (MQ=255)      ttGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGt                                                                                                                 >  1:394221/1‑139 (MQ=255)      ttGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGt                                                                                                                 >  2:469072/1‑139 (MQ=255)         tgGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCgg                                                                                                              >  1:807153/1‑139 (MQ=255)           gTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGt                                                                                                             >  1:168668/1‑138 (MQ=255)                ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCCACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                        >  2:907774/1‑138 (MQ=255)                ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                        >  2:372026/1‑138 (MQ=255)                ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                        >  1:554468/1‑138 (MQ=255)                     tCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAAtt                                                                                                   <  1:402316/138‑1 (MQ=255)                         tttGTAGTTCCATCAATAAACGCACCTCATCTTCGACAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCACAATTTtctc                                                                                              <  2:763554/139‑1 (MQ=255)                            gTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCCACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCt                                                                                            <  1:907774/138‑1 (MQ=255)                                       aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                <  2:251838/139‑1 (MQ=255)                                       aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                <  1:708314/139‑1 (MQ=255)                                                        gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcg                                                               <  2:848971/139‑1 (MQ=255)                                                          cgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgcg                                                             >  1:485696/1‑139 (MQ=255)                                                           ggCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgcg                                                             <  1:552490/138‑1 (MQ=255)                                                            gCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCgg                                                            >  1:476083/1‑138 (MQ=255)                                                                   tatTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                    <  2:674119/139‑1 (MQ=255)                                                                       cAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGat                                                 <  1:818145/138‑1 (MQ=255)                                                                       cAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGat                                                 <  2:550942/138‑1 (MQ=255)                                                                       cAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGat                                                 <  2:476083/138‑1 (MQ=255)                                                                        aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                <  1:532508/138‑1 (MQ=255)                                                                        aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                <  2:143784/138‑1 (MQ=255)                                                                        aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                <  2:211896/138‑1 (MQ=255)                                                                          cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCg                                                                                                                  <  1:112148/70‑1 (MQ=255)                                                                          cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCg                                                                                                                  >  2:112148/1‑70 (MQ=255)                                                                          cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGat                                                 >  1:106635/1‑135 (MQ=255)                                                                          cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGat                                                 <  2:106635/135‑1 (MQ=255)                                                                            tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                              >  2:516836/1‑104 (MQ=255)                                                                            tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                              <  1:516836/104‑1 (MQ=255)                                                                              gcgcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCt                                          <  1:714827/138‑1 (MQ=255)                                                                                         aGAACCACGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                <  2:832464/89‑1 (MQ=255)                                                                                         aGAACCACGGCTGAAATGAATTAACCAACCATACTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                >  1:832464/1‑89 (MQ=255)                                                                                                cGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                >  2:371991/1‑114 (MQ=255)                                                                                                cGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                <  1:371991/114‑1 (MQ=255)                                                                                                cGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                         >  2:901853/1‑137 (MQ=255)                                                                                                cGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                         <  1:901853/137‑1 (MQ=255)                                                                                                cGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCAt                        >  1:568058/1‑138 (MQ=255)                                                                                                cGGCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCAt                        >  1:767260/1‑138 (MQ=255)                                                                                                  gCTGAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATaaa                     <  1:658493/139‑1 (MQ=255)                                                                                                     gAAATGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAAc                    >  2:742564/1‑137 (MQ=255)                                                                                                         tGAATTAACCAACCATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGAt              <  2:56260/139‑1 (MQ=255)                                                                                                               aaccaaccATCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCt         <  1:565252/138‑1 (MQ=255)                                                                                                                   aaccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCt    <  2:568058/139‑1 (MQ=255)                                                                                                                    accaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTa   <  2:55775/139‑1 (MQ=255)                                                                                                                     ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                        <  2:657247/37‑1 (MQ=255)                                                                                                                     ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                        >  1:657247/1‑37 (MQ=255)                                                                                                                     ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                <  2:593647/93‑1 (MQ=255)                                                                                                                     ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGata                                                >  1:593647/1‑93 (MQ=255)                                                                                                                     ccaTCCTGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt  >  2:328218/1‑139 (MQ=255)                                                                                                                           |                                                                                                                                     TACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAT  >  NZ_CP009273/4062288‑4062542

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