Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,646,189 G→C 100% A155G (GCA→GGA)  USA300HOU_1531 ← M24B family Xaa‑Pro dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,646,1890GC100.0% 78.1 / NA 25A155G (GCA→GGA) USA300HOU_1531M24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (13/12);  total (13/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATA  >  CP000730/1646062‑1646323
                                                                                                                               |                                                                                                                                      
tGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTCCTTTTACACCAGt                                                                                                                           >  2:258339/1‑141 (MQ=255)
    acacCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaa                                                                                                                       >  1:139684/1‑141 (MQ=255)
       ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                    >  1:201594/1‑141 (MQ=255)
            tgGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTc                                                                                                               <  1:138783/141‑1 (MQ=255)
            gggTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTc                                                                                                               <  2:220104/140‑1 (MQ=255)
                  tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                         <  1:122815/141‑1 (MQ=255)
                    tctctctATGACCAGAGTCCACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGCTAAAATATATTCATATGttt                                                                                                       <  2:197221/138‑1 (MQ=255)
                            tGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcacca                                                                                               >  1:18994/1‑141 (MQ=255)
                                       tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATCTATTCATATGTTTCATCAACAATATTAtctgc                                                                                    >  2:100601/1‑141 (MQ=255)
                                        aCAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgct                                                                                   <  2:239799/141‑1 (MQ=255)
                                           aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                                                                <  2:37615/141‑1 (MQ=255)
                                                  cGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                                         >  1:176379/1‑141 (MQ=255)
                                                             ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg                                                              >  2:34235/1‑141 (MQ=255)
                                                              cATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGt                                                             <  2:44128/141‑1 (MQ=255)
                                                                tCTGCTCCTAATTATAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCa                                                           <  2:240719/141‑1 (MQ=255)
                                                                        tAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCttt                                                   <  1:216753/141‑1 (MQ=255)
                                                                             cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGt                                              >  1:281037/1‑141 (MQ=255)
                                                                                gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtctc                                           >  1:319105/1‑141 (MQ=255)
                                                                                 cATTTTGCTTTCTAATATTGCCTTTAATTCTTTCTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                          <  1:74433/141‑1 (MQ=255)
                                                                                                        ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCTTCAACAATATTAGCTGTTTTTTGTTTCCAATTAAGGGCGGGTTCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGaaa                   >  2:235046/1‑141 (MQ=255)
                                                                                                        ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAATCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACGGTATTAGaaa                   >  1:26880/1‑141 (MQ=255)
                                                                                                            ttCTCTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGCTTCATCAACAATGTTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTCATCTACTGTATTAGAAATGCt               <  2:228691/141‑1 (MQ=255)
                                                                                                            tgcTTTTTCAGTCATGCCGCCTTTTACAATCTTTAAAATTTATTCATATGTTTCATCAATAATATTAGCTGCTTTTTGACTTAAAGCTATTTCTTCAGCATCTTTGACGTCTCTAATTTTATCTACCGTATTAGAAATGCt               <  2:253443/139‑1 (MQ=255)
                                                                                                                  ttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTaa         >  2:106313/1‑141 (MQ=255)
                                                                                                                         aTGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGAtata  >  2:44462/1‑141 (MQ=255)
                                                                                                                               |                                                                                                                                      
TGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATA  >  CP000730/1646062‑1646323

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