Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,011,856	A→G	T261A (ACG→GCG)	bioP →	biotin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,011,856	0	A	G	100.0%	183.2 / NA	51	T261A (ACG→GCG)	bioP	biotin transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/29);  total (22/29)

GGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAA  >  NC_000913/4011727‑4011987                                                                                                                                  |                                                                                                                                    ggCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTAc                                                                                                                           >  2:482543/1‑140 (MQ=255) ggCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTAc                                                                                                                           >  2:383924/1‑140 (MQ=255)    tACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAgg                                                                                                                        <  1:157288/140‑1 (MQ=255)    tACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAgg                                                                                                                        <  2:115437/140‑1 (MQ=255)       tttATGTGGAACTACGGTGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGgcgc                                                                                                                     >  1:112139/1‑140 (MQ=255)       tttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGgcgc                                                                                                                     <  2:670954/140‑1 (MQ=255)       tttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGgcgc                                                                                                                     >  1:55555/1‑140 (MQ=255)            gtgGAACTACGGCGCGATGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTgg                                                                                                                <  2:469043/140‑1 (MQ=255)            gtgGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTgg                                                                                                                <  1:383924/140‑1 (MQ=255)              ggAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTg                                                                                                              <  1:623688/140‑1 (MQ=255)                aaCTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGAt                                                                                                            >  1:373090/1‑140 (MQ=255)                      ggCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGc                                                                                                      <  1:94120/140‑1 (MQ=255)                      ggCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGc                                                                                                      <  2:663170/140‑1 (MQ=255)                      ggCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGc                                                                                                      <  1:717287/140‑1 (MQ=255)                         gcgACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTc                                                                                                   <  2:373090/140‑1 (MQ=255)                         gcgACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTc                                                                                                   <  2:758881/140‑1 (MQ=255)                             cGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACtg                                                                                               >  2:201643/1‑140 (MQ=255)                             cGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACtg                                                                                               <  1:201643/140‑1 (MQ=255)                                     gACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCa                                                                                       <  2:285882/140‑1 (MQ=255)                                          cGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTa                                                                                  >  1:845828/1‑140 (MQ=255)                                             aaCGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGt                                                                               >  2:58218/1‑140 (MQ=255)                                                  tGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTc                                                                          <  2:112139/140‑1 (MQ=255)                                                          aTGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCcgcg                                                                  >  2:401537/1‑140 (MQ=255)                                                                aataTGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTc                                                            >  1:212570/1‑140 (MQ=255)                                                                   atGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTcaca                                                         <  1:58218/140‑1 (MQ=255)                                                                   atGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTcaca                                                         >  1:357716/1‑140 (MQ=255)                                                                   atGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTcaca                                                         <  2:104036/140‑1 (MQ=255)                                                                   atGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTcaca                                                         >  1:457375/1‑140 (MQ=255)                                                                    tGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaa                                                        >  1:602453/1‑140 (MQ=255)                                                                    tGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaa                                                        >  2:430384/1‑140 (MQ=255)                                                                    tGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaa                                                        >  2:308596/1‑140 (MQ=255)                                                                     gCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaaa                                                       <  2:46235/140‑1 (MQ=255)                                                                     gCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaaa                                                       <  1:322129/140‑1 (MQ=255)                                                                     gCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaaa                                                       <  2:602453/140‑1 (MQ=255)                                                                     gCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaaa                                                       <  1:430384/140‑1 (MQ=255)                                                                            ccGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAga                                                >  2:459347/1‑140 (MQ=255)                                                                                  gggCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCg                                          <  1:308596/140‑1 (MQ=255)                                                                                   ggCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGc                                         >  1:445840/1‑140 (MQ=255)                                                                                                  tGGCTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTgcgcg                          <  1:750460/140‑1 (MQ=255)                                                                                                     cTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTg                              >  1:47315/1‑133 (MQ=255)                                                                                                     cTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTg                              <  2:47315/133‑1 (MQ=255)                                                                                                     cTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTg                       >  2:347417/1‑140 (MQ=255)                                                                                                     cTATCTGGCACCAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTAATTGCGCGCTg                       >  1:478303/1‑140 (MQ=255)                                                                                                               ccAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATaaa             <  1:300001/140‑1 (MQ=255)                                                                                                                cAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAAc            <  2:845828/140‑1 (MQ=255)                                                                                                                cAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAAc            <  2:457375/140‑1 (MQ=255)                                                                                                                cAACAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAAc            <  1:398905/140‑1 (MQ=255)                                                                                                                 aaCAGCCGCACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACg           <  2:55555/140‑1 (MQ=255)                                                                                                                         cACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTa   >  1:321566/1‑140 (MQ=255)                                                                                                                         cACTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTa   >  2:164706/1‑140 (MQ=255)                                                                                                                          aCTGGCCAGCGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTaa  <  1:459347/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                    GGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAA  >  NC_000913/4011727‑4011987

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: ‑