Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,227,278	A→G	S65S (AGT→AGC)	pliG ←	inhibitor of g‑type lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,227,278	0	A	G	100.0%	157.6 / NA	43	S65S (AGT→AGC)	pliG	inhibitor of g‑type lysozyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/19);  total (24/19)

CTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAG  >  NC_000913/1227142‑1227413                                                                                                                                         |                                                                                                                                        cTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTc                                                                                                                                      >  1:725145/1‑140 (MQ=255)    gCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTcaca                                                                                                                                   >  1:710690/1‑140 (MQ=255)           atatTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACtt                                                                                                                            >  1:257550/1‑140 (MQ=255)                           gggAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCtttttggca                                                                                                            >  2:628862/1‑136 (MQ=255)                            ggAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATTAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCa                                                                                                           <  2:409746/140‑1 (MQ=255)                              agtgagtATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACC‑TTTTTGGCATAg                                                                                                        <  1:628862/140‑1 (MQ=255)                               gtgagtATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAg                                                                                                        <  1:688993/140‑1 (MQ=255)                               gtgagtATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAg                                                                                                        <  1:281876/140‑1 (MQ=255)                                tgagtATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGa                                                                                                       <  2:725145/140‑1 (MQ=255)                                  agtATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGa                                                                                                       <  2:548549/138‑1 (MQ=255)                                  agtATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGa                                                                                                       >  1:548549/1‑138 (MQ=255)                                         tCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATg                                                                                              >  1:771159/1‑140 (MQ=255)                                          ccGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTatat                                                                                               >  1:174051/1‑138 (MQ=255)                                          ccGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTatat                                                                                               <  2:174051/138‑1 (MQ=255)                                             tGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTatc                                                                                          >  1:397509/1‑140 (MQ=255)                                                           tGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATCAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCt                                                                            >  1:553886/1‑140 (MQ=255)                                                             gggAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCttt                                                                          >  2:466327/1‑140 (MQ=255)                                                                 atatCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTAttt                                                                      <  1:304217/140‑1 (MQ=255)                                                                 atatCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTAttt                                                                      <  2:771159/140‑1 (MQ=255)                                                                 atatCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTAttt                                                                      <  1:105518/140‑1 (MQ=255)                                                                             aTCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACt                                                          >  1:55572/1‑140 (MQ=255)                                                                             aTCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACt                                                          >  1:658232/1‑140 (MQ=255)                                                                             aTCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTCATTTCGCCGGAGTACt                                                          >  1:624087/1‑140 (MQ=255)                                                                                aaCTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAg                                                       <  2:397509/140‑1 (MQ=255)                                                                                 aCTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAgc                                                      <  2:710690/140‑1 (MQ=255)                                                                                  cTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCGCTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAgcg                                                     <  2:658232/140‑1 (MQ=255)                                                                                  cTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTCATTTCGCCGGAGTACTGAgcg                                                     <  2:624087/140‑1 (MQ=255)                                                                                       tcatcGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCgctgct                                                >  2:40684/1‑140 (MQ=255)                                                                                           cGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATgg                                            >  1:492055/1‑140 (MQ=255)                                                                                            gATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGc                                           <  2:492055/140‑1 (MQ=255)                                                                                               tCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCtt                                        >  1:279978/1‑140 (MQ=255)                                                                                                  tGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGT‑ATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCt                                    >  2:78529/1‑140 (MQ=255)                                                                                                       cGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGT‑ATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACt                               <  1:78529/140‑1 (MQ=255)                                                                                                         aaGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTc                              <  2:369304/140‑1 (MQ=255)                                                                                                         aaGAGGTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTc                              >  1:369304/1‑140 (MQ=255)                                                                                                              gTAAGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACAt                         >  1:520969/1‑140 (MQ=255)                                                                                                                 aGTATCAGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGa                      <  1:466327/140‑1 (MQ=255)                                                                                                                       aGCGCCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCt                >  2:594766/1‑140 (MQ=255)                                                                                                                        gcgcCTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCtt               >  2:392627/1‑140 (MQ=255)                                                                                                                           ccTTCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTAcc            >  1:125226/1‑140 (MQ=255)                                                                                                                             ttCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCtt               <  2:226481/135‑1 (MQ=255)                                                                                                                             ttCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCtt               >  1:226481/1‑135 (MQ=255)                                                                                                                             ttCATTAGAGATGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCg          >  1:282947/1‑140 (MQ=255)                                                                                                                                     agaTGCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAg  >  1:663336/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        CTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAG  >  NC_000913/1227142‑1227413

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