Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	570,833	G→A	intergenic (+405/‑60)	ybcK → / → ybcL	DLP12 prophage; putative recombinase/DLP12 prophage; periplasmic protein YbcL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	570,833	0	G	A	100.0%	167.5 / NA	47	intergenic (+405/‑60)	ybcK/ybcL	DLP12 prophage; putative recombinase/DLP12 prophage; periplasmic protein YbcL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/17);  total (30/17)

ACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCATTTCAGGT  >  NC_000913/570697‑570963                                                                                                                                         |                                                                                                                                   acacCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACAt                                                                                                                                 >  1:686321/1‑140 (MQ=255)                    ttACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATc                                                                                                              >  2:298229/1‑140 (MQ=255)                      aCTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCaa                                                                                                            >  2:28495/1‑140 (MQ=255)                      aCTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCaa                                                                                                            >  2:1037184/1‑140 (MQ=255)                       cTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCaaa                                                                                                           >  1:515021/1‑140 (MQ=255)                        taataataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCaaat                                                                                                          <  1:298229/140‑1 (MQ=255)                            aataaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCaaataaaa                                                                                                      <  2:279154/140‑1 (MQ=255)                              taaAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCaaataaaata                                                                                                    >  2:50982/1‑140 (MQ=255)                                  gCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGa                                                                                                <  2:515021/140‑1 (MQ=255)                                  gCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGa                                                                                                <  1:1037184/140‑1 (MQ=255)                                  gCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGa                                                                                                >  1:1047690/1‑140 (MQ=255)                                  gCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGa                                                                                                <  1:113645/140‑1 (MQ=255)                                       gTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGt                                                                                           <  1:508800/140‑1 (MQ=255)                                        tttttAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTc                                                                                          >  2:26117/1‑140 (MQ=255)                                         ttttAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCa                                                                                         <  1:28495/140‑1 (MQ=255)                                             aGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAAc                                                                                     >  2:324297/1‑140 (MQ=255)                                              gATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACa                                                                                    <  1:26117/140‑1 (MQ=255)                                                  gTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATa                                                                                       <  1:709149/133‑1 (MQ=255)                                                  gTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATa                                                                                       >  2:709149/1‑133 (MQ=255)                                                  gTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAaga                                                                                >  1:586312/1‑140 (MQ=255)                                                  gTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAaga                                                                                >  2:895477/1‑140 (MQ=255)                                                    aTCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAgg                                                                              >  1:432854/1‑140 (MQ=255)                                                             cTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGa                                                                     >  1:132927/1‑140 (MQ=255)                                                             cTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGa                                                                     >  2:655030/1‑140 (MQ=255)                                                             cTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGa                                                                     >  1:354475/1‑140 (MQ=255)                                                              tttGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGaa                                                                    <  1:895477/140‑1 (MQ=255)                                                                       atTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTAt                                                           >  1:188573/1‑140 (MQ=255)                                                                            ggCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACtt                                                             >  2:938022/1‑133 (MQ=255)                                                                            ggCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACtt                                                             <  1:938022/133‑1 (MQ=255)                                                                            ggCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGttt                                                      >  2:146232/1‑140 (MQ=255)                                                                            ggCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGttt                                                      >  1:1059823/1‑140 (MQ=255)                                                                            ggCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGttt                                                      >  2:435936/1‑140 (MQ=255)                                                                            ggCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGttt                                                      >  1:314428/1‑140 (MQ=255)                                                                                     ttgggttgggGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTAt                                             >  2:792122/1‑140 (MQ=255)                                                                                            gggggTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCAtt                                      >  1:745837/1‑140 (MQ=255)                                                                                             ggggTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTc                                     >  2:518890/1‑140 (MQ=255)                                                                                                 tACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATaa                                 >  1:989520/1‑140 (MQ=255)                                                                                                           aTTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCa                                        <  2:364546/123‑1 (MQ=255)                                                                                                           aTTGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCa                                        >  1:364546/1‑123 (MQ=255)                                                                                                            ttGCTTTTGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTgcg                      <  2:1059823/140‑1 (MQ=255)                                                                                                                ttttGCATGAACATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAgg                  >  2:743774/1‑140 (MQ=255)                                                                                                                         aaCATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCAt         >  2:592523/1‑140 (MQ=255)                                                                                                                          aCATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCAtt        >  2:713368/1‑140 (MQ=255)                                                                                                                           cATTGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCAttt       <  2:132927/140‑1 (MQ=255)                                                                                                                              tGCGCCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCATTTCAg    <  1:518890/140‑1 (MQ=255)                                                                                                                                cgcCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCATTTCAGGt  <  1:655030/140‑1 (MQ=255)                                                                                                                                cgcCTTTATACATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCATTTCAGGt  <  2:989520/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   ACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAATTGGGTTGGGGGTACAATTGTGATTGCTTTTGCATGAACATTGCGCCTTTATGCATAATGAGATAAAGGAATATCAAATAAAATAACGATAGGTCATAACAAAGAGGTTTTTATGAAAACACTTATCGTTTCAACTGTATTGGCATTCATAACATTTTCTGCGCAGGCTGCAGCATTTCAGGT  >  NC_000913/570697‑570963

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