Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,777,363	T→C	F83F (TTT→TTC)	yfjZ →	CP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,777,363	0	T	C	100.0%	248.0 / NA	68	F83F (TTT→TTC)	yfjZ	CP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/27);  total (41/27)

TCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAAC  >  NC_000913/2777270‑2777458                                                                                              |                                                                                                tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGc                                                                                           >  2:1444976/1‑100 (MQ=255)  ttccgatctATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCg                                                                                          <  2:467734/92‑1 (MQ=37)        tttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGACCTGAATCCCCGCCATGCCCACTGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCg                                                                                    <  2:707571/100‑1 (MQ=255)        tttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCg                                                                                    <  1:208928/100‑1 (MQ=255)         ttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGa                                                                                   <  2:1346301/100‑1 (MQ=255)           aTCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGagatcggaa                                                                                 >  1:467735/1‑92 (MQ=37)            tCAGCCAGATAGAGTCGAGGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGACACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATAc                                                                                <  2:1161853/100‑1 (MQ=255)             cAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATAct                                                                               >  2:562038/1‑100 (MQ=255)               gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCACTGCGTCACCATGTACAACAACGGTATCACCTGCGAAGCCGATActct                                                                             >  2:1067937/1‑100 (MQ=255)               gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATActct                                                                             >  2:884969/1‑100 (MQ=255)                ccAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCtt                                                                            <  2:798303/100‑1 (MQ=255)                ccAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCtt                                                                            >  2:286293/1‑100 (MQ=255)                 cAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTg                                                                           <  1:1203712/100‑1 (MQ=255)                 cAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTg                                                                           <  2:1427953/100‑1 (MQ=255)                     tAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGTCATGCCCAATGCGTCACTCTGTATCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGt                                                                         <  2:546591/98‑1 (MQ=255)                     tAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAg                                                                       >  2:309210/1‑100 (MQ=255)                       gcgTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGtt                                                                     >  1:1073734/3‑100 (MQ=255)                        agTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTg                                                                    >  1:1103129/1‑100 (MQ=255)                        agTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTg                                                                    >  2:90440/1‑100 (MQ=255)                         gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGc                                                                   >  1:950664/1‑100 (MQ=255)                                 tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCGCAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGt                                                           <  2:1197739/100‑1 (MQ=255)                                  gACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTCCGta                                                          >  2:471856/1‑100 (MQ=255)                                  gACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTAAGta                                                          >  1:1476954/1‑100 (MQ=255)                                   aCCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                         <  2:1570896/100‑1 (MQ=255)                                   aCCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtat                                                         <  1:825664/100‑1 (MQ=255)                                    ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCTTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                        <  2:765655/100‑1 (MQ=255)                                    ccACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGtata                                                        >  2:1101187/1‑100 (MQ=255)                                      aCTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCCCCTGAGAAGCCGATCATCTTGGTAGTTGCGGCTACGtgt                                                         >  2:968997/1‑95 (MQ=21)                                      aCTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACa                                                      >  1:1330385/1‑100 (MQ=255)                                       cTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACAt                                                     >  2:420852/1‑100 (MQ=255)                                       cTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACAt                                                     >  2:1257742/1‑100 (MQ=255)                                         ggTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCg                                                   >  1:1351673/1‑100 (MQ=255)                                         ggTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCg                                                   >  1:953067/1‑100 (MQ=255)                                            gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                                >  1:1062197/1‑100 (MQ=255)                                            gAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTg                                                >  1:403834/1‑100 (MQ=255)                                             aaCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                               >  1:712393/1‑100 (MQ=255)                                             aaCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                               >  2:363560/1‑100 (MQ=255)                                              aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                              >  2:1441160/1‑100 (MQ=255)                                              aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGt                                               >  1:772212/1‑99 (MQ=255)                                               cTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                             <  1:49877/100‑1 (MQ=255)                                                  aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                          <  1:584139/100‑1 (MQ=255)                                                  aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                          >  2:927606/1‑100 (MQ=255)                                                  aaTCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcc                                          <  1:1284962/100‑1 (MQ=255)                                                     tcccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAcccc                                        <  2:959261/98‑1 (MQ=255)                                                     ccccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCa                                       <  2:955521/100‑1 (MQ=255)                                                           cATGCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAAc                                 <  1:1257747/100‑1 (MQ=255)                                                              gCCCAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGct                              >  2:593219/1‑100 (MQ=255)                                                               cccAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGcta                             >  2:712641/1‑100 (MQ=255)                                                                 cAATGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaac                           <  2:772210/100‑1 (MQ=255)                                                                  aaTGCGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaact                          <  1:980179/100‑1 (MQ=255)                                                                      cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                      >  2:25467/1‑100 (MQ=255)                                                                      cGTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAtt                      >  1:1089650/1‑100 (MQ=255)                                                                       gTCACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAttt                     >  1:1527298/1‑100 (MQ=255)                                                                         cACCCTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                    <  1:813955/99‑1 (MQ=255)                                                                            ccTGTACCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACga                >  1:251335/1‑100 (MQ=255)                                                                                 aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca           >  1:1326378/1‑100 (MQ=255)                                                                                 aCCACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagca           >  1:840581/1‑100 (MQ=255)                                                                                  ccACAACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaa          >  1:1580891/1‑100 (MQ=255)                                                                                   cacaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaag         >  1:1372340/1‑100 (MQ=255)                                                                                    acaACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagc        >  2:1326037/1‑100 (MQ=255)                                                                                     caACGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagcaagca       <  1:1326044/100‑1 (MQ=255)                                                                                      aaCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCAt      <  1:712643/100‑1 (MQ=255)                                                                                       aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATg     >  2:682272/1‑100 (MQ=255)                                                                                       aCGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATg     >  1:666737/1‑100 (MQ=255)                                                                                        cGGTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGa    <  1:309210/100‑1 (MQ=255)                                                                                         ggTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGaa   >  1:982675/1‑100 (MQ=255)                                                                                          gTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAac  <  1:1085052/100‑1 (MQ=255)                                                                                          gTTTCACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAac  <  1:997069/100‑1 (MQ=255)                                                                                              |                                                                                                TCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAAC  >  NC_000913/2777270‑2777458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑