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

CCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATA  >  CP000730/1646069‑1646323
                                                                                                                        |                                                                                                                                      
ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                    >  1:197456/1‑141 (MQ=255)
     tgGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTc                                                                                                               <  1:156081/141‑1 (MQ=255)
                    aTGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaac                                                                                                <  1:141121/141‑1 (MQ=255)
                     tGACCAGATGCTACACTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATGTGTTTCATcaaca                                                                                               >  1:124917/1‑141 (MQ=255)
                             tGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACGACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                                                                       <  2:225990/141‑1 (MQ=255)
                             tGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                                                                       <  1:75015/141‑1 (MQ=255)
                                              atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATCAGCTGCTTTTTGAATTAAAg                                                                      >  2:43259/1‑141 (MQ=255)
                                                gatgGTGCATCTGGTCCTAATTCTAGCATTTTGCTTTCTAATATGGCCGTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACGGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCa                                                                    <  2:7869/141‑1 (MQ=255)
                                                 atgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCaa                                                                   <  2:197456/141‑1 (MQ=255)
                                                           tGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc                                                         <  2:217229/141‑1 (MQ=255)
                                                               ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATGGCCAATATTAGCGGCTTGTTGAATTAAAGCAATTTCGTCAGCATCt                                                     <  2:172479/141‑1 (MQ=255)
                                                               ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt                                                     <  1:224731/141‑1 (MQ=255)
                                                                         gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtctc                                           <  2:11458/141‑1 (MQ=255)
                                                                            ttttGCTTTCTATTATTGTCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTaa                                        <  1:126882/141‑1 (MQ=255)
                                                                              ttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTCTATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                                      <  2:164200/141‑1 (MQ=255)
                                                                               tGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAttt                                     <  1:165732/141‑1 (MQ=255)
                                                                                    tCTAATGTTGCCTATCATTCTTGTTCTGTCGTGCCTCCTTTTACATCAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCCATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCt                                <  2:186213/141‑1 (MQ=255)
                                                                                                 ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCCGCATCTTTGACGTCTCTCATTTTCTCTCCAGTATTAGaac                   >  1:75399/1‑140 (MQ=255)
                                                                                                 ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTAGTCAGCATCTTTGACGTCGCTAATTTTATCTACAGTATAAGaaa                   >  1:51739/1‑141 (MQ=255)
                                                                                                    aTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGACTCTAATTTTATCTACAGTATTAGAAATGc                >  1:137951/1‑141 (MQ=255)
                                                                                                                  aTGCCTCCTTCTACCACATTTAAAACATATTCAGATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATCTCCTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAAAAATGCTTATTAATGAtata  <  2:39066/141‑1 (MQ=255)
                                                                                                                        |                                                                                                                                      
CCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATA  >  CP000730/1646069‑1646323

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: