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% 122.7 / NA 35L138P (CTT→CCT) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/20);  total (15/20)

CAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTC  >  NC_000913/3466178‑3466446
                                                                                                                                       |                                                                                                                                     
cagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                    >  1:134172/1‑139 (MQ=255)
 agcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                   <  1:531669/139‑1 (MQ=255)
             cTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcca                                                                                                                       <  2:493207/139‑1 (MQ=255)
              tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                      <  1:137732/139‑1 (MQ=255)
                  gTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGa                                                                                                                  <  2:134172/139‑1 (MQ=255)
                    aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                <  2:728578/139‑1 (MQ=255)
                      cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCAGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATg                                                                                                              <  1:143298/139‑1 (MQ=255)
                        ggCCTGGATCTCTTATTAACCGGGCGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGc                                                                                                            <  2:449116/139‑1 (MQ=255)
                                   ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                                                                                 >  1:727821/1‑139 (MQ=255)
                                   ctTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGAAAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTtcatca                                                                                                 >  2:609313/1‑139 (MQ=255)
                                       ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                             <  1:409232/139‑1 (MQ=255)
                                           ccGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAt                                                                                         >  1:559270/1‑139 (MQ=255)
                                              ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAACTTTc                                                                                      >  1:143876/1‑139 (MQ=255)
                                              ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                      >  1:686667/1‑139 (MQ=255)
                                                               gACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                     <  1:624261/139‑1 (MQ=255)
                                                                  tcGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTTAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa                                                                  <  1:704147/138‑1 (MQ=255)
                                                                    ggCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg                                                                <  1:620540/139‑1 (MQ=255)
                                                                     gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt                                                               >  1:594819/1‑139 (MQ=255)
                                                                      cTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTa                                                              >  2:396368/1‑139 (MQ=255)
                                                                              tcttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                      <  2:686667/139‑1 (MQ=255)
                                                                              tcttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                      <  2:727821/139‑1 (MQ=255)
                                                                                    cgGATCTGTGCTTGCAGAAAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTAt                                                <  1:609313/139‑1 (MQ=255)
                                                                                     ggATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATc                                               >  2:168329/1‑139 (MQ=255)
                                                                                                 gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATg                                                        <  1:532831/118‑1 (MQ=255)
                                                                                                 gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATg                                                        >  2:532831/1‑118 (MQ=255)
                                                                                                   aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCTATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTg                                 >  1:249719/1‑139 (MQ=255)
                                                                                                               aGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGACCTTCTTAATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAt                     >  2:60444/1‑139 (MQ=255)
                                                                                                                  cccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttct                  <  2:559270/139‑1 (MQ=255)
                                                                                                                   ccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                >  2:225630/1‑44 (MQ=255)
                                                                                                                   ccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                <  1:225630/44‑1 (MQ=255)
                                                                                                                       cttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg             <  1:396368/139‑1 (MQ=255)
                                                                                                                       cttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg             <  2:594819/139‑1 (MQ=255)
                                                                                                                                  gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTc  >  2:784089/1‑139 (MQ=255)
                                                                                                                                   aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg                                                                <  1:774632/76‑1 (MQ=255)
                                                                                                                                   aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg                                                                >  2:774632/1‑76 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTC  >  NC_000913/3466178‑3466446

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