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_0009133,466,3130AG100.0% 207.0 / NA 57L138P (CTT→CCT) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/24);  total (33/24)

GGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466193‑3466447
                                                                                                                        |                                                                                                                                      
ggggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGt                                                                                                                      <  2:1302017/139‑1 (MQ=255)
 gggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTc                                                                                                                     >  1:170808/1‑139 (MQ=255)
 gggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTc                                                                                                                     >  2:642883/1‑139 (MQ=255)
      tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGGATCAGATCAGGTTCCGTTTccagccag                                                                                                                       >  1:843830/1‑132 (MQ=255)
      tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCGTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  2:843830/132‑1 (MQ=255)
           ccTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCct                                                                                                           >  1:1091338/1‑139 (MQ=255)
                    ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTC‑CGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                     >  2:1219741/1‑103 (MQ=255)
                    ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTC‑CGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                     <  1:1219741/103‑1 (MQ=255)
                         taACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:897641/119‑1 (MQ=255)
                         taACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 >  1:897641/1‑119 (MQ=255)
                             cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTGGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAtt                                                                                         <  2:331612/139‑1 (MQ=255)
                               ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              <  1:466427/132‑1 (MQ=255)
                               ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              >  2:466427/1‑132 (MQ=255)
                               ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  1:943777/1‑139 (MQ=255)
                               ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  1:281535/1‑139 (MQ=255)
                                                gACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                      <  1:473917/139‑1 (MQ=255)
                                                    cGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGAc                                                                  >  2:1271268/1‑139 (MQ=255)
                                                     ggCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg                                                                 <  1:512667/139‑1 (MQ=255)
                                                     ggCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg                                                                 <  2:943777/139‑1 (MQ=255)
                                                      gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg                                                                       <  2:268506/132‑1 (MQ=255)
                                                      gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg                                                                       >  1:268506/1‑132 (MQ=255)
                                                      gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt                                                                >  2:432709/1‑139 (MQ=255)
                                                      gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt                                                                >  1:721458/1‑139 (MQ=255)
                                                       cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATACTTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCt                                                                                      <  2:776102/116‑1 (MQ=255)
                                                       cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  2:703153/83‑1 (MQ=255)
                                                       cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       >  1:703153/1‑83 (MQ=255)
                                                       cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                <  2:174238/106‑1 (MQ=255)
                                                       cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                >  1:174238/1‑106 (MQ=255)
                                                       cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCt                                                                                      >  1:776102/1‑116 (MQ=255)
                                                        tCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                <  1:951596/105‑1 (MQ=255)
                                                        tCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                >  2:951596/1‑105 (MQ=255)
                                                          aaCCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc                                                            >  1:334684/1‑139 (MQ=255)
                                                               tcttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGt                                                                                                                      <  1:1152130/76‑1 (MQ=255)
                                                               tcttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGt                                                                                                                      >  2:1152130/1‑76 (MQ=255)
                                                               tcttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  2:1188250/139‑1 (MQ=255)
                                                               tcttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  1:642883/139‑1 (MQ=255)
                                                                cttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 >  2:1283979/1‑80 (MQ=255)
                                                                cttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  1:1283979/80‑1 (MQ=255)
                                                                    gcgGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTa                                                  <  1:1271268/139‑1 (MQ=255)
                                                                        aTCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCgg                                              <  1:567663/139‑1 (MQ=255)
                                                                             tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATaa                                         >  2:208443/1‑139 (MQ=255)
                                                                                  gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGATGTAATCATGAACGCTATCGGCCTAACCAAt                                    >  2:567067/1‑139 (MQ=255)
                                                                                   cAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAAtt                                   >  1:845650/1‑139 (MQ=255)
                                                                                    aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTg                                  >  1:591883/1‑139 (MQ=255)
                                                                                       tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              <  1:566564/76‑1 (MQ=255)
                                                                                       tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              >  2:566564/1‑76 (MQ=255)
                                                                                               cAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCaaa                       >  2:44663/1‑139 (MQ=255)
                                                                                                  gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                    >  2:388504/1‑139 (MQ=255)
                                                                                                  gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                    >  1:447191/1‑139 (MQ=255)
                                                                                                  gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                    >  2:927459/1‑139 (MQ=255)
                                                                                                  gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                    >  2:1134409/1‑139 (MQ=255)
                                                                                                       tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 >  1:1024907/1‑41 (MQ=255)
                                                                                                       tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:1024907/41‑1 (MQ=255)
                                                                                                       tcttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGCTATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGACAGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgc               >  2:145615/1‑139 (MQ=255)
                                                                                                                   gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTc   >  1:1046758/1‑139 (MQ=255)
                                                                                                                    aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  <  2:334684/139‑1 (MQ=255)
                                                                                                                    aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCATCACGCAAATTCTTCGCGCCATCCAGCTCa  <  1:208443/139‑1 (MQ=255)
                                                                                                                        |                                                                                                                                      
GGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466193‑3466447

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