Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,419,344	G→A	A38V (GCT→GTT)	ydaG ←	Rac prophage, uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,419,344	0	G	A	100.0%	160.8 / NA	50	A38V (GCT→GTT)	ydaG	Rac prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/27);  total (23/27)

TATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCG  >  NC_000913/1419224‑1419478                                                                                                                         |                                                                                                                                       tataCGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                      <  2:683380/139‑1 (MQ=255)    aCGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGAt                                                                                                                   >  1:626077/1‑139 (MQ=255)    aCGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGAt                                                                                                                   >  1:372804/1‑139 (MQ=255)          tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGt                                                                                                             >  1:975794/1‑139 (MQ=255)                tGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAAATCCCAGTCGTtgc                                                                                                       >  2:1299507/1‑139 (MQ=255)                  aaCGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgctg                                                                                                     >  2:158288/1‑139 (MQ=255)                          aTCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGc                                                                                             <  2:1321699/139‑1 (MQ=255)                             cTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTg                                                                                          <  2:38416/139‑1 (MQ=255)                               ggTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                        <  1:504095/139‑1 (MQ=255)                                gTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  >  1:1032615/1‑128 (MQ=255)                                gTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                  <  2:1032615/128‑1 (MQ=255)                                 ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCTATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                      <  2:1283835/107‑1 (MQ=255)                                 ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                      >  1:1283835/1‑107 (MQ=255)                                 ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGcaca                                                                                      >  2:1326137/1‑139 (MQ=255)                                 ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGcaca                                                                                      >  2:1184399/1‑139 (MQ=255)                                   ggggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                    <  1:301705/139‑1 (MQ=255)                                    gggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGcc                                                                                   <  2:502236/139‑1 (MQ=255)                                       aaCACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTgg                                                                                <  2:626077/139‑1 (MQ=255)                                         cacGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAAATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGtt                                                                              <  1:1299507/139‑1 (MQ=255)                                                     cAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGAc                                                                                                                                    >  2:413368/1‑73 (MQ=255)                                                     cAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGAc                                                                                                                                    <  1:413368/73‑1 (MQ=255)                                                           ttCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaaca                                                            <  1:1184399/139‑1 (MQ=255)                                                             cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGAc                                                                                                                                    >  1:126619/1‑65 (MQ=255)                                                             cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGAc                                                                                                                                    <  2:126619/65‑1 (MQ=255)                                                             cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAAt                                                          <  2:372804/139‑1 (MQ=255)                                                                  ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATc                                                     <  2:975794/139‑1 (MQ=255)                                                                  ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATc                                                     <  1:1326137/139‑1 (MQ=255)                                                                     ccGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGcc                                                                                   <  2:64535/106‑1 (MQ=255)                                                                     ccGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGcc                                                                                   >  1:64535/1‑106 (MQ=255)                                                                          gTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGATTTCTGCCTGCCaga                                                                                                                      <  1:817468/66‑1 (MQ=255)                                                                          gTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                      >  2:817468/1‑66 (MQ=255)                                                                             gATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                  <  1:276418/99‑1 (MQ=255)                                                                             gATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                  >  2:276418/1‑99 (MQ=255)                                                                                    tGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgc                                                                                                                              <  2:989327/48‑1 (MQ=255)                                                                                    tGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgc                                                                                                                              >  1:989327/1‑48 (MQ=255)                                                                                           tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGTGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                            <  1:158288/139‑1 (MQ=255)                                                                                           tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGaa                                           <  2:18295/124‑1 (MQ=255)                                                                                           tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGaa                                           >  1:18295/1‑124 (MQ=255)                                                                                                     acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGAACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATa                                                <  2:351575/109‑1 (MQ=255)                                                                                                     acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGAACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATa                                                >  1:351575/1‑109 (MQ=255)                                                                                                            gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGa                                                                      >  2:731272/1‑80 (MQ=255)                                                                                                            gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGa                                                                      <  1:731272/80‑1 (MQ=255)                                                                                                             ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                            <  2:1236002/121‑1 (MQ=255)                                                                                                             ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                            >  1:1236002/1‑121 (MQ=255)                                                                                                             ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCt          >  1:675641/1‑139 (MQ=255)                                                                                                                   cGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTg                                                                                          <  2:23076/53‑1 (MQ=255)                                                                                                                   cGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTg                                                                                          >  1:23076/1‑53 (MQ=255)                                                                                                                   cGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                            <  2:923890/115‑1 (MQ=255)                                                                                                                   cGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                            >  1:923890/1‑115 (MQ=255)                                                                                                                     tataACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGctg  >  2:1206285/1‑137 (MQ=255)                                                                                                                         |                                                                                                                                       TATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCG  >  NC_000913/1419224‑1419478

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