Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,030,815 G→T A90S (GCA→TCA)  USA300HOU_1885 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,030,8150GT96.8% 219.1 / 0.9 63A90S (GCA→TCA) USA300HOU_1885hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (29/32);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CAATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGCA  >  CP000730/2030681‑2030952
                                                                                                                                       |                                                                                                                                         
caATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGt                                                                                                                                     >  1:614259/1‑141 (MQ=255)
  aTTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTc                                                                                                                                   >  2:196771/1‑141 (MQ=255)
   ttGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCa                                                                                                                                  >  1:97898/1‑141 (MQ=255)
      aaCCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGa                                                                                                                               <  1:116670/141‑1 (MQ=255)
        ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             <  1:483673/141‑1 (MQ=255)
        ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             <  1:438726/141‑1 (MQ=255)
        ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             >  2:551511/1‑141 (MQ=255)
        ccAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGAtt                                                                                                                             >  1:485291/1‑141 (MQ=255)
           aCTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTa                                                                                                                          >  1:137392/1‑141 (MQ=255)
              gCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTc                                                                                                                       <  1:266970/141‑1 (MQ=255)
                    gCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACa                                                                                                                 >  1:513598/1‑141 (MQ=255)
                      aTCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACaaa                                                                                                               <  2:631015/141‑1 (MQ=255)
                        cACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAAc                                                                                                             >  1:5/1‑141 (MQ=255)
                           ttattaAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTc                                                                                                           <  1:451568/141‑1 (MQ=255)
                               taAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAAtt                                                                                                      <  2:338027/141‑1 (MQ=255)
                                aaaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGa                                                                                                                               <  2:613098/115‑1 (MQ=255)
                                aaaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGa                                                                                                                               >  1:613098/1‑115 (MQ=255)
                                aaaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTg                                                                                                     >  1:462976/1‑141 (MQ=255)
                                 aaGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGc                                                                                                    <  2:429991/141‑1 (MQ=255)
                                         cGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGt                                                                                            >  1:77703/1‑141 (MQ=255)
                                         cGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGt                                                                                            >  1:387145/1‑141 (MQ=255)
                                            aGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACg                                                                                         <  1:445498/141‑1 (MQ=255)
                                             gATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGa                                                                                        <  2:177438/141‑1 (MQ=255)
                                             gATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGa                                                                                        >  1:125047/1‑141 (MQ=255)
                                                    tGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGc                                                                                 >  1:449007/1‑141 (MQ=255)
                                                    tGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTGATTCTGCGATCAAGCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGc                                                                                 <  1:383833/141‑1 (MQ=255)
                                                     gCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCt                                                                                <  2:338057/141‑1 (MQ=255)
                                                     gCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCt                                                                                >  2:312930/1‑141 (MQ=255)
                                                         cGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTAc                                                                            <  1:278184/141‑1 (MQ=255)
                                                            tAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTaa                                                                         <  2:133876/141‑1 (MQ=255)
                                                              aTTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAca                                                                       >  2:241280/1‑141 (MQ=255)
                                                               ttGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAcat                                                                      <  1:96232/141‑1 (MQ=255)
                                                                  gTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAAcatcat                                                                   <  1:287274/141‑1 (MQ=255)
                                                                         cTTTGCAGTATTATTTACATATTTATTCGGCGAACAATCACCATTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTAt                                                            <  1:125432/141‑1 (MQ=255)
                                                                            ttGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGa                                                          >  2:402826/1‑141 (MQ=255)
                                                                               cAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTc                                                       <  1:141990/141‑1 (MQ=255)
                                                                               cAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTc                                                       <  2:355769/141‑1 (MQ=255)
                                                                                 gTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCa                                                     >  2:519088/1‑141 (MQ=255)
                                                                                      attTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAGCATCATTAGCTATGATTCCAGGTAt                                                <  1:59634/141‑1 (MQ=255)
                                                                                       tttACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTAtt                                               >  1:344464/1‑141 (MQ=255)
                                                                                          aCATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTCAAATCATTAGCTATGATTCCAGGTATTCat                                            <  2:485291/141‑1 (MQ=255)
                                                                                           cATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCatg                                           >  2:502817/1‑141 (MQ=255)
                                                                                                ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTTTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      >  1:579734/1‑141 (MQ=255)
                                                                                                ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      >  2:307917/1‑141 (MQ=255)
                                                                                                ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      <  2:467794/141‑1 (MQ=255)
                                                                                                ttattCGGCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGcc                                      >  1:48789/1‑141 (MQ=255)
                                                                                                       gCGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACatat                               >  2:462245/1‑141 (MQ=255)
                                                                                                        cGATCAATCACCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATAtt                              <  1:402826/141‑1 (MQ=255)
                                                                                                          atcaatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCa                            >  1:322180/1‑141 (MQ=255)
                                                                                                           tcaatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCaa                           >  2:281894/1‑141 (MQ=255)
                                                                                                            caatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAAc                          >  1:516926/1‑141 (MQ=255)
                                                                                                             aatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACt                         <  1:295314/141‑1 (MQ=255)
                                                                                                             aatcaCCTTTTACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACt                         <  2:5/141‑1 (MQ=255)
                                                                                                                    ttttACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGtc                  <  1:241280/141‑1 (MQ=255)
                                                                                                                     tttACATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGtct                 <  2:619413/141‑1 (MQ=255)
                                                                                                                        aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              <  1:298260/141‑1 (MQ=255)
                                                                                                                        aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              >  2:169359/1‑141 (MQ=255)
                                                                                                                        aCATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGa              <  2:125047/141‑1 (MQ=255)
                                                                                                                         cATATGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCGGTACTAACATCATTAGCTATGATACCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGaa             >  2:383523/1‑141 (MQ=255)
                                                                                                                            atGCATTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACAt          <  1:635465/141‑1 (MQ=255)
                                                                                                                                aTTGAGTTCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGAc      <  1:312930/141‑1 (MQ=255)
                                                                                                                                   gAGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACTTATTCAACTTTTTGTCTCGAACATTGACAGc   <  2:282661/141‑1 (MQ=255)
                                                                                                                                    aGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGCa  <  2:273075/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CAATTGAACCAACTGCCAAGGCATCACTTATTAAAGGTTATCGTAGATTACCTGCTACGGTAATTGGTGCAGGCTTTGCAGTATTATTTACATATTTATTCGGCGATCAATCACCTTTTACATATGCATTGAGTGCAACGTTCACGATTTTATTCTGTACAAAACTCAAATTGCAAGTTGGTACGAATGTCGCTGTACTAACATCATTAGCTATGATTCCAGGTATTCATGATGCCTACATATTCAACTTTTTGTCTCGAACATTGACAGCA  >  CP000730/2030681‑2030952

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