Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,466,888	A→G	V35A (GTC→GCC)	bfd ←	bacterioferritin‑associated ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,466,888	0	A	G	100.0%	176.5 / NA	50	V35A (GTC→GCC)	bfd	bacterioferritin‑associated ferredoxin
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/29);  total (21/29)

CTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGAAAATAGTT  >  NC_000913/3466741‑3467034                                                                                                                                                    |                                                                                                                                                   cTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTgccggccg                                                                                                                                                 >  2:627010/1‑151 (MQ=255)  ttCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGccggccgg                                                                                                                                                >  2:229063/1‑151 (MQ=255)     cACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAAt                                                                                                                                             >  1:190851/1‑151 (MQ=255)     cACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAAt                                                                                                                                             <  1:673449/151‑1 (MQ=255)         ctTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATaa                                                                                                                                           <  1:60192/149‑1 (MQ=255)          tttGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAAtt                                                                                                                                        >  1:511073/1‑151 (MQ=255)            tGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAAttt                                                                                                                                       <  2:620760/150‑1 (MQ=255)                gagTAGAACGGATTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAtttttttt                                                                                                                                   >  2:344501/1‑150 (MQ=255)                    aGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTg                                                                                                                              <  2:70797/151‑1 (MQ=255)                             tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAAc                                                                                                                       <  2:526816/149‑1 (MQ=255)                             tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAAc                                                                                                                       <  2:271606/149‑1 (MQ=255)                             tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACg                                                                                                                      >  1:419447/1‑150 (MQ=255)                                 aGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                                                                                                  >  1:323292/1‑150 (MQ=255)                                      gTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGa                                                                                                            <  2:616781/151‑1 (MQ=255)                                        cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaa                                                                                                           >  2:67897/1‑150 (MQ=255)                                           aaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTGACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACt                                                                                                       <  2:374241/151‑1 (MQ=255)                                               cacacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                                                                                                   >  1:231652/1‑151 (MQ=255)                                                   cTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcac                                                                                                >  1:81735/1‑150 (MQ=255)                                                    tCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcac                                                                                                >  2:527865/1‑149 (MQ=255)                                                     ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcaca                                                                                               <  2:561977/149‑1 (MQ=255)                                                     ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAg                                                                                              <  2:69548/150‑1 (MQ=255)                                                     ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAg                                                                                              <  1:503463/150‑1 (MQ=255)                                                       gTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCt                                                                                           >  1:417274/1‑151 (MQ=255)                                                       gTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCt                                                                                           <  2:50795/151‑1 (MQ=255)                                                              cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAt                                                                                    <  1:648812/151‑1 (MQ=255)                                                              cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAt                                                                                    >  2:206363/1‑151 (MQ=255)                                                                         aaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGc                                                                         <  2:275483/151‑1 (MQ=255)                                                                             ccGGCAGCTGCATTAATTCATCCTCCATCCCTCCACGCGCGGCCCGGACACATTTCCCGCACTGATTGCCGGCCGGATTAAA‑TTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGAt                                                                     <  2:599218/150‑1 (MQ=255)                                                                             ccGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGAt                                                                     >  2:347370/1‑151 (MQ=255)                                                                              cGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATa                                                                    <  2:538665/151‑1 (MQ=255)                                                                                 cAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCa                                                                 <  1:89611/151‑1 (MQ=255)                                                                                       cATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACa                                                            >  2:3651/1‑150 (MQ=255)                                                                                               cATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACaaa                                                    >  2:494297/1‑150 (MQ=255)                                                                                                aTCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACg                                                  >  2:116912/1‑151 (MQ=255)                                                                                                  cctccATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGt                                                 >  1:389659/1‑150 (MQ=255)                                                                                                             tCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTc                                     >  1:442043/1‑151 (MQ=255)                                                                                                              cACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTca                                    <  2:355860/151‑2 (MQ=255)                                                                                                               aCGCGCGGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTccc                                   >  1:196859/1‑151 (MQ=255)                                                                                                                 gcgcgGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGt                                 >  1:210196/1‑151 (MQ=255)                                                                                                                  cgcgGCCCGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGta                                <  1:116914/151‑1 (MQ=255)                                                                                                                         cGGACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTttc                         <  2:412264/151‑1 (MQ=255)                                                                                                                          ggACACATTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTttct                        >  1:125543/1‑151 (MQ=255)                                                                                                                              acaTTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCaa                    <  2:602426/151‑1 (MQ=255)                                                                                                                              acaTTTACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCa                     >  2:478829/1‑150 (MQ=255)                                                                                                                                  ttACCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCaaaaa                 <  2:564568/150‑1 (MQ=255)                                                                                                                                    aCCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAaaataa              <  2:134543/151‑1 (MQ=255)                                                                                                                                    aCCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAaaata               <  2:201387/150‑1 (MQ=255)                                                                                                                                    aCCGCACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAaaata               <  1:484296/150‑1 (MQ=255)                                                                                                                                        cACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGa          <  2:106528/151‑1 (MQ=255)                                                                                                                                        cACTGATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGa          <  2:92591/151‑1 (MQ=255)                                                                                                                                            gATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGaaaa       <  2:47414/150‑1 (MQ=255)                                                                                                                                            gATTGCCGGCCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGaaa        <  1:676279/149‑1 (MQ=255)                                                                                                                                                gccggccgGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGAAAATAGtt  <  1:671771/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACATAGAGACTCCCGTATACTTTCTTCAAAAATAAATGAAAATAGTT  >  NC_000913/3466741‑3467034

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