Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	7,878	A→G	F28L (TTC→CTC)	yaaJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	7,878	0	A	G	100.0%	206.9 / NA	57	F28L (TTC→CTC)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/29);  total (28/29)

CAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA  >  NC_000913/7768‑8011                                                                                                               |                                                                                                                                       caATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTaa                                                                                                            >  2:945073/1‑139 (MQ=255)      aCTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAgg                                                                                                               >  1:396473/1‑131 (MQ=255)      aCTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAgg                                                                                                               <  2:396473/131‑1 (MQ=255)      aCTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCAt                                                                                                       >  1:597245/1‑139 (MQ=255)      aCTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCAt                                                                                                       >  2:1050251/1‑139 (MQ=255)         gAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTAc                                                                                                    >  2:130732/1‑139 (MQ=255)          aaaTGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTAcc                                                                                                   >  2:251247/1‑139 (MQ=255)          aaaTGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTAcc                                                                                                   <  1:445255/139‑1 (MQ=255)           aaTGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCg                                                                                                  <  2:70965/139‑1 (MQ=255)                    aaaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGt                                                                                                              <  2:323676/118‑1 (MQ=255)                    aaaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGt                                                                                                              >  1:323676/1‑118 (MQ=255)                        cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaagga                                                                                    <  1:1128935/139‑5 (MQ=255)                          ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                   <  1:251247/139‑1 (MQ=255)                            tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGt                                                                                                              >  1:191097/1‑110 (MQ=255)                            tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGt                                                                                                              <  2:191097/110‑1 (MQ=255)                              gCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGt                                                                               >  1:308216/1‑139 (MQ=255)                                        aaTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCg                                                                                                                        >  2:701205/1‑88 (MQ=255)                                        aaTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCg                                                                                                                        <  1:701205/88‑1 (MQ=255)                                         aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAagag                                                                                                                   >  2:115869/1‑92 (MQ=255)                                         aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAagag                                                                                                                   <  1:115869/92‑1 (MQ=255)                                         aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGgag                                                                    >  1:933220/1‑139 (MQ=255)                                            cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaagga                                                                                    >  1:893375/1‑115 (MQ=255)                                            cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaagga                                                                                    <  2:893375/119‑5 (MQ=255)                                                    aaGACTTTTGCCAAGCTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCTCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCt                                                                                 >  2:740575/1‑115 (MQ=255)                                                    aaGACTTTTGCCAAACTGGTGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                         <  2:184788/139‑1 (MQ=255)                                                    aaGACTTTTGCCAAACTGGTGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                         >  1:184788/1‑139 (MQ=255)                                                    aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCt                                                                                 <  1:740575/115‑1 (MQ=255)                                                      gACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                       <  2:933220/139‑1 (MQ=255)                                                          tttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                   <  1:1198036/139‑1 (MQ=255)                                                          tttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                   <  2:758408/139‑1 (MQ=255)                                                           ttGCCAAACTGGGGGTTGATGCGACACTGCACAAATCCGGTGCGAAAATTGAGCTAACAACCTGCGCCGAAGCGCAGGTAGATCATTACAGATCCCCAAAAGGAAGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                  >  1:927737/1‑139 (MQ=255)                                                           ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAgg                                                                                                               >  2:1170660/1‑78 (MQ=255)                                                           ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAgg                                                                                                               <  1:1170660/78‑1 (MQ=255)                                                           ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                  <  1:707640/139‑1 (MQ=255)                                                                                cGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAgg                                                                                                               >  2:257233/1‑57 (MQ=255)                                                                                cGAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAgg                                                                                                               <  1:257233/57‑1 (MQ=255)                                                                                 gAAACTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgcg                            >  2:423775/1‑139 (MQ=255)                                                                                     cTGCACAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGa                        <  2:1077664/139‑1 (MQ=255)                                                                                          caAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                   <  1:130732/139‑1 (MQ=255)                                                                                          caAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                   <  1:1050251/139‑1 (MQ=255)                                                                                             aTCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGaa                                                                                                                      >  2:549861/1‑37 (MQ=255)                                                                                             aTCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGaa                                                                                                                      <  1:549861/37‑1 (MQ=255)                                                                                             aTCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAg                                                                       <  2:297013/84‑1 (MQ=255)                                                                                             aTCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAg                                                                       >  1:297013/1‑84 (MQ=255)                                                                                             aTCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                   >  1:1083782/1‑136 (MQ=255)                                                                                             aTCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                   <  2:1083782/136‑1 (MQ=255)                                                                                                  gTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa           >  1:337345/1‑139 (MQ=255)                                                                                                  gTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa           >  2:745030/1‑139 (MQ=255)                                                                                                    gCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaaggac                                                                                   >  2:338468/1‑59 (MQ=255)                                                                                                    gCGAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaaggac                                                                                   <  1:338468/64‑6 (MQ=255)                                                                                                      gAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCa                                                                                                        <  2:328941/42‑1 (MQ=255)                                                                                                      gAAAAGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCa                                                                                                        >  1:328941/1‑42 (MQ=255)                                                                                                       aaaaGTGAGCCAACAACCTGCGCCGAAGAGCAGGTa                                                                                                             >  2:393291/1‑36 (MQ=255)                                                                                                       aaaaGTGAGCCAACAACCTGCGCCGAAGAGCAGGTa                                                                                                             <  1:393291/36‑1 (MQ=255)                                                                                                       aaaaGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAAc         <  2:843603/136‑1 (MQ=255)                                                                                                       aaaaGTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAAc         >  1:843603/1‑136 (MQ=255)                                                                                                           gTGAGCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa  <  2:778855/139‑1 (MQ=255)                                                                                                               |                                                                                                                                       CAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCA  >  NC_000913/7768‑8011

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: ‑