Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,612,289	G→A	*54* (TAG→TAA)	ytjA →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,612,289	0	G	A	100.0%	201.6 / NA	62	*54* (TAG→TAA)	ytjA	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/32);  total (30/32)

ATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGA  >  NC_000913/4612158‑4612422                                                                                                                                    |                                                                                                                                      aTCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCa                                                                                                                                 >  2:336372/1‑138 (MQ=255) aTCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCa                                                                                                                                 <  1:336372/138‑1 (MQ=255)                gccgcACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAAtt                                                                                                                 >  2:849439/1‑139 (MQ=255)                           gTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGc                                                                                                      <  1:350687/139‑1 (MQ=255)                           gTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGc                                                                                                      <  1:1241468/139‑1 (MQ=255)                             ttGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCa                                                                                                    <  1:849439/139‑1 (MQ=255)                                     cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCg                                                                                            >  2:455252/1‑139 (MQ=255)                                       ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          <  1:455252/139‑1 (MQ=255)                                        gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCa                                                                                                                                 >  1:1408235/1‑99 (MQ=255)                                        gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCa                                                                                                                                 <  2:1408235/99‑1 (MQ=255)                                                  cTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                        >  1:1138543/1‑130 (MQ=255)                                                  cTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                        <  2:1138543/130‑1 (MQ=255)                                                  cTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                               >  2:498454/1‑139 (MQ=255)                                                       ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGttt                                                                          >  1:753743/1‑139 (MQ=255)                                                           cAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                      <  2:1416722/139‑1 (MQ=255)                                                            aGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAAtt                                                                                                                 >  1:1628133/1‑95 (MQ=255)                                                            aGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAAtt                                                                                                                 <  2:1628133/95‑1 (MQ=255)                                                             gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                      <  1:705867/137‑1 (MQ=255)                                                             gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                      >  2:705867/1‑137 (MQ=255)                                                             gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGa                                                                    >  2:287980/1‑139 (MQ=255)                                                                     tGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                                            <  1:498454/139‑1 (MQ=255)                                                                        cTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                         >  1:1719977/1‑139 (MQ=255)                                                                           tGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          >  1:845303/1‑103 (MQ=255)                                                                           tGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          <  2:845303/103‑1 (MQ=255)                                                                            gtcgtcGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATCGCGtt                                                     >  1:397544/1‑139 (MQ=255)                                                                                      ttattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                           <  1:287980/139‑1 (MQ=255)                                                                                              ttCCTGGTTATTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGATTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcac                                   <  2:753743/139‑1 (MQ=255)                                                                                              ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                        <  2:483636/86‑1 (MQ=255)                                                                                              ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                        >  1:483636/1‑86 (MQ=255)                                                                                              ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTg                                                                                <  2:47991/94‑1 (MQ=255)                                                                                              ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTg                                                                                >  1:47991/1‑94 (MQ=255)                                                                                                ccTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                         <  1:1079802/115‑1 (MQ=255)                                                                                                ccTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                         >  2:1079802/1‑115 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCa                                                                                                                                 >  1:679318/1‑40 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCa                                                                                                                                 <  2:679318/40‑1 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAAttt                                                                                                                <  2:440670/57‑1 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAAttt                                                                                                                >  1:440670/1‑57 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                                            <  2:1491373/109‑1 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                                            >  1:1491373/1‑109 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                     >  1:333025/1‑132 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                     <  2:333025/132‑1 (MQ=255)                                                                                                   ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATCTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCAcaacaa                              >  1:207845/1‑139 (MQ=255)                                                                                                    gTGAGTTTGTTCATGGGCCGAAAACTACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTg                                                                                <  2:1132929/88‑1 (MQ=255)                                                                                                    gTGAGTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACTATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTg                                                                                >  1:1132929/1‑88 (MQ=255)                                                                                                        gTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                  <  1:454456/130‑1 (MQ=255)                                                                                                        gTTTGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                  >  2:454456/1‑130 (MQ=255)                                                                                                          ttGTTCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAgg                       <  2:1719977/139‑1 (MQ=255)                                                                                                             ttCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          >  2:771218/1‑69 (MQ=255)                                                                                                             ttCATGGGCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          <  1:771218/69‑1 (MQ=255)                                                                                                                  gggCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                        >  1:1102714/1‑66 (MQ=255)                                                                                                                  gggCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                        <  2:1102714/66‑1 (MQ=255)                                                                                                                   ggCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAg                                                                                              <  1:263009/59‑1 (MQ=255)                                                                                                                   ggCCGAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAg                                                                                              >  2:263009/1‑59 (MQ=255)                                                                                                                       gAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                      <  2:1221234/79‑1 (MQ=255)                                                                                                                       gAAAACGACCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                      >  1:1221234/1‑79 (MQ=255)                                                                                                                              aCCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          <  2:1638719/52‑1 (MQ=255)                                                                                                                              aCCCTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                          >  1:1638719/1‑52 (MQ=255)                                                                                                                               cccTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTAc                                                   >  1:498024/1‑90 (MQ=255)                                                                                                                               cccTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTAc                                                   <  2:498024/90‑1 (MQ=255)                                                                                                                               cccTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAgg                       >  2:443555/1‑118 (MQ=255)                                                                                                                               cccTAAATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAgg                       <  1:443555/118‑1 (MQ=255)                                                                                                                               cccTAAATTTCAGACCATTAAGAAAATCTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGa  <  2:207845/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      ATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGA  >  NC_000913/4612158‑4612422

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