Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,466,313	A→G	L138P (CTT→CCT)	bfr ←	bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,466,313	0	A	G	100.0%	278.2 / NA	75	L138P (CTT→CCT)	bfr	bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/36);  total (39/36)

AACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466188‑3466447                                                                                                                              |                                                                                                                                       aaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           <  2:576085/139‑1 (MQ=255)    cTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcca                                                                                                                        <  2:678149/139‑1 (MQ=255)           aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:830908/139‑1 (MQ=255)               ggCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                             <  1:1570407/139‑1 (MQ=255)                          ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTc                                                                                                                              <  1:1162270/111‑1 (MQ=255)                          ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTc                                                                                                                              >  2:1162270/1‑111 (MQ=255)                          ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                                                                                  >  1:196078/1‑139 (MQ=255)                                    gggAGGGTTCTCCCTCCCGACACGGCTCAACATTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCCGGTTCCGTTTCCAGCCAGTCGATATGACCTTCTTCATCACGCAAAAttt                                                                                        >  1:542281/1‑139 (MQ=255)                                     ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                    >  2:972625/1‑94 (MQ=255)                                     ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                    <  1:972625/94‑1 (MQ=255)                                     ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:273846/113‑1 (MQ=255)                                     ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 >  1:294617/1‑113 (MQ=255)                                     ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 >  1:273846/1‑113 (MQ=255)                                     ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:294617/113‑1 (MQ=255)                                            tctcCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     >  2:699552/1‑134 (MQ=255)                                            tctcCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     <  1:699552/134‑1 (MQ=255)                                                     cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                                                                           >  1:1640117/1‑119 (MQ=255)                                                     cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                                                                           <  2:1640117/119‑1 (MQ=255)                                                     cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg                                                                       >  2:922359/1‑139 (MQ=255)                                                      ggcacGGGTCAACCTGCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCTATCTTCTGAAGCAGATCAGTTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                      <  2:542281/137‑1 (MQ=255)                                                      gACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                      <  2:196078/139‑1 (MQ=255)                                                        cacGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg                                                                    >  1:1500554/1‑139 (MQ=255)                                                          cGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATc                                                                           >  2:803069/1‑130 (MQ=255)                                                          cGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATc                                                                           <  1:803069/130‑1 (MQ=255)                                                            gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg                                                                    <  1:1153974/135‑1 (MQ=255)                                                            gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg                                                                    >  2:1153974/1‑135 (MQ=255)                                                                 aCCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCa                                                           >  1:1359797/1‑139 (MQ=255)                                                                       ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  1:828698/73‑1 (MQ=255)                                                                       ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       >  2:828698/1‑73 (MQ=255)                                                                         cgcgGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCt                                                   >  2:54046/1‑139 (MQ=255)                                                                            ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              >  1:719974/1‑93 (MQ=255)                                                                            ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              <  2:719974/93‑1 (MQ=255)                                                                                 tgtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCg                                                                                                                    >  1:1586245/1‑66 (MQ=255)                                                                                 tgtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCg                                                                                                                    <  2:1586245/66‑1 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                     <  1:1513287/47‑1 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                     >  2:1513287/1‑47 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  2:115086/1‑69 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  1:115086/69‑1 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     >  2:1390771/1‑95 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     <  1:1390771/95‑1 (MQ=255)                                                                                   tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATaa                                         >  1:500466/1‑139 (MQ=255)                                                                                        gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  1:1353595/120‑1 (MQ=255)                                                                                        gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       >  2:1353595/1‑120 (MQ=255)                                                                                        gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAAt                                    >  2:1409947/1‑139 (MQ=255)                                                                                         cAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  1:980539/1‑87 (MQ=255)                                                                                         cAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       <  2:980539/87‑1 (MQ=255)                                                                                          aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                                                                                  >  1:1472497/1‑91 (MQ=255)                                                                                          aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatca                                                                                  <  2:1472497/91‑1 (MQ=255)                                                                                          aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTg                                  >  1:429823/1‑139 (MQ=255)                                                                                          aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTg                                  >  2:37195/1‑139 (MQ=255)                                                                                             tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcg                                                                          <  1:753318/96‑1 (MQ=255)                                                                                             tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcg                                                                          >  2:753318/1‑96 (MQ=255)                                                                                             tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCt                               <  1:37195/139‑1 (MQ=255)                                                                                              aaTTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTc                              >  1:547990/1‑139 (MQ=255)                                                                                                ttttGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCt                               >  1:1173899/1‑136 (MQ=255)                                                                                                ttttGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCt                               <  2:1173899/136‑1 (MQ=255)                                                                                                   tGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTa                                                  >  2:41579/1‑114 (MQ=255)                                                                                                   tGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTa                                                  <  1:41579/114‑1 (MQ=255)                                                                                                         cccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttct                   <  2:500466/139‑1 (MQ=255)                                                                                                          ccATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       >  1:890790/1‑38 (MQ=38)                                                                                                          ccATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  2:890790/38‑1 (MQ=38)                                                                                                          ccATCTTCTGAATCAGATCAGGTTCCGTTTccagcc                                                                                                                         <  2:710623/36‑1 (MQ=255)                                                                                                          ccATCTTCTGAATCAGATCAGGTTCCGTTTccagcc                                                                                                                         >  1:710623/1‑36 (MQ=255)                                                                                                             tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgc               <  2:547990/139‑1 (MQ=255)                                                                                                             tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgc               >  2:1478800/1‑139 (MQ=255)                                                                                                                      tCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       >  1:1338901/1‑90 (MQ=255)                                                                                                                      tCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  2:1338901/90‑1 (MQ=255)                                                                                                                         gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTGTTCATCACGCAAAAtt                                                                                         <  2:1487736/53‑1 (MQ=39)                                                                                                                         gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTGTTCATCACGCAAAAtt                                                                                         >  1:1487736/1‑53 (MQ=39)                                                                                                                         gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                <  2:395440/46‑1 (MQ=255)                                                                                                                         gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                >  1:395440/1‑46 (MQ=255)                                                                                                                         gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              <  1:520247/48‑1 (MQ=255)                                                                                                                         gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              >  2:520247/1‑48 (MQ=255)                                                                                                                          aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  >  2:706251/1‑139 (MQ=255)                                                                                                                          aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  <  2:1500554/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                       AACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466188‑3466447

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