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

ACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466189‑3466447                                                                                                                             |                                                                                                                                       aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagc                                                                                                                          <  1:357863/139‑1 (MQ=255)    tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  2:68557/139‑1 (MQ=255)                            attaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           <  1:634399/111‑1 (MQ=255)                            attaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           >  2:634399/1‑111 (MQ=255)                               aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaa                                                                                            >  2:221337/1‑139 (MQ=255)                                aCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                                                                           >  1:413201/1‑139 (MQ=255)                                 ccGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  2:335636/1‑118 (MQ=255)                                 ccGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  2:34517/1‑118 (MQ=255)                                 ccGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  1:335636/118‑1 (MQ=255)                                 ccGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  1:34517/118‑1 (MQ=255)                                    ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                    >  2:284749/1‑94 (MQ=255)                                    ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                    <  1:284749/94‑1 (MQ=255)                                    ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  2:878728/1‑139 (MQ=255)                                            ctccctccCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           >  2:170674/1‑95 (MQ=255)                                            ctccctccCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           <  1:170674/95‑1 (MQ=255)                                                  cccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgc                                                                         >  1:831006/1‑139 (MQ=255)                                                   ccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgcg                                                                        >  2:203974/1‑139 (MQ=255)                                                                 ccTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     >  1:504871/1‑112 (MQ=255)                                                                 ccTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     <  2:504871/112‑1 (MQ=255)                                                                         gcgGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTa                                                  <  2:624025/139‑1 (MQ=255)                                                                           ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGt                                                                                                                                      >  2:986838/1‑53 (MQ=255)                                                                           ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGAc                                                                  <  1:361122/121‑1 (MQ=255)                                                                           ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGAc                                                                  >  2:361122/1‑121 (MQ=255)                                                                               cTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCa                                            <  2:988019/139‑1 (MQ=255)                                                                                tgtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCAt                                           <  1:203974/139‑1 (MQ=255)                                                                                  tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           <  2:722869/57‑1 (MQ=255)                                                                                  tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  2:1014701/69‑1 (MQ=255)                                                                                  tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  1:1014701/1‑69 (MQ=255)                                                                                  tgCTTGCACATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                           >  1:722869/1‑57 (MQ=255)                                                                                       gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAAt                                    >  1:940213/1‑139 (MQ=255)                                                                                       gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACACTATCGGCATAACCAAt                                    >  1:880554/1‑139 (MQ=255)                                                                                         aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  2:461897/62‑1 (MQ=255)                                                                                         aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  1:461897/1‑62 (MQ=255)                                                                                          gATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  1:650387/61‑1 (MQ=255)                                                                                          gATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  2:650387/1‑61 (MQ=255)                                                                                            tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     <  2:1105870/85‑1 (MQ=255)                                                                                            tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     >  1:1105870/1‑85 (MQ=255)                                                                                                  tGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCa                         >  2:1015785/1‑139 (MQ=255)                                                                                                            tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               >  1:107148/1‑43 (MQ=255)                                                                                                            tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                               <  2:107148/43‑1 (MQ=255)                                                                                                             cttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg              <  1:221337/139‑1 (MQ=255)                                                                                                                 tGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                >  2:632988/1‑53 (MQ=255)                                                                                                                 tGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                <  1:632988/53‑1 (MQ=255)                                                                                                                   aaTCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       >  1:957877/1‑92 (MQ=255)                                                                                                                   aaTCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  2:957877/92‑1 (MQ=255)                                                                                                                        gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTc   >  1:48371/1‑139 (MQ=255)                                                                                                                        gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTc   >  2:493365/1‑139 (MQ=255)                                                                                                                         aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  >  1:1099912/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                       ACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466189‑3466447

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