Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,622,890	C→T	S374S (AGC→AGT)	yjjJ →	toxin YjjJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,890	0	C	T	100.0%	150.0 / NA	45	S374S (AGC→AGT)	yjjJ	toxin YjjJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/24);  total (21/24)

AAGTGATCTGGGCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCG  >  NC_000913/4622754‑4623025                                                                                                                                         |                                                                                                                                        aaGTGATCTGGGCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCt                                                                                                                                      >  2:570799/1‑140 (MQ=255)       tCTGGGCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAAtg                                                                                                                                <  1:850482/140‑1 (MQ=255)         tGGGCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAAtgct                                                                                                                              >  1:828474/1‑140 (MQ=255)         tGGGCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAAtgct                                                                                                                              >  1:811821/1‑140 (MQ=255)           ggCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAAtgctgc                                                                                                                            >  2:87942/1‑140 (MQ=255)           ggCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGGCTAGGCACCAAACTGTGCTGGAAtgctgc                                                                                                                            >  1:226042/1‑140 (MQ=255)           ggCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTACCGATGGGCTAGGCACCAAACGGTGCTGGAAtgctgc                                                                                                                            >  1:261364/1‑140 (MQ=255)            gCATTTGGGCGACTTATTGCCACCAGCGATAGGCAGGCATGTAATTTACCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCg                                                                                                                           <  2:261364/140‑1 (MQ=255)            gCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCg                                                                                                                           <  2:811821/140‑1 (MQ=255)            gCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCg                                                                                                                           <  1:79876/140‑1 (MQ=255)               tttGGGCGACTTATCGCCAACAGCGATATGCACGCGGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGa                                                                                                                        >  2:706324/1‑140 (MQ=255)                 tGGGCGACTTATCGCCAACAGCTATATGCCCGCAGGTCATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGAtg                                                                                                                      <  1:223764/140‑1 (MQ=255)                     cGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCg                                                                                                                           <  2:857554/131‑1 (MQ=255)                     cGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCg                                                                                                                           >  1:857554/1‑131 (MQ=255)                            tCGCCAACAGCGATATGAACGCCGGTAATTTCTCGATTTCTTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAg                                                                                                           <  2:226042/140‑1 (MQ=255)                               ccAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTg                                                                                                        <  1:711521/140‑1 (MQ=255)                                   cAGCGATATGCACGCAGGTATTTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGt                                                                                                    <  2:578245/140‑1 (MQ=255)                                    aGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGtt                                                                                                   <  2:828474/140‑1 (MQ=255)                                    aGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGtt                                                                                                   <  2:568982/140‑1 (MQ=255)                                    aGCGATATGCACGCAGGTAATTGATCGTTTTATTTATTTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGtt                                                                                                   <  1:333617/140‑1 (MQ=255)                                              aCGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAAc                                                                                         <  1:731476/140‑1 (MQ=255)                                               cGCAGGTAATTTATCCTTTTATGTAGCTGCACCGCCAGTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACg                                                                                        <  1:305207/140‑1 (MQ=255)                                                       aTTTATCGTTTTATTTGCCTGATCCGCCATTGGGGTTGAGGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTaaa                                                                                <  2:406900/140‑1 (MQ=255)                                                              gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAgcgc                                                                           >  1:639896/1‑138 (MQ=255)                                                              gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAgcgc                                                                           <  1:170607/138‑1 (MQ=255)                                                              gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAgcgc                                                                           <  2:639896/138‑1 (MQ=255)                                                              gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAgcgc                                                                           >  2:170607/1‑138 (MQ=255)                                                              gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTCGAGCTTAACGTCAGTAAAAGCGCtt                                                                         >  1:406900/1‑140 (MQ=255)                                                                                 ccATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCt                                                      >  1:700303/1‑140 (MQ=255)                                                                                   aTTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTgg                                                    <  1:706324/140‑1 (MQ=255)                                                                                      tGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTCATCTTAACGTCAATAAAAGCGCTTGATTAACGGCGATCCCGCTGgcgc                                                 >  2:361876/1‑140 (MQ=255)                                                                                         gcTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagc                                              >  2:156583/1‑140 (MQ=255)                                                                                           tGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGGATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGcagcag                                            >  2:486490/1‑140 (MQ=255)                                                                                             aCGCCCGTCTACGACATGCTGCGGACGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGtt                                          <  1:361876/140‑1 (MQ=255)                                                                                                       aCGACATGCTGCCGATGTTCTATGCACCAAACATTTCTGGAATGCTGCCTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACg                                <  2:374455/140‑1 (MQ=255)                                                                                                         gACATGCTGCCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGt                              <  1:570799/140‑1 (MQ=255)                                                                                                              gctgcCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGt                              <  1:137798/135‑1 (MQ=255)                                                                                                              gctgcCGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGt                              >  2:137798/1‑135 (MQ=255)                                                                                                                  ccGATGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGTAGTTCTGGCAAACGGTCGCCAgaga                     >  1:374455/1‑140 (MQ=255)                                                                                                                      tGGTCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCcg                 >  1:206813/1‑140 (MQ=255)                                                                                                                         tCTATGCACCAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGt              >  2:144317/1‑140 (MQ=255)                                                                                                                                  cAAACAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTACTGGCAACCGGGCGCCAGAGAGCCGCGTATCAGCGAg     >  1:334142/1‑140 (MQ=255)                                                                                                                                    aaCAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGACGCCAGAGATCCGCGTATCAGCGAGGc   >  2:667359/1‑140 (MQ=255)                                                                                                                                     aCAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCg  <  1:144317/140‑1 (MQ=255)                                                                                                                                     aCAGTGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCg  <  1:156583/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        AAGTGATCTGGGCATTTGGGCGACTTATCGCCAACAGCGATATGCACGCAGGTAATTTATCGTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCG  >  NC_000913/4622754‑4623025

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