Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,181,782	A→G	intergenic (‑57/+1)	ymfA ← / ← potD	DUF3592 family inner membrane protein/spermidine/putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,181,782	0	A	G	100.0%	166.8 / NA	45	intergenic (‑57/+1)	ymfA/potD	DUF3592 family inner membrane protein/spermidine/putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/22);  total (23/22)

AGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCC  >  NC_000913/1181664‑1181906                                                                                                                       |                                                                                                                             aGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCt                                                                                                          <  1:186183/139‑1 (MQ=255) aGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCAATAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCt                                                                                                          <  2:803360/139‑1 (MQ=255)        gTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGata                                                                                                   >  1:732745/1‑139 (MQ=255)           ccAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                <  2:892219/139‑1 (MQ=255)           ccAGACCGATACCGAGGAATATACGAAAGAAAACCTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                <  1:864676/139‑1 (MQ=255)            cAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACGAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATAc                                                                                               <  2:732745/139‑1 (MQ=255)              gACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCt                                                                                                          >  1:805540/1‑126 (MQ=255)              gACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCt                                                                                                          <  2:805540/126‑1 (MQ=255)                              tataCGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGctgctg                                                                             >  2:1030841/1‑139 (MQ=255)                                     aagaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTgcgc                                                                      >  2:1057487/1‑139 (MQ=255)                                      agaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAAc                                                                                                                         >  2:515272/1‑87 (MQ=255)                                      agaaaACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAAc                                                                                                                         <  1:515272/87‑1 (MQ=255)                                        aaaaCTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCaa                                                                   <  1:1030841/139‑1 (MQ=255)                                                              ttAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                >  1:811450/1‑88 (MQ=255)                                                              ttAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                <  2:811450/88‑1 (MQ=255)                                                              ttAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAAc                                                                  >  2:448433/1‑118 (MQ=255)                                                              ttAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAAc                                                                  <  1:448433/118‑1 (MQ=255)                                                              ttAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCAtt                                             >  1:617320/1‑139 (MQ=255)                                                                aGCCACATCCTTGCTAACTAAATAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                <  1:768362/86‑1 (MQ=255)                                                                aGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGataata                                                                                                >  2:768362/1‑86 (MQ=255)                                                                   cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGata                                                                                                   >  1:1157787/1‑80 (MQ=255)                                                                   cacaTCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGata                                                                                                   <  2:1157787/80‑1 (MQ=255)                                                                       tCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGt                                    >  2:800979/1‑139 (MQ=255)                                                                       tCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGt                                    >  2:1121240/1‑139 (MQ=255)                                                                           tGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCa                                >  2:976333/1‑139 (MQ=255)                                                                           tGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCa                                >  1:302564/1‑139 (MQ=255)                                                                            gCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAg                               >  1:849652/1‑139 (MQ=255)                                                                              taactaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCa                             <  1:1057487/139‑1 (MQ=255)                                                                                actaaAAAACGGGCGTTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATc                           <  2:849652/139‑1 (MQ=255)                                                                                actaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATA‑‑CTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCg                         >  1:1108297/1‑139 (MQ=255)                                                                                 ctaaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATcc                          <  1:469044/139‑1 (MQ=255)                                                                                  taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATaaa                                                                                     <  2:439354/79‑1 (MQ=255)                                                                                  taaAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATaaa                                                                                     >  1:439354/1‑79 (MQ=255)                                                                                    aaaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCggg                       <  2:617320/139‑1 (MQ=255)                                                                                     aaaaCGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGt                      <  2:302564/139‑1 (MQ=255)                                                                                      aaaCGGGCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTa                     <  1:800979/139‑1 (MQ=255)                                                                                          gggCGGTAATACCACCGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAgt                 >  2:576638/1‑139 (MQ=255)                                                                                                       accGCCCGCTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCTGGTAGAGTGTTTTATCGTTCg    >  2:694124/1‑139 (MQ=255)                                                                                                        ccgcccgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTc                                                                                                        >  1:1174533/1‑38 (MQ=255)                                                                                                        ccgcccgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTc                                                                                                        <  2:1174533/38‑1 (MQ=255)                                                                                                         cgcccgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGcc  >  2:415317/1‑139 (MQ=255)                                                                                                          gcccgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAA‑GCTGCTGGCTGCGCCAACg                                                                 >  1:877675/1‑74 (MQ=255)                                                                                                          gcccgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAA‑GCTGCTGGCTGCGCCAACg                                                                 <  2:877675/74‑1 (MQ=255)                                                                                                             cgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt              <  1:365381/123‑1 (MQ=255)                                                                                                             cgcTTGCTGAGTTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGtt              >  2:365381/1‑123 (MQ=255)                                                                                                                       |                                                                                                                             AGGATCAGTACCAGACCGATACCGAGGAATATACGAAAGAAAACTTTACTGTCCTGACTCATTAGCCACATCCTTGCTAACTAAAAAACGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCC  >  NC_000913/1181664‑1181906

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