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,1890GC97.2% 139.1 / ‑5.3 45A155G (GCA→GGA) USA300HOU_1531M24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (24/18);  minor base T (1/1);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACG  >  CP000730/1646054‑1646325
                                                                                                                                       |                                                                                                                                        
ttATCACTGGCCACACCATGTGGTAATGCACCTCTATGACCAGATGCTCAAATCGTATCTAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCATTCATGCCTCCtttt                                                                                                                                     <  1:293387/141‑1 (MQ=255)
       ttGCAACACCATGTGGTAATGCACCTCTATGACCAGATTCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTATTTTTataacat                                                                                                                              >  2:465673/1‑140 (MQ=255)
       ttGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                                                                              >  1:201600/1‑141 (MQ=255)
            acacCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCAATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaa                                                                                                                         >  2:246996/1‑141 (MQ=255)
               ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                      <  1:553786/141‑1 (MQ=255)
               ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                      <  1:438036/141‑1 (MQ=255)
                cATGTGGTAATGCACCCCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTACTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTCAAAattt                                                                                                                     >  2:32258/1‑139 (MQ=255)
                          tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACGGTTAAAATATATTCATATGt                                                                                                           >  1:5106/1‑141 (MQ=255)
                              ccTCTATGACCAGATGCTACAATCGTATCGAACGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtca                                                                                                       <  1:106754/141‑1 (MQ=255)
                                tctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCATTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                                                                     <  1:375289/141‑1 (MQ=255)
                                       ccAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCACCTTTTCCAACAGTTAAACTATATTCATATGTTTCAGTGACAata                                                                                              >  1:138098/1‑141 (MQ=255)
                                            tGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                                                                         <  1:278638/141‑1 (MQ=255)
                                            tGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgc                                                                                         <  1:278829/141‑1 (MQ=255)
                                               tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                                      >  2:563184/1‑141 (MQ=255)
                                                   aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                                                                  <  1:97683/141‑1 (MQ=255)
                                                    tCGTATCAAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCttttt                                                                                 <  2:204557/141‑1 (MQ=255)
                                                      gTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTAATCAACAATATTAGCTGCCTTTTGa                                                                               >  2:400147/1‑141 (MQ=255)
                                                        aTCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAACATTAGCTGCTTTTTGAAt                                                                             >  1:251450/1‑141 (MQ=255)
                                                             atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAg                                                                        <  1:454811/141‑1 (MQ=255)
                                                                     ccATCTGCTCATAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCCACACGATTACCTGCTTTTTGAATTAAAGCAATGTCg                                                                >  1:406769/1‑141 (MQ=255)
                                                                      cATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGt                                                               >  1:421035/1‑141 (MQ=255)
                                                                      cATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGAAATTTCGt                                                               >  1:245030/1‑141 (MQ=255)
                                                                            atCCTAATTCTAGAATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTCAAATATATTCATATGTTTCATCAACAATATTCACTGCTTTTTGAATTAAAGCAATATTGTCAGCAt                                                         >  2:516019/2‑141 (MQ=255)
                                                                                tAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTTTTTTTTTTTTAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAATCATTTTCTTCAGCTTCttt                                                     >  2:282632/1‑141 (MQ=255)
                                                                                   ttCTAGCATCTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCCTCAGCATCTTTGAc                                                  <  1:81509/141‑1 (MQ=255)
                                                                                     cttgCATTTTCCTTTCTAATAGTGCCTTTAATTCTTTTTCACTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATGAGCTGCGTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGt                                                <  2:530555/138‑1 (MQ=255)
                                                                                     cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTTTTTTTTTTTTTTTTTTTTTTCTTTTAAAATATATTTATATTTTTCATCAACAATATTAGCTGCTTTTTGAATTTAAGCAATTTCGTCAGCATCTTTTActt                                                >  2:125210/1‑139 (MQ=255)
                                                                                      tAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtc                                               <  1:563184/141‑1 (MQ=255)
                                                                                        gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtctc                                             >  1:122106/1‑141 (MQ=255)
                                                                                         cATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                            >  1:399163/1‑141 (MQ=255)
                                                                                             ttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                                        >  1:118895/1‑141 (MQ=255)
                                                                                               gCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtttt                                      >  2:511516/1‑141 (MQ=255)
                                                                                                 tttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACCGTTCACAAATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTAt                                    <  1:128925/141‑1 (MQ=255)
                                                                                                           tGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAtt                          <  2:5106/141‑1 (MQ=255)
                                                                                                            gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTa                         <  1:295154/141‑1 (MQ=255)
                                                                                                               tttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTAAAGTATTAGaa                      >  2:437589/1‑141 (MQ=255)
                                                                                                                      cTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCtta               >  2:484455/1‑141 (MQ=255)
                                                                                                                        ttttCGGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatt             >  2:322811/1‑141 (MQ=255)
                                                                                                                         tttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatta            <  1:246996/141‑1 (MQ=255)
                                                                                                                           tCAGTGATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACACTATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAAt          >  2:288096/1‑141 (MQ=255)
                                                                                                                            cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATg         >  1:329062/1‑141 (MQ=255)
                                                                                                                             aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGa        >  1:551302/1‑141 (MQ=255)
                                                                                                                                cATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGatat     >  2:46003/1‑141 (MQ=255)
                                                                                                                                  tGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATAc   <  2:551302/141‑1 (MQ=255)
                                                                                                                                   gCCTTCTTCTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACg  <  1:222972/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACG  >  CP000730/1646054‑1646325

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: ‑