Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,769,486	G→A	*274* (TAG→TAA)	yfjQ →	CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,769,486	0	G	A	100.0%	189.7 / NA	59	*274* (TAG→TAA)	yfjQ	CP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/29);  total (30/29)

AAGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAA  >  NC_000913/2769352‑2769606                                                                                                                                       |                                                                                                                         aaGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGt                                                                                                                      <  2:389147/139‑1 (MQ=255)        gACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGt                                                                                                               >  2:371464/1‑139 (MQ=255)            cAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                           <  2:555720/139‑1 (MQ=255)              gTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                           <  2:984587/137‑1 (MQ=255)              gTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                           >  1:984587/1‑137 (MQ=255)               tGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATg                                                                                                                >  2:493019/1‑131 (MQ=255)               tGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATg                                                                                                                <  1:493019/131‑1 (MQ=255)               tGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        >  2:1265259/1‑139 (MQ=255)               tGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        <  1:1265259/139‑1 (MQ=255)                           tAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCg                                                                                            <  1:371464/139‑1 (MQ=255)                                  ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                     >  2:981064/1‑139 (MQ=255)                                  ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGACAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                     >  1:1357/1‑139 (MQ=255)                                          tCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCg                                                                             <  2:730926/139‑1 (MQ=255)                                             tACTCGTGACGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAAAACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                          <  2:956380/139‑1 (MQ=255)                                                  gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        >  2:1532397/1‑104 (MQ=255)                                                  gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        <  1:1532397/104‑1 (MQ=255)                                                  gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTa                                                                     >  1:331808/1‑139 (MQ=255)                                                   tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                         >  1:1497686/1‑102 (MQ=255)                                                   tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                         <  2:1497686/102‑1 (MQ=255)                                                   tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCg                                                                                            <  2:1422447/115‑1 (MQ=255)                                                   tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCg                                                                                            >  1:1422447/1‑115 (MQ=255)                                                        tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        <  1:927285/98‑1 (MQ=255)                                                        tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        >  2:927285/1‑98 (MQ=255)                                                        tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCt                                                               >  2:891723/1‑139 (MQ=255)                                                           gtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCgt                                                            <  2:739543/139‑1 (MQ=255)                                                             ggCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCgtgt                                                          >  1:217211/1‑139 (MQ=255)                                                              gCATTGACGGGGATGTGAAGCTTAACCGTGCGCCGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGtgtg                                                         >  2:612232/1‑139 (MQ=255)                                                                    aCGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCaa                                                   >  1:1506546/1‑139 (MQ=255)                                                                                aGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTtat                                       <  2:331808/139‑1 (MQ=255)                                                                                aGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTtat                                       <  1:891723/139‑1 (MQ=255)                                                                                   ttAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        >  1:281816/1‑71 (MQ=255)                                                                                   ttAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                        <  2:281816/71‑1 (MQ=255)                                                                                          tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                          >  1:392343/1‑94 (MQ=255)                                                                                          tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                          <  2:392343/94‑1 (MQ=255)                                                                                          tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                             >  2:1198038/1‑139 (MQ=255)                                                                                           gcgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAcc                            <  2:1029533/139‑1 (MQ=255)                                                                                            cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                               <  2:869866/71‑1 (MQ=255)                                                                                            cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                               >  1:869866/1‑71 (MQ=255)                                                                                                 tgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGtt                      >  2:682188/1‑139 (MQ=255)                                                                                                    gTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGa                                                                                          <  2:482623/68‑1 (MQ=255)                                                                                                    gTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGa                                                                                          >  1:482623/1‑68 (MQ=255)                                                                                                     tGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACAt                  <  2:217211/139‑1 (MQ=255)                                                                                                         ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                          >  2:33621/1‑79 (MQ=255)                                                                                                         ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                          <  1:33621/79‑1 (MQ=255)                                                                                                         ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACAt                  <  1:566386/135‑1 (MQ=255)                                                                                                         ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACAt                  >  2:566386/1‑135 (MQ=255)                                                                                                           cAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                           <  2:744159/44‑1 (MQ=255)                                                                                                           cAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                           >  1:744159/1‑44 (MQ=255)                                                                                                              aaaCACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgt                                                                                                       >  1:706332/1‑45 (MQ=255)                                                                                                              aaaCACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgt                                                                                                       <  2:706332/45‑1 (MQ=255)                                                                                                                   cTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                         >  1:620270/1‑38 (MQ=255)                                                                                                                   cTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                         <  2:620270/38‑1 (MQ=255)                                                                                                                     gCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGaaa  >  2:138691/1‑139 (MQ=255)                                                                                                                           gCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccccc                                             >  2:1344526/1‑90 (MQ=255)                                                                                                                           gCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccccc                                             <  1:1344526/90‑1 (MQ=255)                                                                                                                            cAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCa                                                                           <  2:753924/59‑1 (MQ=255)                                                                                                                            cAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCa                                                                           >  1:753924/1‑59 (MQ=255)                                                                                                                                   aGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCt                                                      >  2:1295824/1‑73 (MQ=255)                                                                                                                                   aGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCt                                                      <  1:1295824/73‑1 (MQ=255)                                                                                                                                       |                                                                                                                         AAGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAA  >  NC_000913/2769352‑2769606

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