Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,598,557	A→G	*368Q (TAA→CAA)	livJ ←	branched‑chain amino acid ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,598,557	0	A	G	100.0%	159.2 / NA	45	*368Q (TAA→CAA)	livJ	branched‑chain amino acid ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/20);  total (25/20)

TACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTAT  >  NC_000913/3598423‑3598684                                                                                                                                       |                                                                                                                                tACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTc                                                                                                                             >  2:266028/1‑139 (MQ=255)             tGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCgcc                                                                                                                 >  2:1244569/1‑139 (MQ=255)                 aTTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCgccg                                                                                                                >  2:687884/1‑136 (MQ=255)                 aTTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCgccg                                                                                                                <  1:687884/136‑1 (MQ=255)                              tttGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCa                                                                                                >  2:127956/1‑139 (MQ=255)                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTg                                                                                                   <  1:1032293/126‑1 (MQ=255)                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTg                                                                                                   >  2:1032293/1‑126 (MQ=255)                                               ggATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACt                                                                               >  2:1204041/1‑139 (MQ=255)                                                gATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTg                                                                                                   <  2:1256330/118‑1 (MQ=255)                                                gATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTg                                                                                                   >  1:1256330/1‑118 (MQ=255)                                                 aTAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCa                                                                             >  2:188387/1‑139 (MQ=255)                                                           tACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGc                                                                                                  >  2:903211/1‑108 (MQ=255)                                                           tACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGc                                                                                                  <  1:903211/108‑1 (MQ=255)                                                           tACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGGCGAACTCAAAGCCTTTCa                                                                   >  1:904107/1‑139 (MQ=255)                                                           tACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCa                                                                   >  2:931113/1‑139 (MQ=255)                                                             cgccgcATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGa                                                                 <  1:931113/139‑1 (MQ=255)                                                                 gcATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGc                                                             <  1:1244569/139‑1 (MQ=255)                                                                  cATCCGGCATGGAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTCATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGt                                                                                                    <  1:855977/99‑1 (MQ=255)                                                                  cATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGt                                                                                                    >  2:855977/1‑99 (MQ=255)                                                                       ggCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTtct                                                       <  2:1017256/139‑1 (MQ=255)                                                                            gAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACg                                                                                      <  1:61090/103‑1 (MQ=255)                                                                            gAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACg                                                                                      >  2:61090/1‑103 (MQ=255)                                                                                   cGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGt                                             >  1:301076/1‑137 (MQ=255)                                                                                   cGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGt                                             <  2:301076/137‑1 (MQ=255)                                                                                     tACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGc                                         >  2:892265/1‑139 (MQ=255)                                                                                             atTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCgccg                                                                                                                >  2:312929/1‑60 (MQ=255)                                                                                             atTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCgccg                                                                                                                <  1:312929/60‑1 (MQ=255)                                                                                              ttAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCAtt                                <  1:188387/139‑1 (MQ=255)                                                                                               tAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTa                               >  1:933774/1‑139 (MQ=255)                                                                                                aGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGc                                                                                                  <  1:1122767/71‑1 (MQ=255)                                                                                                aGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGc                                                                                                  >  2:1122767/1‑71 (MQ=255)                                                                                                  cAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACgg                            >  1:534479/1‑139 (MQ=255)                                                                                                    atttggtagcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGa                                                                 <  2:194641/92‑1 (MQ=255)                                                                                                    atttggtagcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGa                                                                 >  1:194641/9‑100 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCTCATCGGTCGCCGTGCCGTTg                                                                                                        <  2:792815/59‑1 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTg                                                                                                        >  1:792815/1‑59 (MQ=255)                                                                                                          cggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACg                             >  2:1151385/1‑130 (MQ=255)                                                                                                          cggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACg                             <  1:1151385/130‑1 (MQ=255)                                                                                                           ggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACgg                   <  1:892265/139‑1 (MQ=255)                                                                                                                         gCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACgg                            <  2:222296/116‑1 (MQ=255)                                                                                                                         gCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACgg                            >  1:222296/1‑116 (MQ=255)                                                                                                                          ccGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAg                                                                          <  1:953892/69‑1 (MQ=255)                                                                                                                          ccGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGt    >  2:1212292/1‑139 (MQ=255)                                                                                                                          ccGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGACAGTCAAATACGCCGAACTCAAAg                                                                          >  2:953892/1‑69 (MQ=255)                                                                                                                            gATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTAt  <  2:534479/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                TACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTAT  >  NC_000913/3598423‑3598684

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