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.9 / NA 27A155G (GCA→GGA) USA300HOU_1531M24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (15/12);  total (15/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATAT  >  CP000730/1646059‑1646322
                                                                                                                                  |                                                                                                                                     
aCTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                                                                             >  1:281759/1‑141 (MQ=255)
        cacCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                     <  1:258516/141‑1 (MQ=255)
              gtgGTAATGCACCTCTATGACGAGATGCTACAATCGAATCGAAGGATGGTCTAGCTGATCCTCATGCTAGAGTTTTGCTTGCTAATATTGCCTTTAATTCTTTTTCAGTGGTGCCTCCTTTTACAACAGTTgaagtctctt                                                                                                               >  2:139711/1‑134 (MQ=255)
                ggTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCa                                                                                                             <  1:88829/141‑1 (MQ=255)
                     tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                        <  1:275389/141‑1 (MQ=255)
                                 aCCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAat                                                                                            >  1:49098/1‑141 (MQ=255)
                                    ctaTGATACAGTCGTCGCGAATGATGGTCCATCTGCTCCTGATTCTAGTACTTTGCTTTCGAATATTGACTTTAATGCTCTTTCAGTCACGCCTCCTTTTACAACAATTAAAATATATTCATATGTTTCATCAACAATAtt                                                                                         <  2:98938/139‑1 (MQ=255)
                                            cAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtt                                                                                 >  2:226379/1‑141 (MQ=255)
                                              aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                                                               <  1:210620/141‑1 (MQ=255)
                                                     cGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCACCAGTATTAGCTGCTTTTTGAATTa                                                                        >  2:266215/1‑141 (MQ=255)
                                                      ggatgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                       <  1:46728/139‑1 (MQ=255)
                                                         ttatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTGTCTAATACTGCTCTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACATTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGc                                                                    <  2:224965/139‑1 (MQ=255)
                                                          gatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCa                                                                   <  1:71946/141‑1 (MQ=255)
                                                             ggTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAtt                                                                >  2:9410/1‑141 (MQ=255)
                                                             ggTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAtt                                                                >  1:114309/1‑141 (MQ=255)
                                                             ggTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGACTTTAATTGTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAtt                                                                >  1:66398/1‑141 (MQ=255)
                                                                    ccgCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAg                                                         <  1:337207/139‑1 (MQ=255)
                                                                            aaTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTg                                                 >  1:110633/1‑141 (MQ=255)
                                                                                         tGCTTTCTAATATTGCCTTTAATTCTTTTTTTTTTTTTTTTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGTATTAAATCAATTTCGTCAGCTTCTTTGATGTCTCTAAttt                                    >  2:201666/1‑141 (MQ=255)
                                                                                                         cTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATAAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATATTTGACGTCTCTAATTTTATTTACAGTATTAGa                    >  1:159062/1‑141 (MQ=255)
                                                                                                              aTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGc               >  2:265372/1‑141 (MQ=255)
                                                                                                                 cTTTGTCAGTACTACCTCCTTTTCCACCAGTTAAAATATATTCGTATGTTTTATCAACAATTTTAGCTGATTCTTGTATTAAAGCAATTTCGTCATGACGTTGGACGTCTCTAATTTTATCTACAGTCTTATAAAGGCtta            >  2:104325/1‑141 (MQ=255)
                                                                                                                   ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatt          >  1:299761/1‑141 (MQ=255)
                                                                                                                    tttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatta         <  2:66398/141‑1 (MQ=255)
                                                                                                                    tttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatta         <  2:114309/141‑1 (MQ=255)
                                                                                                                      tCAGTCATTCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAAt       <  2:63357/141‑1 (MQ=255)
                                                                                                                           cATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGatat  >  1:98031/1‑141 (MQ=255)
                                                                                                                                  |                                                                                                                                     
ACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATAT  >  CP000730/1646059‑1646322

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