Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	524,866	A→G	T536A (ACC→GCC)	rhsD →	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	524,866	0	A	G	100.0%	187.0 / NA	54	T536A (ACC→GCC)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/23);  total (31/23)

ATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTC  >  NC_000913/524769‑524961                                                                                                  |                                                                                                aTCGCCCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGcc                                                                                               >  2:421024/1‑100 (MQ=255)  tCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCt                                                                                              >  1:349636/1‑100 (MQ=255)   cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTg                                                                                             >  1:455688/1‑100 (MQ=255)   cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTg                                                                                             >  1:254435/1‑100 (MQ=255)   cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAAAGGAAGCGCAGGGCAGCACCCGGCAGATGGCCTg                                                                                             >  2:387916/1‑100 (MQ=255)    gcgcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTgg                                                                                            <  2:370595/100‑1 (MQ=255)     cgcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGCTGCGACGGGCAGCACCAGGCAGATGGCCTGGa                                                                                           >  2:100252/1‑100 (MQ=255)      gcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAg                                                                                          <  1:463934/100‑1 (MQ=255)      gcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAg                                                                                          >  1:279389/1‑100 (MQ=255)         gCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCg                                                                                       >  2:944145/1‑100 (MQ=255)          cGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGc                                                                                      <  2:3305/100‑1 (MQ=255)            gggAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTa                                                                                    <  2:1049132/100‑1 (MQ=255)              gagaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACg                                                                                  >  1:42923/1‑100 (MQ=255)               agaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACgg                                                                                 >  1:1090242/1‑100 (MQ=255)                gaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACggg                                                                                >  2:791224/1‑100 (MQ=255)                 aCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGc                                                                               >  1:429632/1‑100 (MQ=255)                     tACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGtt                                                                           <  1:1303511/100‑1 (MQ=255)                      aCGCTACCGCTACGATGACGCGCACAGTGAGATACCGTCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTg                                                                          <  2:881241/100‑1 (MQ=255)                       cgctaccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGc                                                                         >  2:1068904/1‑100 (MQ=255)                          taccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTgg                                                                      <  2:1297391/100‑1 (MQ=255)                          taccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTgg                                                                      <  1:998714/100‑1 (MQ=255)                           accgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGc                                                                     >  1:144490/1‑100 (MQ=255)                            ccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCg                                                                    >  2:561082/1‑100 (MQ=255)                               ctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTc                                                                 <  2:508548/100‑1 (MQ=255)                                        cgcgCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGc                                                        >  1:433192/1‑100 (MQ=255)                                           gcACAGTGAATTACCGGCGACGACAACGGATCCGACGGCCAGCACCCGGCAGATAGCATGGAGCCGCAACGGGCACTTGCTGGCGTTCACCGACTGCt                                                       >  2:379004/1‑98 (MQ=255)                                             acaGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCggg                                                   >  1:1162041/1‑100 (MQ=255)                                                  gAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTAcc                                              <  1:866068/100‑1 (MQ=255)                                                  gAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTAcc                                              >  1:1209878/1‑100 (MQ=255)                                                  gAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTAcc                                              <  1:163237/100‑1 (MQ=255)                                                    gTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAg                                            >  1:1081085/1‑100 (MQ=255)                                                       aCCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAcc                                         >  2:47084/1‑100 (MQ=255)                                                        ccGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAcc                                         <  2:429628/99‑1 (MQ=255)                                                         cGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCg                                       <  2:1264222/100‑1 (MQ=255)                                                         cGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCg                                       >  2:371474/1‑100 (MQ=255)                                                            cgacgaCAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTa                                    <  1:561086/100‑1 (MQ=255)                                                             gacgacAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGtt                                     >  2:1183928/1‑98 (MQ=255)                                                             gacgacAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTAt                                   >  2:830815/1‑100 (MQ=255)                                                              acgacAACGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATg                                  <  2:636679/100‑1 (MQ=255)                                                                 actaCGGATGCGACGGGCAGCCCCCGGCAGGTGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTAt                                   <  2:1147164/93‑1 (MQ=255)                                                                   aaCGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATAc                             <  1:945307/100‑1 (MQ=255)                                                                    aCGGATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACg                            >  1:788496/1‑100 (MQ=255)                                                                       gATGCGACGGGCTGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGAcc                         >  1:856955/1‑100 (MQ=255)                                                                       gATGCGACGGGCCGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGAcc                         >  1:856724/1‑100 (MQ=255)                                                                       gATGCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGAcc                         <  2:279386/100‑1 (MQ=255)                                                                          gCGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCt                      <  1:443008/100‑1 (MQ=255)                                                                           cGACGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCtt                     >  1:647323/1‑100 (MQ=255)                                                                              cGGGCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCgg                  >  2:282920/1‑100 (MQ=255)                                                                               gggCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGc                 <  2:1081075/100‑1 (MQ=255)                                                                                ggCAGCACCCGGCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGcc                <  1:791229/100‑1 (MQ=255)                                                                                   agcaCCCGGCAGATGGCCTGGAGCCGCTACGGTCAGTTGTTGCCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGa             <  2:915845/100‑1 (MQ=255)                                                                                          ggCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcggc      <  1:937034/100‑1 (MQ=255)                                                                                          ggCAGATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAAACCCGTTAGGAATACGACCGCTTCGGCCAGATGAcggc      >  2:710797/1‑100 (MQ=255)                                                                                              gATGGCCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTc  >  1:1162384/1‑100 (MQ=255)                                                                                                  |                                                                                                ATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTC  >  NC_000913/524769‑524961

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