Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,842,707 T→C 6.9% I204V (ATT→GTT)  tyrS ← Tyrosine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,842,7070TC 6.9% 184.8 / 5.0 58I204V (ATT→GTT) tyrSTyrosine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base C (2/2);  total (33/25)
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.60e-01

CTTGTTTTTTCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCTAAACGAC  >  USA300TCH1516_ALE/1842574‑1842847
                                                                                                                                     |                                                                                                                                            
cTTGTTTTTTCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACtt                                                                                                                                       >  2:346817/1‑141 (MQ=255)
    ttttttCAGCATCTAACCAAACAGCACCTGACTCAGTCTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGa                                                                                                                                   >  2:154677/1‑141 (MQ=255)
    ttttttCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGa                                                                                                                                   <  2:219422/141‑1 (MQ=255)
    ttttttCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGa                                                                                                                                   >  1:271446/1‑141 (MQ=255)
    ttttttCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGa                                                                                                                                   <  1:266136/141‑1 (MQ=255)
    ttttttCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGa                                                                                                                                   <  2:390742/141‑1 (MQ=255)
          cAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTAc                                                                                                                             >  1:102059/1‑141 (MQ=255)
          cAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAACACCACTTGTGATATTAc                                                                                                                             >  1:226280/1‑141 (MQ=255)
          cAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAACACCACTTGTGATATTAc                                                                                                                             >  1:226287/1‑141 (MQ=255)
              aTCTAACCAACCAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCa                                                                                                                         <  2:414646/141‑1 (MQ=255)
              aTCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCa                                                                                                                         <  2:279984/141‑1 (MQ=255)
                cTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCAtt                                                                                                                       >  2:66761/1‑141 (MQ=255)
                             aCCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCAcc                                                                                                          >  1:168880/1‑141 (MQ=255)
                                 gACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACt                                                                                                      <  1:452589/141‑1 (MQ=255)
                                   cTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTg                                                                                                    <  2:316387/141‑1 (MQ=255)
                                       gACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATc                                                                                                >  2:403059/1‑141 (MQ=255)
                                       gACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCATATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATc                                                                                                >  1:234968/1‑141 (MQ=255)
                                                       tACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATtctc                                                                                >  2:227661/1‑141 (MQ=255)
                                                       tACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATtctc                                                                                >  1:312490/1‑141 (MQ=255)
                                                        aCCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTctct                                                                               >  1:429856/1‑141 (MQ=255)
                                                        aCCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTctct                                                                               >  2:349051/1‑141 (MQ=255)
                                                        aCCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTctct                                                                               >  2:315648/1‑141 (MQ=255)
                                                              tGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCaa                                                                         <  2:141890/141‑1 (MQ=255)
                                                               gATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCaaa                                                                        >  1:121711/1‑141 (MQ=255)
                                                                 tttAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATg                                                                      >  2:427856/1‑141 (MQ=255)
                                                                  ttAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGa                                                                     <  2:323718/141‑1 (MQ=255)
                                                                  ttAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGa                                                                     <  2:265531/141‑1 (MQ=255)
                                                                                   aGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTg                                                    >  1:318538/1‑141 (MQ=255)
                                                                                   aGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTg                                                    >  1:42601/1‑141 (MQ=255)
                                                                                   aGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTg                                                    >  2:285455/1‑141 (MQ=255)
                                                                                   aGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTg                                                    >  2:28515/1‑141 (MQ=255)
                                                                                     ttAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTa                                                  >  2:316067/1‑141 (MQ=255)
                                                                                     ttAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTa                                                  >  2:436838/1‑141 (MQ=255)
                                                                                     ttAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTa                                                  >  1:117701/1‑141 (MQ=255)
                                                                                     ttAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTa                                                  >  1:353350/1‑141 (MQ=255)
                                                                                     ttAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTa                                                  >  2:141225/1‑141 (MQ=255)
                                                                                     gtAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTa                                                  >  1:436454/2‑141 (MQ=255)
                                                                                         aCCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAAt                                              <  1:154677/141‑1 (MQ=255)
                                                                                         aCCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAAt                                              <  2:271446/141‑1 (MQ=255)
                                                                                           cGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCg                                            >  2:408101/1‑141 (MQ=255)
                                                                                             tATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCgtg                                          <  1:190330/141‑1 (MQ=255)
                                                                                             tATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCgtg                                          <  1:201456/141‑1 (MQ=255)
                                                                                                 cGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTAtg                                      <  1:7673/141‑1 (MQ=255)
                                                                                                 cGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACGTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATGCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTAtg                                      <  1:193952/141‑1 (MQ=255)
                                                                                                  gTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTAtgt                                     <  1:28296/141‑1 (MQ=255)
                                                                                                   tCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTAtgtg                                    >  2:324510/1‑141 (MQ=255)
                                                                                                     tGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGaa                                  >  1:458120/1‑141 (MQ=255)
                                                                                                     tGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGaa                                  >  2:34955/1‑141 (MQ=255)
                                                                                                       tttGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTCACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAAtt                                >  2:424850/1‑141 (MQ=255)
                                                                                                                     cgacgCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATGAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATa                  <  1:249686/141‑1 (MQ=255)
                                                                                                                     cgacgCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATa                  <  1:211410/141‑1 (MQ=255)
                                                                                                                       acgCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATAcc                <  1:346817/141‑1 (MQ=255)
                                                                                                                       acgCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATAcc                <  2:194484/141‑1 (MQ=255)
                                                                                                                           aTTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGt            >  1:319998/1‑141 (MQ=255)
                                                                                                                              aaTTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCt         <  2:68880/141‑1 (MQ=255)
                                                                                                                              aaTTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCt         <  2:231975/141‑1 (MQ=255)
                                                                                                                                ttCAACACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCTaa       <  2:226280/141‑1 (MQ=255)
                                                                                                                                ttCAACACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCTaa       <  2:226287/141‑1 (MQ=255)
                                                                                                                                     tACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCTAAACGAc  >  2:210107/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                            
CTTGTTTTTTCAGCATCTAACCAAACAGCACCTGACTCAGACTTACCAAATTTCTTACCATCTGATTTAGTTACAAGCGGAATAGTTAAACCGTATGCGTCTGTTTGACCATACATACGACGCATTAATTCAATACCACTTGTGATATTACCCCATTGATCTGATCCACCTACTTGAATCTTACAATTCAATTCTCTATTCAAATGACCGAAATCAATAGCTTGTAAAATCGTGTATGTGAATTCTGTATATGAAATACCATGTTCTAAACGAC  >  USA300TCH1516_ALE/1842574‑1842847

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