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%	195.3 / NA	55	L138P (CTT→CCT)	bfr	bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/23);  total (32/23)

CTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTC  >  NC_000913/3466191‑3466446                                                                                                                           |                                                                                                                                      cTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcca                                                                                                                       >  1:1148211/1‑139 (MQ=255)  tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTTCATATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                      <  1:321305/139‑1 (MQ=255)     ggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                          <  1:475373/132‑1 (MQ=255)     ggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                          >  2:475373/1‑132 (MQ=255)        aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                <  1:104277/139‑1 (MQ=255)           tGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCg                                                                                                                   >  1:256164/1‑133 (MQ=255)           tGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCg                                                                                                                   <  2:256164/133‑1 (MQ=255)           tGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATgg                                                                                                             >  2:1773029/1‑139 (MQ=255)            ggCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                            <  1:1773029/139‑1 (MQ=255)                tGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCcttc                                                                                                        >  2:131997/1‑139 (MQ=255)                       ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                                                                                 <  2:1078135/139‑1 (MQ=255)                                cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAtt                                                                                        >  1:1120137/1‑139 (MQ=255)                                  ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcca                                                                                                                       >  2:762878/1‑106 (MQ=255)                                  ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcca                                                                                                                       <  1:762878/106‑1 (MQ=255)                                  ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                      <  2:1148211/139‑1 (MQ=255)                                  ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                      >  1:835852/1‑139 (MQ=255)                                    aGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                    <  2:559439/139‑1 (MQ=255)                                                  cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg                                                                      >  1:1711347/1‑139 (MQ=255)                                                     cacGGCTCAACCATCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg                                                                   <  2:1670608/139‑1 (MQ=255)                                                         gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg                                                                      >  1:476356/1‑132 (MQ=255)                                                         gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg                                                                      <  2:476356/132‑1 (MQ=255)                                                           tCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                   <  1:1610571/69‑1 (MQ=255)                                                           tCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                   >  2:1610571/1‑69 (MQ=255)                                                                    ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCTATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                >  2:1027302/1‑79 (MQ=255)                                                                    ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCTATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                <  1:1027302/79‑1 (MQ=255)                                                                    ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg                                                                   >  1:613972/1‑124 (MQ=255)                                                                    ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTg                                                                   <  2:613972/124‑1 (MQ=255)                                                                    ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACg                                                    >  2:1401799/1‑139 (MQ=255)                                                                         ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATc                                               >  2:1395118/1‑139 (MQ=255)                                                                                tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                                                                                <  2:1752272/99‑1 (MQ=255)                                                                                tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                                                                                >  1:1752272/1‑99 (MQ=255)                                                                                      cAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                             <  2:885231/80‑1 (MQ=255)                                                                                      cAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                             >  1:885231/1‑80 (MQ=255)                                                                                          tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTcc                                                                                                                            >  1:1247432/1‑45 (MQ=255)                                                                                          tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTcc                                                                                                                            <  2:1247432/45‑1 (MQ=255)                                                                                               ttGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                      <  2:1293521/110‑1 (MQ=255)                                                                                               ttGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                      >  1:1293521/1‑110 (MQ=255)                                                                                                tGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCa                        >  1:1309621/1‑139 (MQ=255)                                                                                                 gCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                            <  2:1470477/54‑1 (MQ=255)                                                                                                 gCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                            >  1:288420/1‑54 (MQ=255)                                                                                                 gCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                            >  1:1470477/1‑54 (MQ=255)                                                                                                 gCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                            <  2:288420/54‑1 (MQ=255)                                                                                                   aGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcc                                                                                                                        <  2:1458589/40‑1 (MQ=38)                                                                                                   aGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcc                                                                                                                        >  1:1458589/1‑40 (MQ=38)                                                                                                     gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                >  2:1564624/1‑46 (MQ=255)                                                                                                     gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                <  1:1564624/46‑1 (MQ=255)                                                                                                     gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                   >  1:1035764/1‑139 (MQ=255)                                                                                                     gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                   >  1:1700085/1‑139 (MQ=255)                                                                                                      cccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttct                  >  2:481901/1‑139 (MQ=255)                                                                                                        cATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACACAAAttcttc                >  2:1826398/1‑139 (MQ=255)                                                                                                           cttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg             <  2:1035764/139‑1 (MQ=255)                                                                                                                 aaTCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATcc       >  2:1314167/1‑139 (MQ=255)                                                                                                                 aaTCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATcc       >  2:1313171/1‑139 (MQ=255)                                                                                                                      gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCCTCTTTATCACGCAAAATTTCTATCATCATATCGCCGCTGACGTAATCATGAACGCTATCGGCATTACCAATTGTCTCACGCAAATTCTTCGCGACATCCAGCTc  >  1:877995/1‑139 (MQ=255)                                                                                                                        tCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCAt                                                         >  1:942845/1‑82 (MQ=255)                                                                                                                           |                                                                                                                                      CTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTC  >  NC_000913/3466191‑3466446

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