Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,043,171	A→G	intergenic (‑23/+141)	yqfA ← / ← yqfB	hemolysin III family HylIII inner membrane protein/UPF0267 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,043,171	0	A	G	100.0%	171.2 / NA	50	intergenic (‑23/+141)	yqfA/yqfB	hemolysin III family HylIII inner membrane protein/UPF0267 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/24);  total (26/24)

CCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACACGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCAATTTAAAG  >  NC_000913/3043028‑3043317                                                                                                                                                |                                                                                                                                                   ccTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACAt                                                                                                                                             <  1:572380/151‑1 (MQ=255)          gTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACctc                                                                                                                                    >  2:703868/1‑151 (MQ=255)            aGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACctct                                                                                                                                   <  2:740153/150‑1 (MQ=255)              cAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTcc                                                                                                                                 >  1:318556/1‑150 (MQ=255)               aaCACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTcc                                                                                                                                 >  1:323964/1‑149 (MQ=255)                          aaCGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTg                                                                                                                       >  1:23841/1‑148 (MQ=255)                            cGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAAt                                                                                                                  <  1:736720/151‑1 (MQ=255)                             gATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAAt                                                                                                                  >  2:742646/1‑150 (MQ=255)                             gATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAAt                                                                                                                  >  1:93051/1‑150 (MQ=255)                             gATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAAt                                                                                                                  >  1:608376/1‑150 (MQ=255)                                   aaaCACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTcc                                                                                                            >  1:525724/1‑150 (MQ=255)                                       acCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTg                                                                                                        >  1:595375/1‑150 (MQ=255)                                          aaCCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAaca                                                                                                    <  1:348833/151‑1 (MQ=255)                                                    cGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCt                                                                                          >  2:305881/1‑151 (MQ=255)                                                     gTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTa                                                                                         <  1:706966/151‑1 (MQ=255)                                                        actgacGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCt                                                                                          <  1:205032/147‑1 (MQ=255)                                                        actgacGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTaaa                                                                                       <  2:547977/150‑1 (MQ=255)                                                          tgacGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAaaat                                                                                     >  1:670990/1‑150 (MQ=255)                                                                 gTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTg                                                                             <  2:124660/151‑1 (MQ=255)                                                                 gTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTg                                                                             <  1:42371/151‑1 (MQ=255)                                                                  ttGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTg                                                                             >  2:99901/1‑150 (MQ=255)                                                                   tGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGa                                                                            <  2:525776/150‑1 (MQ=255)                                                                    ggCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGAtt                                                                          <  2:743747/151‑1 (MQ=255)                                                                                 gCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAAtct                                                              <  2:662589/150‑1 (MQ=255)                                                                                   cAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCt                                                           >  1:352758/1‑151 (MQ=255)                                                                                     ttGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTa                                                          >  2:250025/2‑150 (MQ=255)                                                                                       gAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTaaaa                                                       <  2:730610/151‑1 (MQ=255)                                                                                         atatCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAAt                                                     <  1:99014/151‑1 (MQ=255)                                                                                          tatCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATc                                                    >  2:124736/1‑151 (MQ=255)                                                                                           atCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATc                                                    <  2:209959/150‑1 (MQ=255)                                                                                           atCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATc                                                    <  2:313784/150‑1 (MQ=255)                                                                                                    tAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTa                                           >  1:294710/1‑150 (MQ=255)                                                                                                         aGGGGCTTCTGAACCATAACTTACTCCGGATAAACGTACGCGCACATGTATACCTCACCAGCGTAACTGAGAATAGTTCCAGTGAACTCCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGct                                     <  1:595298/151‑2 (MQ=255)                                                                                                          ggggCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCg                                     >  2:613720/1‑150 (MQ=255)                                                                                                           gggCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGa                                    >  1:97761/1‑150 (MQ=255)                                                                                                             gCTTCTGAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAAt                                 >  2:372328/1‑151 (MQ=255)                                                                                                                   gAACCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAAc                           >  2:10400/1‑151 (MQ=255)                                                                                                                     aCCATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACt                          >  2:423034/1‑150 (MQ=255)                                                                                                                       cATAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAAc                       <  2:361779/151‑1 (MQ=255)                                                                                                                        aTAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAAc                       >  2:155815/1‑150 (MQ=255)                                                                                                                         tAACTTACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAAc                       >  1:400601/1‑149 (MQ=255)                                                                                                                              tACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCg                  <  2:624710/149‑1 (MQ=255)                                                                                                                              tACTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGt                 <  1:773714/150‑1 (MQ=255)                                                                                                                               aCTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTg                >  2:317712/1‑150 (MQ=255)                                                                                                                                cTCCGGAGAAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGa               <  2:616877/150‑1 (MQ=255)                                                                                                                                     gagaAACGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCaa         <  1:209500/151‑1 (MQ=255)                                                                                                                                        aaaCGTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCAAttt      <  2:208074/151‑1 (MQ=255)                                                                                                                                            gTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCAATTTaa    >  1:523193/1‑149 (MQ=255)                                                                                                                                            gTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCAATTTa     >  1:416943/1‑148 (MQ=255)                                                                                                                                            gTACGCGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCAATTTAAAg  <  1:232508/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACGCTGTTGGCTATTTCCTCTGCCAGTGAATATCCCTGCTTAATGAGGGGCTTCTGAACCATAACTTACTCCGGAGAAACGTACACGCACATGTATACCTCTCCAGCGTAACTGAGAATAGTTCCAGTGAACACCTGTTAGCTAAAATAAATTCTGATTTGAGACTAATCTCCTAAAAATCATGAAATTAAATGCGAAATTTCAACTAACAGGCGTGAGTTCAATTTAAAG  >  NC_000913/3043028‑3043317

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