Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,359,666	A→G	G122G (GGT→GGC)	cadB ←	putative lysine/cadaverine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,666	0	A	G	100.0%	166.4 / NA	47	G122G (GGT→GGC)	cadB	putative lysine/cadaverine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/20);  total (27/20)

CATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATACCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCT  >  NC_000913/4359521‑4359807                                                                                                                                                  |                                                                                                                                              cATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCccgc                                                                                                                                          >  1:386574/1‑151 (MQ=255)      aaCCAATGCCAGCCAACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGaa                                                                                                                                     <  2:83219/151‑1 (MQ=255)       aCCAATGCCAGCCAACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAAc                                                                                                                                    <  1:32827/151‑1 (MQ=255)            tGCCAGCCAACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGAt                                                                                                                               >  2:4522/1‑151 (MQ=255)             gccagccaACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGAt                                                                                                                               <  2:543200/150‑1 (MQ=255)                agccaACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCAtt                                                                                                                           >  1:88333/1‑151 (MQ=255)                   caacaaTAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTa                                                                                                                         >  2:231981/1‑150 (MQ=255)                       aaTAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACt                                                                                                                    >  1:416276/1‑151 (MQ=255)                       aaTAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACt                                                                                                                    >  2:423504/1‑151 (MQ=255)                         tAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTgg                                                                                                                  >  2:577743/1‑151 (MQ=255)                            cAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGga                                                                                                                <  2:374495/150‑1 (MQ=255)                             aGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGgaag                                                                                                              <  2:430492/151‑1 (MQ=255)                              gTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGgaaga                                                                                                             <  1:28640/151‑1 (MQ=255)                                                aGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAagataagata                                                                                           <  1:100621/151‑1 (MQ=255)                                                  aaCCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATAc                                                                                          >  2:627443/1‑150 (MQ=255)                                                        cACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGt                                                                                    <  1:516453/150‑1 (MQ=255)                                                                  aTAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCg                                                                         >  1:667919/1‑151 (MQ=255)                                                                      tGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGcc                                                                       >  2:450/1‑149 (MQ=255)                                                                        gTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAgg                                                                    >  1:291831/1‑150 (MQ=255)                                                                         ttAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGtt                                                                  >  1:297236/1‑151 (MQ=255)                                                                         ttAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGtt                                                                  >  1:635757/1‑151 (MQ=255)                                                                            aaCGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTAcc                                                               <  1:570109/151‑1 (MQ=255)                                                                                  ttACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATcc                                                          >  2:280646/1‑150 (MQ=255)                                                                                   tACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATcc                                                          <  2:206317/149‑1 (MQ=255)                                                                                   tACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATcc                                                          <  1:306027/149‑1 (MQ=255)                                                                                    aCCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGt                                                       >  1:637754/1‑151 (MQ=255)                                                                                    aCCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGt                                                       >  2:96951/1‑151 (MQ=255)                                                                                        aaGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGc                                                   >  2:432100/1‑151 (MQ=255)                                                                                        aaGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGc                                                   >  2:415338/1‑151 (MQ=255)                                                                                            aCCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATgg                                               >  2:470533/1‑151 (MQ=255)                                                                                              cgccgAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGta                                             <  2:3957/151‑1 (MQ=255)                                                                                                          ttACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACAcc                                 >  1:174569/1‑151 (MQ=255)                                                                                                           tACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACAcc                                 <  1:388140/150‑1 (MQ=255)                                                                                                             cAAAGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGt                              <  2:677318/151‑1 (MQ=255)                                                                                                              aaaGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTc                             <  1:422819/151‑1 (MQ=255)                                                                                                               aaGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCt                            >  2:108348/1‑151 (MQ=255)                                                                                                               aaGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCt                            >  1:394686/1‑151 (MQ=255)                                                                                                                aGGTAAATACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTg                           <  2:165346/151‑1 (MQ=255)                                                                                                                      aTACCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAAcc                     >  2:62659/1‑151 (MQ=255)                                                                                                                        aCCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCaa                   <  2:247301/151‑1 (MQ=255)                                                                                                                        aCCCAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCa                    <  1:402415/150‑1 (MQ=255)                                                                                                                         cccAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCa                    <  1:225148/149‑1 (MQ=255)                                                                                                                          ccAGACGATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAAt                 >  1:361495/1‑151 (MQ=255)                                                                                                                                gATAGCAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAggg           >  2:581899/1‑151 (MQ=255)                                                                                                                                     cAATACAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAAt      >  2:253811/1‑151 (MQ=255)                                                                                                                                        tACGGGCGATGCCCACCGGAACAGGATCATTTAATACCGGGAAGAAGGCGGAAAGATAAGATACAGCGTTAATGCCAATCGCCAGGTTACCAATCCAGTTAGAATGGTAATAAAGAACACCTGCTTGAAAGCCAAAGGCAGGGGGAAtgt    <  2:539601/150‑3 (MQ=255)                                                                                                                                         aCAGGCGATGCCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTtct  >  2:146546/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                              CATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCACAGGAATAAGAACCAGCACCAGACCAATAGTGGTTAAACGGCTTACCCAAGTACCGCCGAGCATATTTACAAAGGTAAATACCCAGACGATAGCAATACAGGCGATACCCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCT  >  NC_000913/4359521‑4359807

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