Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,314,744 T→C T161A (ACT→GCT)  yciE ← putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,314,7440TC100.0% 147.2 / NA 39T161A (ACT→GCT) yciEputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/16);  total (23/16)

AACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTG  >  NC_000913/1314597‑1314887
                                                                                                                                                   |                                                                                                                                               
aaCACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTc                                                                                                                                              >  1:663347/1‑151 (MQ=255)
     gATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGc                                                                                                                                         >  2:47047/1‑151 (MQ=255)
         gCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAAt                                                                                                                                     <  1:627782/151‑1 (MQ=255)
             aCGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTaa                                                                                                                                  <  1:5280/150‑1 (MQ=255)
                     aTCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTtctc                                                                                                                         <  1:194933/151‑1 (MQ=255)
                         gTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAgtt                                                                                                                     >  2:276024/1‑151 (MQ=255)
                            ggTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAgttgtt                                                                                                                  >  2:626995/1‑151 (MQ=255)
                                          cagcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTg                                                                                                    <  2:295207/151‑1 (MQ=255)
                                               ggATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATc                                                                                                >  2:561471/1‑150 (MQ=255)
                                                 atatCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGc                                                                                             <  1:275374/151‑1 (MQ=255)
                                                  tatCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGcc                                                                                            >  1:140536/1‑151 (MQ=255)
                                                   atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGcc                                                                                            >  1:307077/1‑150 (MQ=255)
                                                     cGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAg                                                                                          <  1:3321/150‑1 (MQ=255)
                                                      gTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTc                                                                                        >  2:99151/1‑151 (MQ=255)
                                                      gTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTc                                                                                        <  1:623376/151‑1 (MQ=255)
                                                          ttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGcc                                                                                    >  1:43828/1‑151 (MQ=255)
                                                          ttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGc                                                                                     <  2:486606/150‑1 (MQ=255)
                                                            aaTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCa                                                                                   <  1:171770/150‑1 (MQ=255)
                                                                ccGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGTCCATTTg                                                                               <  1:488680/150‑1 (MQ=255)
                                                                 cGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCt                                                                             >  1:41856/1‑151 (MQ=255)
                                                                              ccccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaa                                                                <  2:113479/151‑1 (MQ=255)
                                                                               cccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaaa                                                               <  1:561466/151‑1 (MQ=255)
                                                                                     atatACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCt                                                         <  2:43827/151‑1 (MQ=255)
                                                                                      tataCGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCtt                                                        >  2:29819/1‑151 (MQ=255)
                                                                                          cGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGAt                                                    >  1:450410/1‑151 (MQ=255)
                                                                                               tCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTg                                               >  2:358742/1‑151 (MQ=255)
                                                                                                     aGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTg                                         >  1:407786/1‑151 (MQ=255)
                                                                                                       gACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGaa                                       <  2:550247/151‑1 (MQ=255)
                                                                                                           tGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTg                                   >  1:547535/1‑151 (MQ=255)
                                                                                                              tataAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTAt                                >  2:228976/1‑151 (MQ=255)
                                                                                                              tataAAGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTAt                                >  2:94606/1‑151 (MQ=255)
                                                                                                                  aaGGGTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCAcc                            >  1:590026/1‑151 (MQ=255)
                                                                                                                      gTATTATTTCTTCGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGc                         >  2:284814/1‑150 (MQ=255)
                                                                                                                                 tcGCTTCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTtgctgc             >  2:232425/1‑151 (MQ=255)
                                                                                                                                      tCTACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTtgctgctgc          <  1:94604/149‑1 (MQ=255)
                                                                                                                                        tACGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAAt      >  2:610370/1‑151 (MQ=255)
                                                                                                                                         aCGCCATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAAt      <  1:590177/150‑1 (MQ=255)
                                                                                                                                             cATCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTg  >  2:412179/1‑150 (MQ=255)
                                                                                                                                              aTCAGCTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTg  >  2:249404/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTG  >  NC_000913/1314597‑1314887

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