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

TTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGAT  >  CP000730/1646053‑1646320
                                                                                                                                        |                                                                                                                                   
tttATCACTTGCAACTCCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATTGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTATTATTTCCTTTAATTCTTTTTCAGTCATCCCCTCttt                                                                                                                                 >  2:791407/1‑141 (MQ=255)
tttATCACCTGCACCACCATGTGGTAATGCACCCCTATGACCAGATGCTACAATCGCATCGAATGATGGGCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                                                                 <  2:711346/141‑1 (MQ=255)
 ttATCACTTGGAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCGCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtttt                                                                                                                                <  2:263008/141‑1 (MQ=255)
 ttATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtttt                                                                                                                                <  1:268675/141‑1 (MQ=255)
   aTCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCTGTCATGCCTCCTTTTac                                                                                                                              >  1:824858/1‑141 (MQ=255)
    tCACTTGCAACACCATGTGGTAATGCACCTCTATGACCCGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTaca                                                                                                                             >  1:526280/1‑141 (MQ=255)
     cACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa                                                                                                                            >  1:629220/1‑141 (MQ=255)
           cAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTCGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTTAGTCATGCCTCCTTTTACAACAGTTa                                                                                                                      <  2:564679/141‑1 (MQ=255)
           cAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTCCCGTCCGGCCTCCTTTTACAACAGTTa                                                                                                                      <  1:226701/141‑1 (MQ=255)
               acCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACGACAGTTAAAat                                                                                                                  >  1:340026/1‑141 (MQ=255)
               acCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAat                                                                                                                  >  1:14319/1‑141 (MQ=255)
                 cATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCTGTTTTTCTTCCTTTTACAACAGTTAAAatat                                                                                                                >  2:387124/1‑141 (MQ=255)
                 cATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCTGTCATGCCTCATTTTACCACAGTTacatcgt                                                                                                                >  2:738840/1‑135 (MQ=255)
                   tgtgGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTGTTTCAGTCATGCCTCCTTTTACAACGGTTAAAatatat                                                                                                              >  2:871362/1‑141 (MQ=255)
                           tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                      <  1:579206/141‑1 (MQ=255)
                           tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                      <  1:728221/141‑1 (MQ=255)
                           tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                      <  2:247380/141‑1 (MQ=255)
                                 ttcaTGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                                                                <  2:257555/138‑1 (MQ=255)
                                 tctATAACCGGTTGCCACAATCCTGTCGAATGATTGTCCACTTGCTCCTAATTCTATCATTTTGCTTTTTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcagc                                                                                                <  2:175695/141‑3 (MQ=255)
                                 cctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                                                                <  1:381228/140‑1 (MQ=255)
                                  ctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatca                                                                                               <  2:750756/141‑1 (MQ=255)
                                      gACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaacaa                                                                                           >  1:635540/1‑141 (MQ=255)
                                                tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                                 <  1:495367/141‑1 (MQ=255)
                                                tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                                 <  2:478916/141‑1 (MQ=255)
                                                tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                                 <  1:349962/141‑1 (MQ=255)
                                                   aaTCGTATTGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCttt                                                                              <  1:169092/141‑1 (MQ=255)
                                                    aTCGTATCGAATGATGGTCCATCTGCTCCTTATTTTAGCATTTTGCTTTCTAATCTTGCCTTTAATTCTTTTTCAGTCATGCCTCCTATTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                                                             <  1:546379/141‑1 (MQ=255)
                                                    aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                                                             >  2:930870/1‑141 (MQ=255)
                                                          tCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAtt                                                                       <  1:732844/141‑1 (MQ=255)
                                                               tgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGc                                                                  <  2:252614/141‑1 (MQ=255)
                                                               tgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAACATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGc                                                                  <  2:3833/141‑1 (MQ=255)
                                                                 atgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCaa                                                                >  2:65669/1‑141 (MQ=255)
                                                                        ctCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTc                                                         <  1:387124/140‑1 (MQ=255)
                                                                        aTCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTc                                                         <  1:738840/141‑1 (MQ=255)
                                                                          cTGCTTCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAg                                                       <  1:342500/141‑1 (MQ=255)
                                                                           tGCTCCTAATTTTAGCATGTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCGTCTTTTTGACTTAAAGCAATTTCGTCAGc                                                      <  2:391889/141‑1 (MQ=255)
                                                                           tGCTCCTAAATCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCTTTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc                                                      <  1:358286/141‑1 (MQ=255)
                                                                             cTCCTAATTCTAGCATGTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCAt                                                    >  1:526416/1‑141 (MQ=255)
                                                                               ccTAATTCTAGTATTTTGCTTTCTAATATCAACTTATAGAATTTTCTTTTCTTGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAACGACATTTCGCCAGCAtat                                                  <  2:902828/141‑3 (MQ=255)
                                                                                     tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACg                                            >  1:877937/1‑141 (MQ=255)
                                                                                          cATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                       <  2:842621/141‑1 (MQ=255)
                                                                                            ttttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTaa                                     <  1:476164/141‑1 (MQ=255)
                                                                                                    tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACCATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCt                             >  1:635851/1‑141 (MQ=255)
                                                                                                    tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGCATTAAAGCAATTTCGTATTCATCTTTGANGTCTCTAATTTTATCt                             >  1:726455/1‑141 (MQ=255)
                                                                                                      tAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTAc                           >  2:521973/1‑141 (MQ=255)
                                                                                                          atTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTa                       <  2:823966/141‑1 (MQ=255)
                                                                                                          atTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTa                       >  1:86468/1‑141 (MQ=255)
                                                                                                          atTGCCTTTAATTCTTTTTCAGTCATGCCTCCTCTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTCCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTa                       <  2:526280/141‑1 (MQ=255)
                                                                                                               cTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGa                  <  1:37831/141‑1 (MQ=255)
                                                                                                                tttAATTCTTTTTCAGTCATGCCTCCTTTTACATCAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGaa                 <  1:405271/141‑1 (MQ=255)
                                                                                                                tttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAACTGCTCATTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACATTATTCGaa                 >  2:335216/1‑141 (MQ=255)
                                                                                                                 ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGaaa                >  1:900668/1‑141 (MQ=255)
                                                                                                                 ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCAACAGTATTAGaaa                >  1:131715/1‑141 (MQ=255)
                                                                                                                    aTTCTGTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCCTATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGc             <  2:365745/141‑1 (MQ=255)
                                                                                                                        tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttat         <  2:735368/141‑1 (MQ=255)
                                                                                                                         ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACACTATTTGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTCCAGTATTAGCAATGcttctt        >  1:129705/1‑141 (MQ=255)
                                                                                                                              aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATg  >  2:801388/1‑140 (MQ=255)
                                                                                                                               gTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGat  >  2:587789/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGAT  >  CP000730/1646053‑1646320

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: