Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 177,662 C→T *267* (TAG→TAA btuF ← vitamin B12 ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,6620CT100.0% 105.0 / NA 34*267* (TAG→TAAbtuFvitamin B12 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/19);  total (15/19)

GAACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTA  >  NC_000913/177528‑177784
                                                                                                                                      |                                                                                                                          
gAACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                        >  1:350142/1‑139 (MQ=255)
     ccAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTAcc                                                                                                                   >  2:545625/1‑139 (MQ=255)
                cATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcgc                                                                                                        >  2:47482/1‑139 (MQ=255)
                cATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAgcgc                                                                                                        >  2:1067887/1‑139 (MQ=255)
                    cGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTa                                                                                                    <  2:808810/139‑1 (MQ=255)
                    cGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTa                                                                                                    <  1:206607/139‑1 (MQ=255)
                      cAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACa                                                                                                  <  2:350142/139‑1 (MQ=255)
                       aaTTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACag                                                                                                 <  2:593122/139‑1 (MQ=255)
                            tCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCtg                                                                                            <  1:1067887/139‑1 (MQ=255)
                                     cAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                   <  1:545625/139‑1 (MQ=255)
                                             cGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGataa                                                                           >  1:43739/1‑139 (MQ=255)
                                                 gaTGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACg                                                                       <  2:43739/139‑1 (MQ=255)
                                                                        cTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGt                                                <  1:47482/139‑1 (MQ=255)
                                                                              ccAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCg                                          >  2:1071501/1‑139 (MQ=255)
                                                                                   cAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                        <  2:106450/56‑1 (MQ=255)
                                                                                   cAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                        >  1:106450/1‑56 (MQ=255)
                                                                                     aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                        <  2:225882/54‑1 (MQ=255)
                                                                                     aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAAt                                                                                                                        >  1:225882/1‑54 (MQ=255)
                                                                                     aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                  <  2:112235/92‑1 (MQ=255)
                                                                                     aTAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCg                                                                                  >  1:112235/1‑92 (MQ=255)
                                                                                      tAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                              <  1:581927/63‑1 (MQ=255)
                                                                                      tAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGa                                                                                                              >  2:581927/1‑63 (MQ=255)
                                                                                                     cAATCCCCACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTg                   <  1:503650/139‑1 (MQ=255)
                                                                                                         ccccACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                               >  1:487669/1‑59 (MQ=255)
                                                                                                         ccccACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACagag                                                                                               <  2:487669/59‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt            <  1:775471/139‑1 (MQ=255)
                                                                                                            cACTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt            <  2:819785/139‑1 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCa                                                                                     <  2:946237/65‑1 (MQ=255)
                                                                                                             aCTGAGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCa                                                                                     >  1:946237/1‑65 (MQ=255)
                                                                                                                 aGCACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAccc       >  1:927984/1‑139 (MQ=255)
                                                                                                                  gCACCTGAAAAGGCCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGa                                                                              >  2:498114/1‑67 (MQ=255)
                                                                                                                  gCACCTGAAAAGGCCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGa                                                                              <  1:498114/67‑1 (MQ=255)
                                                                                                                   cACCTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCa     >  2:523671/1‑139 (MQ=255)
                                                                                                                      cTGAAAAGGTCAGGCGTTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTa  <  1:1071501/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                          
GAACACCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTA  >  NC_000913/177528‑177784

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