Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,644,099 T→G A516A (GCA→GCC yfhM ← alpha2‑macroglobulin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,644,0990TG100.0% 46.4 / NA 14A516A (GCA→GCCyfhMalpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/6);  total (8/6)

AACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGCTATATCGCCGAATTCAAAACCAGGTAAGGCTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGA  >  NZ_CP009273/2643965‑2644231
                                                                                                                                      |                                                                                                                                    
aaCTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGCTATATCGCCGAATTCAAAACCAGGTAAGGCTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCAc                                                                                                                                  <  1:13687/139‑1 (MQ=255)
     aaCTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGCTATATCGCCGAATTCAAAACCAGGTAAGGCTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCACCATa                                                                                                                              >  1:109465/1‑138 (MQ=255)
                                                       cAAAACCAGGTAAGGCTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCg                                                                            <  1:227261/138‑1 (MQ=255)
                                                                               cTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAg                                                   >  2:181219/1‑139 (MQ=255)
                                                                                    cGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCa                                              >  2:167688/1‑139 (MQ=255)
                                                                                         tGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCAc                                                                                                                                  >  1:388288/1‑50 (MQ=255)
                                                                                         tGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCAc                                                                                                                                  <  2:388288/50‑1 (MQ=255)
                                                                                         tGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATgcgc                                         >  1:70974/1‑139 (MQ=255)
                                                                                           ggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGctc                                       <  2:109465/139‑1 (MQ=255)
                                                                                                                gCAAAGTATTACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTg                  >  1:36478/1‑139 (MQ=255)
                                                                                                                        ttACCATTAGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATccc          <  2:70974/139‑1 (MQ=255)
                                                                                                                                aGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATccc          <  1:341704/131‑1 (MQ=255)
                                                                                                                                aGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATccc          >  2:341704/1‑131 (MQ=255)
                                                                                                                                aGCAGGGGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGa  >  1:118907/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                    
AACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGCTATATCGCCGAATTCAAAACCAGGTAAGGCTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGA  >  NZ_CP009273/2643965‑2644231

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